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THERMO-MECHANICAL  AND  THERMAL 
BEHAVIOR  OF  HIGH-TEMPERATURE  STRUCTURAL 
MATERIALS 


PREFACE 


Technical  ceramics  because  of  their  chemical  inertness,  high  melting 
point,  good  wear  resistance,  excellent  mechanical  stability  at  high  tempera¬ 
ture  and  other  unique  properties,  represent  a  class  of  materials  eminently 
suited  for  many  critical  engineering  applications.  Unfortunately,  because 
of  their  brittleness  and  unfavorable  combination  of  pertinent  material  pro¬ 
perties,  technical  ceramics  generally  are  highly  susceptible  to  catastrophic 
failure  in  non-uniform  thermal  environments,  which  give  rise  to  thermal 
stresses  of  high  magnitude. 

Thermal  stress  failure  analysis  of  structural  materials  represents  a 
multi-disciplinary  problem  which  involves  the  principles  of  heat  transfer, 
mechanics  and  materials  engineering.  Over  the  last  few  decades  much  general 
understanding  of  the  nature  of  thermal  stress  failure  of  brittle  materials 
has  been  generated.  However,  due  to  the  multi-disciplinary  nature  of  the 
problem,  the  ability  to  predict  thermal  stress  failure  quantitatively  for 
design  or  other  purposes  has  lagged  behind  the  progress  made  in  other  engi¬ 
neering  fields.  The  objective  of  the  present  program  is  to  improve  the 
qualitative  and  quantitative  understanding  of  the  nature  of  thermal  stress 
failure  of  brittle  structural  materials.  In  order  to  achieve  this  objective, 
the  participating  investigators  and  scope  of  the  program  are  organized 
such  that  full  advantage  is  taken  of  the  combined  inputs  from  a  number  of 
engineering  disciplines.  This  report,  in  the  form  of  individual  chapters, 
presents  studies  completed  for  publication  during  the  period  1/1/79  to 
12/31/79.  For  one  study,  because  of  its  length  and  to  keep  the  printing 
cost  to  a  minimum,  only  the  abstract  is  given.  Upon  publication,  reprints 
can  be  sent  to  the  readers  interested  in  further  detail.  Together  with  the 
eleven  papers  reported  previously,  the  papers  contained  in  the  present 
report  bring  the  total  number  of  completed  studies  to  twenty-three  for  the 
contract  period  to  date  (A/1/78  -  12/31/79). 


The  readers  of  this  report  are  invited  to  provide  their  comments. 
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CHAPTER  I 


ANALYSIS  OF  EFFECT  OF  HEAT  TRANSFER  VARIABLES  ON  THERMAL  STRESS  RESISTANCE 
OF  BRITTLE  CERAMICS  MEASURED  BY  QUENCHING  EXPERIMENTS 


by 


J.  P.  Singh,  J.  R.  Thomas  and  D.  P.  H.  Hasselman 


Departments  of 

Materials  and  Mechanical  Engineering 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061 


ABSTRACT 

Heat  transfer  theory  was  used  to  analyze  the  variables  which  con¬ 
trol  the  thermal  stress  fracture  of  brittle  ceramic  specimen  subjected 
to  rapid  cooling  by  quenching  into  fluid  media.  Expressions  weje  de¬ 
rived  for  the  maximum  quenching  temperature  difference  to  which  circular 
solid  cylindrical  specimens  can  be  subjected  in  terms  of  the  pertinent 
properties  of  the  material  and  fluid  media.  Appropriate  thermal  stress 
resistance  parameters  were  derived.  The  specimen  density  was  introduced 
as  an  additional  property  which  controls  thermal  stress  resistance.  Re¬ 
commendations  were  made  for  procedures  to  be  followed  in  conducting 
quenching  experiments. 
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I.  Introduction 

Technical  ceramics  due  to  their  unique  physical  properties,  are 
excellent  candidate  materials  for  engineering  applications  involving 
extreme  thermal  and  chemical  environments.  Unfortunately,  due  to  their 
brittle  nature,  ceramic  materials  are  highly  susceptible  to  catastrophic 
failure  due  to  high  thermal  stresses  encountered  during  transient  or 
steady  state  heat  transfer.  For  this  reason,  proper  selection  of 
materials  for  high  temperature  application  is  critical  to  maintain 
structural  integrity.  Such  selection  can  be  based  on  theoretical 
analysis1*2  as  well  as  on  experimental  testing3"7.  A  popular  test 
for  the  thermal  stress  resistance  of  brittle  ceramics  consists  of 
quenching  specimens  from  higher  temperatures  in  a  fluid  bath  at  lower 
temperature.  Such  fluids  may  include  water,  oil,  eutectic  mixtures 
of  salt  and  fluidized  beds.  These  tests  are  useful  for  assessing  re¬ 
lative  difference  in  thermal  stress  resistance  of  different  materials. 
Unfortunately,  however,  at  this  time  due  to  the  complex  nature  of 
heat  transfer,  quenching  experiments  do  not  appear  to  be  very  suitable 
for  a  quantitative  assessment  of  thermal  stress  resistance.  It  is 
the  purpose  of  this  paper  to  discuss  the  variables  which  affect  heat 
transfer  in  quenching  tests,  and  to  derive  appropriate  expression 
for  the  heat  transfer  coefficient  in  order  to  provide  an  improved 
quantitative  basis  for  the  interpretation  of  the  results  for  thermal 
stress  resistance  determined  by  quenching  experiments. 

II.  Analysis 

A.  General . 

In  general,  quenching  experiments  can  be  performed  in  two  different 
ways,  each  corresponding  to  a  different  mode  of  heat  transfer.  In  one 


3 


type  of  experiment3-5,  the  specimen  is  preheated  to  thermal  equilibrium 
at  a  specified  temperature  and  then  quickly  transferred  to  the  fluid 
bath  where  it  is  held  stationary.  In  this  situation  the  heat  is  trans¬ 
ferred  from  the  specimen  surface  by  natural  convection  due  to  the  fluid 
motion  resulting  from  the  density  change  caused  by  the  heating  process®. 
In  the  second  type  of  experiments,  the  preheated  specimen  is  dropped® 
into  and  falls  through  the  stationary  fluid.  In  this  case,  there  is  a 
relative  motion  between  the  fluid  and  the  specimen  and  the  heat  trans¬ 
fer  takes  place  by  forced  convection  as  a  result  of  the  imposed  fluid 
flow  by  a  stream  around  the  heated  object®.  General  equations  for 
heat  transfer  will  be  obtained  for  both  modes  namely  'natural'  and 
'forced'  convection.  For  convenience  of  mathematical  analysis,  the 
simple  specimen  geometry  of  a  circular  horizontal  cylinder  will  be 
considered.  However,  the  analysis  can  easily  be  extended  to  different 
geometries.  Also  the  analysis  will  be  limited  to  the  conditions  where 
nucleate  boiling  and  film  formation  do  not  take  place. 

B.  Natural  Convection  (stationary  specimen). 

The  average  heat  transfer  coefficient  h  from  an  immersed  horizontal 
cylinder  in  a  stationary  fluid  is  given  by®: 

h  '  c<NcVNPr>"  T  (1) 

where  n  and  C  are  constants  which  depend  on  the  value  of  (N^*  Np^) >  d 
is  the  diameter  of  the  cylinder  and  subscript  ' f ’  refers  to  the  fluid 
property  values  at  the  film  temperature,  generally  taken  as  the  mean 
of  the  specimen  surface  and  the  fluid  temperature. 


Using  the  definition  of  NCr  and  Npr,  the  value  of  h  can  be  ob¬ 
tained  from  equation  1  in  terms  of  material  and  environment  properties. 
Thus, 

h  =  C(gypfCpf/vf)n(AT)nd3n_1Kf1~n  (2) 

C.  Forced  Convection  (falling  specimen). 

In  addition  to  many  other  parameters,  the  heat  transfer  by  forced 
convection  also  depends  upon  the  velocity  with  which  the  specimen  falls 
through  the  fluid  in  a  quenching  experiment.  It  is  important,  therefore, 
to  evaluate  the  specimen  velocity  before  the  general  expression  for  h 
can  be  obtained.  In  a  quenching  experiment,  while  the  specimen  is  fal¬ 
ling  through  the  fluid  bath,  the  downward  force  due  to  the  specimen 

•k 

weight  is  opposed  by  the  upward  buoyancy  and  shear  force  .  When  the 
upward  forces  balance  the  weight  of  the  specimen,  the  specimen  attains 
an  equilibrium  constant  velocity  frequently  termed  as  terminal  velocity10. 
For  mathematical  simplicity  it  can  be  assumed  that  the  specimen  attains 
the  terminal  velocity  as  soon  as  it  hits  the  fluid  surface. 

For  a  horizontal  cylinder  with  its  axis  perpendicular  to  flow, 
at  terminal  velocity  V  , 


TT  ,2.  7f  2,  ,  _ 

7  d  Lp  =  7  d  Lp.  +  F 
4  s  4  f 


(3) 


where  F  is  given  by11  f^p^V^  Ap/2g,  friction  coefficient  f  is  a  func¬ 
tion  of  Reynold's  number  and  Ap  is  the  projected  area.  Using  the  appropriate 
values  for  f^  from  Table  I,  the  values  of  V  are  obtained  as: 


N_  <0.5 
Re 


12  Ps  -1  '*-c* 

V  =  (-  d  g(  ~  -  l)  v  '  (2.002  -  In  —  ) 

tJO  Pj-  t  Vj 


(3a) 


*  Shear  force  is  a  function  of  the  specimen  velocity. 


0.5  <  N  <10 

Kfc  v  n  0.714,1. 143, Ps  ,,0.714  -0.429 

Vt  =  0.267  g  d  (- - 1)  vf  (3b) 


10  <  N  <  40 

„  n  *  0. 588  0.765, s  .,0.588„  -0.176 
V  =  0 . 5g  d  (—  -  1) 


(3c) 


40  <  N  <  400 

Re  n  ,,,,  0.6035 ,0.8105,Ps  ,,0.6035  -0.2070 

V  =  0.4513  g  d  ( —  -1)  v c 

z  p .  r 

(3d) 


4xl02  <  Nn  <  4xl05 
Re 


,  0.5 .0.5  ’  s  ,,0.5 

V  =  1.45  g  d  ( —  -1) 

t  Pf 


(3e) 


These  values  of  the  terminal  velocity  were  used  in  the  expression 
for  NRe  and  the  equations  for  heat  transfer  coefficient  were  obtained 
for  the  range  of  Reynold's  number  commonly  encountered  in  the  thermal 
quenching  experiments. 

The  average  heat  transfer  coefficient  h  in  forced  convection  for 
a  circular  cylinder  in  a  fluid  with  its  axis  perpendicular  to  the  flow 
is  given  by14. 


m  I  K 

h  =  Cr(NRe)m  (Npr)3  ^  (4) 

d 

where  m  and  C  are  constants  whose  values  depend  on  and  N_  .  As 
r  K  Pr  Re 

far  as  the  heat  transfer  is  concerned,  this  case  is  identical  to  one 
in  which  a  horizontal  cylinder  falls  through  the  stationary  fluid,  since 
in  both  cases  the  heat  is  transferred  due  to  the  relative  motion  bet¬ 
ween  the  cylinder  and  the  fluid. 

Substituting  for  N0  and  N_  in  terms  of  velocity  and  fluid  pro- 
Re  Pr 

* 

perties  and  using  equation  (3e)  for  velocity,  h  can  be  expressed  as: 


*  Equation  (3e)  has  been  used  for  the  velocity  to  obtain  an  expression 
for  h  because  this  represents  the  range  of  Reynold's  number  typical  for 
laboratory  thermal  quench  experiments.  However,  other  equations  should 
he  used  whenever  appropriate. 
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400  <  N„  <  4000,  N„  >  0.6 
Re  Pr 

m  1  _  1  1  2  3  m_1 

h  =  Cr  f2.1g(~  -  l)]2  vf3  pf3Cpf3Kf3d2  (5) 

where  m  =  0.466  and  C  =  0.683 

r 

It  should  be  pointed  out  that  the  equations  2  and  5  represent  the 
average  value  of  h  over  the  surface  of  the  specimen.  There  are,  in  fact, 
local  variations  in  h  over  the  specimen  surface.  To  incorporate  this 
effect  in  the  quantitative  evaluation  of  thermal  stress  resistance  is  the 
topic  of  future  research. 

D.  Calculation  of  (AT) 

max 

The  maximum  temperature  differenct  (AT)max>  to  which  the  specimen 
can  be  subjected  without  failure  in  a  thermal  quench  experiment,  can 
be  obtained  by  equating  the  maximum  thermal  stress  to  the  tensile 
strength  St  of  the  specimen.  The  exact  solutions15  of  the  thermal 
stresses  developed  in  a  circular  cylinder  in  case  of  transient  cooling 
are  too  complex  to  yield  any  simple  analytical  solution  of  (AT)max  for 
a  given  heat  transfer  coefficient  value.  However,  when  the  maximum 
tangential  stress  o^15  was  plotted  as  a  function  of  Biot's  number  (figure 
not  shown  here),  it  resulted  in  a  straight  line  represented  by: 

4-  =  1.451  (1  +  ~ — )  (6) 

amax  6 

where  a*  =  (l-v)a_/(aEAT) 
max  y 

For  thermal  fracture  during  transient  cooling,  the  tangential 

stress  Og  is  the  relevant  stress  to  consider  since  the  other  stresses 

are  either  equal  to  or  smaller  than  o^.  Thus,  equating  oQ  to  the  tensile 

strength  S.  of  the  specimen,  the  maximum  temperature  difference  (AT)  , 
t  max 


to  which  the  specimen  can  be  subjected  without  fracture,  can  be  given 
by: 


~ 1Jt('1~")  (1  +  6-82  (7) 

Dl.  Natural  Convection  (Stationary  Specimen). 

Substituting  the  value  of  h  from  equation  (2)  in  equation  (7) 
yields: 


(AT). 


max 


1.451  St(l-v) 
oE 


[1  +  6.82  —■(. - i: 

C  vgYP,C 


f  Pf 


d3nKf1-n(AT)n 


(8) 


It  is  important  to  note  that  AT  in  the  right  hand  side  of  the  equa¬ 
tion  (8)  is  the  temperature  difference  between  the  specimen  wall  and  the 

fluid  at  any  time  whereas  (AT)  is  the  maximum  or  initial  temperature 

max 

difference.  However,  for  small  B,  the  relative  surface  temperature  of 
the  specimen  remains  close1  to  the  initial  temperature  during  the  time 
period  when  the  maximum  stress  is  reached  in  the  specimen.  Thus, 
in  this  case,  AT  can  be  assumed  to  be  approximately  equal  to  (AT)max  for 
thermal  fracture  predictions.  Also,  for  a  small  value  of  B,  the  second 
term  in  the  parenthesis  of  equation  (8)  is  much  greater  than  1  and 
therefore,  (AT)max  can  be  approximated  to  be: 


S  fl-v)  K  v.  .  — 

(AT'  -  c  [—  —  — — ( — - — )n — - — ln+1 

' ° 1  'max  l  aE  ,3n  ey  „  1-n1 

d  f  pf  Kj. 


(9) 


where  is  a  constant  which  depends  upon  the  value  of  (N^  •  Npf) .  For 

a  given  quenching  medium  and  specimen  size,  if  6  is  small, 

1 


S  a-v)  K  n+1 

(AT)  -  - "I 

max  2 


(10) 


6 


D2.  Forced  Convection  (falling  specimen) 


From  equations  (5)  and  (7),  (AT)max  given  by: 


(AT) 


max 


1.451St(l-v) 

aE 


6.82K 


m  m- 


[1+__a{2.lg(^  _1)}  ‘Vf 


Jpf 


1  _  2 
3  3 


1  3 

~  3  ~  2 

p,  d 


(li) 


Similar  to  equation  (10),  for  a  given  quenching  medium  and  specimen 
size,  if  6  is  small,  equation  11  yields: 

P. 


S  U-v) 

(AT)  =  C-.  K  (—  -  1) 

max  3  aE  s 


m 

2 


(12) 


III.  Discussion  and  Conclusions 

The  analyses  presented  above  have  a  number  of  implications  for 
the  quantitative  significance  of  the  experimental  results  of  quench¬ 
ing  experiments  and  for  the  procedures  to  be  followed  in  conducting 
such  experiments.  To  date,  the  interpretation  of  the  quenching  ex¬ 
periments  is  based  on  the  assumption  that  the  heat  transfer  coefficient, 
h,  is  constant  for  a  given  quenching  medium  irrespective  of  the  specimen 
size  and  the  quenching  conditions.  With  the  assumption  of  a  constant 
h  for  a  given  material,  the  magnitude  of  thermal  stress  is  uniquely 
defined  by  the  value  of  the  Biot's  number,  B.  However,  equations  (1), 
(2),  (4)  and  (5)  indicate  that  the  heat  transfer  coefficient,  h,  is 
a  function  of  specimen  size,  specimen  velocity  and  the  instantaneous 
temperature  difference  between  the  specimen  surface  and  the  fluid  medium 
in  addition  to  the  pertinent  fluid  properties.  These  factors  render 
the  interpretation  of  quenching  experiments  considerably  more  complex. 
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especially  when  comparing  the  thermal  quench  data  for  specimens  of  dif¬ 


ferent  materials  having  large  differences  in  (AT)  values  as  well  as 

iricix 

for  specimens  of  a  given  material  which  differ  significantly  in  size. 
Throughout  the  following  discussion,  it  will  be  assumed  that  Biot's 
number  is  sufficiently  small  (6  <  1)  for  the  natural  convection  case 
of  heat  transfer  so  that  over  the  period  of  maximum  stress,  the  change 
in  the  relative  specimen  surface  temperature  is  small.  Therefore,  to 
a  good  approximation  AT  =  (AT)max. 

A.  Effect  of  Specimen  Size  on  (AT) 

r  max 

For  a  heat  transfer  coefficient  independent  of  specimen  size, 
thermal  stress  analysis'  predicts  that  (AT)max  in  a  quenching  ex¬ 
periment  is  inversely  proportional  to  the  speciman  size,  which  can  be 
expressed,  in  general,  as: 


<AT>max  =  A  +  Bd"1  (13) 

where  A  and  B  are  constants. 

The  results  of  the  present  analysis  suggest,  however,  that  this  approach 
requires  modification.  In  the  present  analysis,  the  heat  transfer  co¬ 
efficient,  h,  is  found  to  depend  on  the  specimen  size  as  shown  in  equa¬ 
tion  (5).  Thus,  (AT)max,  based  on  the  size-dependent  heat  transfer  is 
given  by  equation  (11)  and  can  be  written  in  the  general  form: 


(AT) max  =  A  +  Bd 


(14) 


For  forced  convection  heat  transfer  conditions  typically  encountered 
in  laboratory  quenching  experiments,  m  ~  0.466,  which  vields: 

<AT>mnx  "  A  +  Bd'0'6"  (15) 
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Comparison  of  equations  (13)  and  (15)  shows  that  in  forced  convection 

heat  transfer,  where  h  is  given  hv  equation  (5),  (AT)  is  less  dependent 

max 

on  specimen  size  than  for  a  size-independent  heat  transfer  coefficient. 

This  conclusion  is  important  for  the  interpretation  of  experimental 

data  for  (AT)  as  a  function  of  specimen  size.  Equation  (13)  indicates 
ma  x 

that  a  plot  of  (AT)  versus  d  should  yield  an  intercept  equal  to 
max 

-1  * 

A  at  d  =0,  from  which  a  value  of  St(l-v)/aE  can  be  obtained.  How¬ 
ever,  if  h  is  a  function  of  specimen  size,  the  value  of  A  obtained  from 
equation  (13)  is  an  overestimate.  The  present  analysis  indicates  that 

for  a  size-dependent  h,  the  appropriate  value  of  A  is  obtained  by 

,  „  •  .-0.699 

plotting  (AT)  against  d 
max 

It  is  critical  to  note  that  in  forced  convection  heat  transfer 
(falling  specimen),  the  value  of  the  exponent  of  the  diameter  (equa¬ 
tion  5)  depends  on  the  value  of  m.  For  this  reason,  especially  for 
large  specimen  sizes  and/or  low  specimen  velocities  in  highly  viscous 
fluid,  the  nature  of  heat  transfer  must  be  analyzed  quantitatively 
in  order  to  establish  the  proper  size  dependence  for  a  correct  inter¬ 
pretation  of  the  experimental  data. 

For  natural  convection  heat  transfer,  the  size-dependence  of 

(AT)  is  more  complex  for  the  reason  that  the  heat  transfer  coef- 
max 

ficient  h  itself  depends  on  AT.  However,  for  B  <  1,  the  second  term  in 

equation  (8)  is  large  as  compared  to  unity  and  the  change  in  the  relative 

specimen  surface  temperature  is  small  before  the  maximum  stress  is 

reached.  Thus,  AT  -  (AT)  so  that  the  size  dependence  of  (AT) 

’  max  max 

*  The  terms  St(l-v)/aE  and  St(l-v)Ks/aE  will  be  frequently  used  through¬ 
out  discussion  because  these  terms  represent  well  known  thermal  stress 
resistance  parameters  R  and  R' lfi  respectively. 
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can  be  expressed  as: 


(AT)  <*  d"3n/n+1  (16) 

inflx 

For  n  =  •-  (for  laboratory  quenching  conditions),  equation  16  yields: 

(AT)  «  d"0,6  (17) 

max 

which  indicates  that  for  natural  convection  heat  transfer,  a  plot  of 

(AT)  versus  d  should  yield  a  straight  line.  Equation  (17) 
max 

indicates  that  similar  to  the  forced  convection  case,  for  natural  con¬ 
vection  heat  transfer  with  a  size-dependent  heat  transfer  coefficient, 

(AT)  is  less  dependent  on  specimen  size  than  for  a  size-independent 
max 

heat  transfer  coefficient  expressed  by  equation  (13). 

For  natural  convection  heat  transfer,  caution  should  be  exercised 

in  extrapolation  of  a  plot  of  (AT)  versus  d  to  high  values  of  d. 

max 

Such  an  extrapolation  includes  values  of  the  Biot's  number  3  >>  1  which 
violates  the  assumption  AT  -  (AT)  Unless  accompanied  by  extensive 

numerical  analysis  for  the  variation  of  h  with  AT  during  the  quench 
test,  such  an  extrapolation  is  not  expected  to  result  in  a  meaningful 
value  of  the  intercept.  For  this  reason,  a  reliable  determination  of 
the  quantity  St(l-v)/aE  should  be  based  on  quenching  tests  using  fal¬ 
ling  specimens  rather  than  stationary  specimens. 

B.  Dependence  of  h  on  AT 

For  heat  transfer  by  natural  convection  (stationary  specimen), 
the  heat  transfer  coefficient  h  is  a  function  of  the  instantaneous 
temperature  difference  between  the  specimen  surface  and  the  quenching 
fluid  (see  equation  2).  This  implies  that  different  materials  with 
varying  values  of  (AT)max  ate  subjected  to  different  heat  transfer 


conditions.  Due  to  this  reason,  for  dissimilar  materials,  no  direct 
proportionality  is  expected  between  (AT)max  and  the  quantities 
S.(l-v)/aE  or  S. (l-v)K  /aE.  To  examine  the  role  of  these  quantities 

t  L  o 

in  establishing  (AT)  ,  equation  (9)  (for  6<1  and  n=%)  can  be  written 

max 

as: 

4/5 


(AT)  “  {S  (l-v)K  /aE)  (18) 

max  t  s 

Equation  18  indicates  that  the  relative  difference  in  thermal  stress 
resistance  of  different  materials  as  measured  bv  (AT)  in  a  quenching 
experiment  is  less  than  the  relative  difference  in  the  respective 
quantities  of  St (l-v)Ks/aE.  In  order  to  obtain  linearity  for  the 
purpose  of  inter-  or  extrapolation,  it  is  recommended  that  a  plot  be  ob- 

5/4 

tained  of  (AT)  versus  S  (l-v)K  /aE.  Conversely,  linearity  can  be 

max  l  s 

also  achieved  by  plotting  (AT)max  against  the  thermal  stress  resistance 


parameter,  R^  given  by: 


Rn  =  {St(l-v)Ks/aE} 


4/5 


(19) 


appropriate  for  the  natural  convection  heat  transfer  with  a  AT~dependent 
heat  transfer  coefficient  and  low  8. 


C.  Effect  of  Specimen  Density  on  (AT) 

max 

As  indicated  by  equations  (3e),  (5)  and  (11),  for  a  quenching  experiment 

involving  forced  convection  (falling  specimen)  the  specimen  velocity  and 

therefore  the  heat  transfer  coefficient  and  (AT)  are  function  of 

max 

the  ratio  of  specimen  to  fluid  density.  For  this  reason,  different 

materials  with  identical  values  for  the  quantity  St(l-v)Kg/aE  can  still 

exhibit  different  values  of  (AT)  due  to  possible  difference  in 

max 

densities.  This  effect  should  be  taken  Into  account  in  the  interpre¬ 
tation  of  experimental  data  for  (AT)  .  The  direct  transfer  of  the 
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relative  thermal  stress  resistance  for  different  materials  as  measured 

by  (AT)  in  a  falling  specimen  method  to  the  experimental  conditions 
max 

where  the  specimen  densities  do  not  affect  the  magnitude  of  thermal 

stress  resistance  could  lead  to  considerable  error.  To  avoid  such 

error,  it  is  proposed  that  the  values  of  (AT)  be  correlated  with 
r  r  max 

the  thermal  stress  resistance  parameter  given  by: 


S 

The  value  of  t 
known  values  of 


S  (l-v)K  p 

R  =  {— - -}{  —  -  1] 

f  aE  pf 


m 

2 


(20) 


(l-v)K 


s  can  be  easily  extracted  from  equation  (20)  for 


aE 

and  m.  It  should  be  noted  that  the  quantity 
Pf 


depends  on  the  value  of  the  exponent  m.  Thus,  the  nature  of  the  heat 
transfer  in  a  thermal  quench  experiment  must  be  assessed  quantitatively 
in  terms  of  the  pertinent  fluid  properties,  specimen  size  and  terminal 
velocity  before  a  valid  correlation  between  the  theory  and  experiments 
can  be  obtained. 

The  effect  of  specimen  density  on  (AT)  also  can  play  a  significant 

role  in  the  effect  of  porosity  on  thermal  stress  resistance  for  a  given 

material.  This  effect  can  be  measured  bv  changes  in  (AT)  in  the 

max 

falling  specimen  method.  An  increase  in  porosity,  (i.e.,  a  decrease  in 
density)  for  a  porous  specimen  whose  density  approaches  the  fluid  den¬ 
sity  will  result  in  an  increase  in  (AT)  due  to  large  increase  in 
„  max 

_  in 

P  ir  ~ 

{— 8  -  ll  L .  In  contrast,  for  specimens  with  densities  much  higher 
Pf 

than  the  fluid  density,  the  increase  in  porosity  is  expected  to  de¬ 
crease  (AT)  due  to  greater  decrease  in  the  value  of  S,(l-v)K  / aE. 

max  L  s 

In  general,  the  effect  of  porisitv  on  thermal  stress  resistance  will 
depend  on  the  method  of  measurement.  For  valid  comparison  of  the 
results  obtained  by  different  testing  methods  including  the  falling 
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specimen  method,  the  effect  of  specimen  density  must  be  taken  into 


account • 

D.  Recommendations  for  Quantitative  Quenching  Experiments 

The  general  analytical  results  of  the  present  study  also  lead  to 
a  number  of  recommendations  for  procedures  to  be  followed  in  quenching 
experiments. 

As  indicated  by  equations  (2)  and  (5)  for  both  natural  convection 

(stationary  specimen)  and  forced  convection  (falling  specimen),  the 

heat  transfer  coefficient  is  a  function  of  properties  of  the  fluid 

medium.  Many  of  these  fluid  properties,  especially  the  viscosity,  are 

strongly  temperature-dependent.  For  this  reason  it  is  critical  to 

define  the  temperature  of  the  quenching  bath  in  reporting  values  for 

(AT)  .  Because  of  the  relatively  rapid  variation  of  fluid  properties 
max 

with  temperature,  reporting  bath  temperature  as  "room  temperature"  which 
can  vary  by  as  much  as  five  to  ten  degrees,  must  not  be  considered  adequate. 
In  order  to  facilitate  comparisons  of  data  obtained  by  different  investi¬ 
gators  or  laboratories,  it  is  recommended  that  bath  temperature  be  con¬ 
trolled  and  maintained  at  an  internationally  agreed-upon  value.  A 
temperature  recommended  here,  which  should  be  easily  accomodated  in  most 
studies,  is  35°C  +  0.5°C.  This  temperature  is  slightly  above  the 
room  temperature  in  the  vast  majority  of  cases  and  can  be  maintained 
with  minimum  requirement  for  heating  and  control  devices. 

Another  recommendation  for  the  quantitative  evaluation  of  (AT) 

'  max 

data  obtained  by  the  falling  specimen  method  concerns  the  velocity  with 
which  the  specimen  falls  through  the  fluid.  The  present  analysis  is 
based  on  the  assumption  that  the  specimen  falls  through  the  fluid  with 
its  terminal  >/elocity.  Therefore,  to  avoid  anv  undue  complexity  in 


n 


the  data  analysis  it  is  desirable  that  the  specimen  should  move  with 
its  terminal  velocity  over  the  time-period  between  entering  the  fluid 
till  the  time  of  failure.  Any  deviation  from  the  terminal  velocity, 
especially  at  the  time  of  specimen  entry  into  the  fluid,  will  introduce 
transient  effects  in  heat  transfer.  Such  transient  effects  can  be 
eliminated  if  the  specimens  were  to  hit  the  fluid  surface  with  its 
terminal  velocity.  This  can  be  accomplished  by  dropping  the  specimen 
from  a  pre-determined  optimum  height  such  that  the  velocity  of  the  free 
fall  (in  air)  corresponds  to  the  terminal  velocity  when  it  hits  the 
fluid  surface.  This  height  can  be  determined  as  follows: 

For  a  free  fall  through  height  H  in  air,  the  specimen  velocity 
Va  (neglecting  air  friction)  just  before  it  touches  the  fluid  surface 
is  given  by, 

Va  =  /2iH  (21) 

The  terminal  velocity  Vt  of  a  circular  rod  specimen  falling  horizontal¬ 
ly  through  t'-e  fluid  can  be  expressed  as, 


V 

t 


2f, 


g  (—  -  l)  d] 


(22) 


where  fp  is  a  function  of  the  Reynold's  number  as  shown  in  Table  I. 

For  the  optimum  height,  V  =  V  ,  and  therefore,  from  equations  (21)  and 

3  t 

(22),  the  optimum  height  H  can  be  given  by: 


H  =  (--a- 

m  4fD  pf 


1  )d 


(23) 


For  fluid  properties  and  specimen  velocities  commonly  encountered  in 

2  5 

laboratory  experiments,  4x10  <  <  4x10  which  corresponds  to  a 


value  of  f  =  0.75  such  that 
D 


H 

m 


1.05  (J?.  _  i)d 


(24) 


t 
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The  value?  of  for  typical  ceramics  and  specimen  sizes  encountered 
in  practice  correspond  to  heights  of  a  few  centimeters.  In  general 
for  a  given  fluid  and  specimen  of  arbitrary  size  and  density,  the 
Reynold's  number  needs  to  be  established  before  the  optimim  height 
can  be  calculated. 

Finally,  it  is  recommended  that  a  careful  consideration  should 
be  given  to  the  depth  of  the  quenching  bath  used  for  the  falling 
specimen.  In  general,  the  thermal  stresses  do  not  reach  their  maximum 
value  at  the  instant  the  specimen  hits  the  fluid  surface.  For  fluids 
with  relatively  low  heat  transfer  coefficient  values  in  combination 
with  high  thermal  conductivity  and  small  size  of  the  specimen  (i.e.. 
low  Biot's  number),  the  time  to  maximum  stress  could  be  on  the  order 
of  a  few  seconds.  Under  these  conditions,  the  specimen  will  take  less  time 
to  hit  the  bottom  of  the  bath  than  the  time  required  to  reach  the 
maximum  stress.  Consequently,  the  heat  transfer  coefficient  would  not 
be  constant  over  the  time  period  required  to  reach  maximum  stress  which 
could  introduce  considerable  error  in  quantitative  analvsis  of  (AT) 

max 

In  addition,  a  meaningful  comparison  cannot  be  made  of  (AT)  data 
obtained  from  different  laboratories  where  the  bath  depts  could  be 
quite  different.  Thus,  the  heat  transfer  properties  of  the  fluid 
medium  and  the  thermo-mechanical  response  of  the  specimen  must  be 
well  understood  for  any  quantitative  interpretation  of  thermal  quench 
data . 
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TABLE  I 


Drag  coefficient  for  laminar  flow  past  immersed  body^2 


Surface 

Reynold's  number 

Drag 

Coefficient 

Nd  <0.5 

Re 

fD  = 

8it 

Circular  cylinder 
with  axis  per¬ 
pendicular  to  flow. 

Re (2 . 002-lnRe) 

10  * 
p  0.6 

Re 

0. 5<N„  <10 

Re 

fD  = 

* 

10<N„  <40 

Re 

Ml 

o 

II 

5 

V  °-3 

Re 

40<N„  <400 

Re 

fD  = 

5.87 

„  0.343 

Re 

400<h'  <4x10^ 

Re 

fD 

J. 

0.75" 

*  These  equations  are  obtained  by  curve  fitting  the  ex- 
permental  data. 
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ABSTRACT 
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The  effect  of  a  spatial  variation  of  the  thermal  conductivity  on  the 
magnitude  of  the  maximum  tensile  thermal  stress  in  a  solid  circular  cylinder 
subjected  to  sudden  convective  heating  was  calculated  by  finite  element 
methods.  The  general  results  show  that  by  lowering  the  thermal  conductivity 
in  the  surface  region  of  the  cylinder,  the  magnitude  of  the  maximum  tensile 
thermal  stress  at  the  center  of  the  cylinder  is  reduced  significantly.  Also, 
the  time  to  maximum  stress  is  increased  substantially  as  well.  It  is  also 
shown,  that  a  change  in  the  thermal  diffusivity  affects  the  time  of  maximum 
tensile  thermal  stress  but  not  its  magnitude. 
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I.  Introduction 

Brittle  ceramic  materials  for  high-temperature  applications  are 
highly  susceptible  to  catastrophic  failure  due  to  thermal  stresses  which 
result  from  constrained  thermal  expansions^.  A  detailed  analysis  of  the 
magnitude  and  distribution  of  the  thermal  stresses  and  selection  of 
materials  with  good  thermal  stress  resistance  is  an  essential  feature 
of  the  design  and  construction  of  structures  or  components  operating 
at  non-ambient  environments  The  selection  of  material  with  optimum 
thermal  stress  resistance  can  be  based  on  "thermal  stress  resistance" 
parameters  or  "f igures-of-merit" ,  derived  for  a  number  of  heat  trans¬ 
fer  conditions,  nerforraance  criteria  and  failure  modes  as  reviewed  re¬ 
cently^. 

Unfortunately,  even  the  best  material,  selected  on  this  basis,  may  be 
unable  to  withstand  thermal  stresses  encountered  in  severe  thermal  environ¬ 
ment.  A  partial  solution  of  this  problem  is  to  develop  new  or  improved 
materials.  Another  approach  is  to  reduce  the  magnitude  of  thermal  stress 
to  which  any  given  material  may  be  subjected.  In  respect  to  this  latter 
approach,  it  should  be  noted  that  the  magnitude  of  maximum  thermal  stress 
in  a  structure  or  component  is  a  function  of  the  maximum  difference  of 
temperature  in  the  structure  and  the  manner  in  which  the  temperature  is 
distributed.  In  principle  then,  the  magnitude  of  thermal  stress  can  be 
reduced  by  modifications  of  the  temperature  distribution.  One  technique 
by  which  this  can  be  accomplished  is  to  introduce  a  spatially  varying  thermal 

Q 

conductivity  in  the  structure.  In  a  recent  study,  Hasselman  and  Youngblood 
demonstrated  that  a  spatially  varying  thermal  conductivity  in  a  hollow  cy¬ 
linder  can  cause  a  significant  reduction  in  the  magnitude  of  the  tensile 
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thermal  stress  under  conditions  of  steady-state  heat  flow.  It  is  the 
purpose  of  the  present  study  to  investigate  whether  a  spatially  varying 
thermal  conductivity  will  also  reduce  the  magnitude  of  the  maximum  tensile 
thermal  stress  under  transient  conditions. 

II.  Analysis 

An  infinitely  long  solid  circular  cylinder  is  selected  as  a  geometry 
appropriate  for  the  present  study.  The  cylinder  is  subjected  to  sudden 
heating  by  convective  Newtonian  heat  transfer,  from  an  initial  temperature, 
Tq  to  the  final  temperature  of  the  surrounding,  T^. 

The  spatial  variation  of  the  thermal  conductivity  is  arbitrarily 
chosen  to  be  of  the  form: 

K<r)  *  KqU  -  C(r/a)n)  (1) 

where  n  and  C  (0  <  C  <  1)  are  constants,  r  is  the  radial  coordinate,  a 
is  the  radius  of  the  cylinder  and  is  the  thermal  conductivity  cor¬ 
responding  to  the  value  for  the  uniformly  distributed  thermal  conducti¬ 
vity.  Eq.  1  indicates  that  at  r  »  0,  K  =  Kq  for  any  value  of  C  and  n. 

For  0  <  r  <  a,  K  <  Kq,  which  implies  that  the  spatial  variation  of  the 
thermal  conductivity  is  such  that  the  thermal  conductivity  is  decreased 
with  respect  to  its  value  at  the  center.  The  lowest  value  of  thermal 
conductivity  K  -  Kq(1  ~  C)  occurs  at  r  =  a.  In  order  to  establish  the 
effect  of  the  spatial  variation  of  the  thermal  conductivity  on  the 
magnitude  of  thermal  stress  alone,  all  other  pertinent  material  properties 
are  assumed  to  be  spatially  and  temperature  invariant. 

It  should  be  noted  that  for  a  spatially  uniform  thermal  conductivity 
closed  form  solutions  for  the  transient  temperatures  and  thermal  stresses 
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in  solid  circular  cylinder  are  available  in  the  literature  ’  .  However, 

calculations  for  the  transient  temperature,  and  thermal  stresses  in  a 
solid  circular  cylinder  with  a  spatially  varying  thermal  conductivity 
distribution  requires  numerical  methods  such  as  finite  element  analysis 
selected  for  the  present  study. 

The  differential  equation  governing  the  present  heat  transfer  problem 
is: 


9  r  t,  /  \  3T,  3T 

-  (rK(r)  =  Pc  ^ 


(2) 


where  p  and  c  are  the  density  and  specific  heat  of  the  cylinder  and  T 
is  the  instantaneous  temperature.  For  the  present  problem  the  boundary 
conditions  are: 


II 

3r 


0  at  r 


0 


(3a) 


{ rK(r)  ^}+  rh(T  -  TJ  =  0  at  r  =  a  (3b) 

Since  the  heat  flow  in  the  cylinder  occurs  in  the  radial  direction 
only,  the  transient  temperatures  can  be  obtained  with  a  consistent  one¬ 
dimensional  (radial)  element  formulation  with  two  nodes  per  element  and 
continuity  of  only  the  nodal  temperature  across  the  inter-element  boundaries. 
Using  a  finite  element  code  compiled  for  this  program,  the  temperatures 
are  obtained  with  a  radial  element  configuration  composed  of  33  elements. 

Because  of  the  radial  and  tangential  symmetry  it  is  sufficient  to  analyze 
a  sector  of  the  circular  cross-section  for  the  calculation  of  the  thermal 
stresses  from  the  temperature  distribution.  This  sector  is  being  modelled 
as  an  assemblage  of  23  quadr ilateral  elements  and  imposing  the  condition 
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that  Che  displacement  field  is  radial  only.  Conditions  of  plane  strain 
and  linear  elastic  behavior  are  assumed  throughout.  For  purposes  of  a 
slight  computational  advantage,  the  quadrilateral  elements  can  be  generated 
from  an  assemblage  of  four  triangular  elements  followed  by  a  static  condensa¬ 
tion  to  eliminate  the  common  internal  node.  The  reader  is  referred  to  the 
text  of  Zienkiewicz^  for  further  details  on  the  finite  element  method. 

The  validity  of  this  finite  element  program  was  established  by  verifying  the 
analytical  results  of  Jaeger^,  which  resulted  in  excellent  agreement. 

For  convenience  and  for  purposes  of  comparison,  the  numerical  results 

•k 

are  reported  in  terms  of  the  non-dimensional  thermal  stress,  o  *  o(l-v)/aEAT, 

*  2 

the  non-dimensional  time,  x  «  cct/a  where  <  is  the  thermal  diffusivity,  iWoc 

and  t  is  the  time,  the  non-dimensional  radial  coordinate,  r/a  and  the  Biot 

number,  8  ■  ah/K  . 

o  o 

As  a  final  remark,  it  should  be  noted  that  only  sufficient  numerical 
results  are  obtained  to  indicate  the  general  trends,  rather  than  making 
calculations  of  the  thermal  stresses  for  specific  thermal  shock  experiments. 

III.  Numerical  Results  and  Discussion 

Since  brittle  ceramic  materials  are  far  more  susceptible  to  failure 
in  tension  than  in  compression,  in  this  study  the  tensile  thermal  stresses 
are  of  primary  interest.  For  a  solid  circular  cylinder  subjected  to 
convective  heating,  the  tensile  thermal  stresses  have  their  maximum  at 
the  center  of  the  cylinder  (r  ■  0).  Figure  1  shows  the  time  dependence 
of  the  maximum  value  of  the  tensile  thermal  stress  for  values  of  C  *  0.9 
and  0.99  and  n  -  5  and  10,  for  the  Biot  number  &Q  -  5.  For  the  purpose 
of  comparison  the  maximum  tensile  stress  for  a  spatially  invariant  thermal 
conductivity  is  included  in  Figure  1  as  well.  Comparison  of  the  curves 
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shows  that  the  effect  of  the  spatially  varying  thermal  conductivity  is 

to  decrease  the  magnitude  of  the  tensile  thermal  stress  at  any  value  of 

* 

t  by  a  significant  fraction.  Noteworthy  also  is  that  the  spatial  varia¬ 
tion  in  thermal  conductivity  leads  to  a  significant  increase  in  the 

* 

values  of  x  at  which  the  tensile  thermal  stresses  reach  their  maximum 
values . 

Figures  2a  and  2b  show  the  temperatures  and  the  thermal  stresses 
as  a  function  of  the  radial  distance  at  the  instant  at  which  the  tensile 
thermal  stresses  reach  their  maximum  values  for  the  values  of  C,  n  and  £ 

o 

shown  in  Fig.  1.  As  indicated  by  Fig.  2a,  the  effect  of  the  spatial 
variation  of  thermal  conductivity  is  to  cause  a  significant  change  in 
the  temperature  distribution  in  support  of  the  hypothesis  on  which  this 
study  is  based.  Figure  2b  indicates  that  the  decrease  in  the  value  of 
maximum  tensile  stress  at  r/a  =  0  is  accompanied  by  an  increase  in  the 
value  of  compressive  stress  in  the  surface  at  r/a  =  1. 

For  8Q  *  5,  Fig.  3  gives  the  values  of  the  maximum  tensile  thermal 
stress  for  a  range  of  values  of  C  and  n.  Generally,  the  higher  the  value 
of  C  (i.e.,  the  lower  the  relative  thermal  conductivity  at  the  surface), 
the  lower  is  the  value  of  tensile  thermal  stress.  However,  for  a  given 
value  of  C,  the  maximum  value  of  thermal  stress  is  a  function  of  n,  with 
a  general  inverse  dependence  between  C  and  n. 

★ 

For  6q  =  5,  Fig.  4  shows  the  values  of  time,  i  as  a  function  of 

C  and  n,  at  which  the  tensile  thermal  stresses  reach  their  maximum  values. 

★ 

For  high  values  of  C  and  lower  values  of  n,  the  values  of  t  are  signifi- 

k 

cantly  higher  than  the  value  of  t  -  0.11  for  constant  thermal  conductivit 

k 

Figure  5  shows  the  maximum  tensile  thermal  stress,  a  for  C  *  0.9 
and  n  ■  3  as  a  function  of  the  Biot  number,  p  .  For  comparison,  the 
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value  of  o  for  constant  thermal  conductivity  is  included  in  Fig.  5  as 

well.  These  results  indicate  that  the  relative  reduction  in  the  maximum 

value  of  the  tensile  thermal  stress  is  greater  for  the  higher  values  of 

the  Biot  number.  The  excessive  computer  time  required  precluded  a  com- 

* 

plete  optimization  of  the  values  of  a  as  a  function  of  SQ  for  all 

* 

possible  variations  of  C  and  n.  As  a  result,  the  values  of  a  as  a 
function  of  Bq  shown  in  Fig.  5  could  even  be  lower  than  indicated. 

Additional  calculations  were  made  for  the  magnitude  of  thermal  stresses 
for  the  same  values  of  C,  n  and  Bq  as  indicated  in  Fig.  1,  but  for  convec¬ 
tive  cooling  rather  than  heating.  The  results  showed  that  the  magnitude 
of  the  stresses  are  identical  to  the  heating  case  but  of  opposite  sign. 

This  result  implies  that  in  case  of  cooling,  for  the  spatial  variation  of 
thermal  stresses  in  fact,  will  be  increased.  It  appears  then  that  for  a 
decrease  in  the  tensile  thermal  stresses  on  cooling,  a  thermal  conductivity 
distribution  must  be  chosen  which  differs  from  the  one  given  by  eq.  1. 

A  further  calculation  concerned  the  time  at  which  the  tensile  stress 
reaches  its  maximum  value.  For  a  constant  thermal  conductivity  distribution, 
the  magnitude  of  the  thermal  diffusivity  controls  the  time  at  which  the 
stresses  reach  their  maximum  value,  but  not  their  magnitude.  To  cneck  whether 
this  still  holds  for  a  spatially  varying  thermal  conductivity,  the  maximum 
tensile  stresses  were  calculated  for  C  =  0.999,  n  »  20  and  two  values  of 
pc  *  0.997  and  1.994.  Here  the  value  of  thermal  diffusivity  was  varied  by 
a  factor  of  2  keeping  K  constant.  The  above  values  of  C  and  n  were  chosen 
deliberately  to  obtain  an  extreme  variation  in  thermal  conductivity.  The 
numerical  results  show  in  Fig.  6  indeed  indicate  that  the  magnitude  of 
thermal  stress  is  not  affected  by  changes  in  the  thermal  diffusivity.  As 
for  the  constant  thermal  conductivity  case,  the  inverse  dependence  of  time 


of  maximum  stress  and  thermal  diffusivity  is  maintained. 


The  general  results  of  this  study  may  have  practical  significance.  For 
constant  thermal  conductivity,  increase  in  thermal  stress  resistance  (for 
all  other  properties  being  constant)  requires  an  increase  in  thermal  con¬ 
ductivity.  For  nominally  pure,  fully-dense  materials,  such  increase  in 
thermal  conductivity  may  be  hard  to  obtain.  Increasing  thermal  stress  re¬ 
sistance  by  lowering  the  thermal  conductivity  in  the  outer  surface  as  sug¬ 
gested  by  the  present  result  should  be  much  easier  to  obtain.  Such  spatial¬ 
ly  varying  thermal  conductivity  may  be  accomplished  by  spatially  varying  pore 
content,  non-uniformly  distributed  impurities  or  solid-solution  content, 
composites  with  graded  compositions,  a  spatially  non-uniform  degree  of  cry¬ 
stallization  in  glass  ceramics  and  perhaps  by  a  host  of  other  special  techni¬ 
ques  which  will  depend  on  the  specific  material  under  consideration. 

The  data  shown  in  Fig.  5  should  be  of  particular  siginficance  to  high- 
temperature  technology.  The  large  decreases  in  tensile  thermal  stresses  are 
obtained  at  the  high  values  of  the  Biot  number,  Bq.  High  values  of  Bq 
generally  imply  high  heat  flux,  large  structure  size,  or  low  thermal  con¬ 
ductivity.  Under  these  conditions,  brittle  ceramic  materials  are  particularly 
susceptible  to  thermal  fracture,  because  of  the  generally  high  magnitude  of 
thermal  stress.  For  this  reason,  the  results  of  this  study  are  most  appropriate 
for  those  environments  for  which  the  thermal  fracture  problem  is  more  severe. 

In  particular,  thermal  fatigue  life,  because  of  its  extreme  stress  dependence, 
should  exhibit  orders  of  magnitude  improvement. 

It  will  be  recalled  that  for  the  purpose  of  reporting  convenience,  the 
values  of  maximum  tensile  thermal  stresses  were  compared  at  a  given  value 
of  the  Biot  number  (BQ)  corresponding  to  the  value  of  thermal  conductivity, 

Kq.  However,  for  a  spatial  variation  of  the  thermal  conductivity  described 
by  eq.  1,  the  relative  thermal  conductivity  at  the  surface  is  decreased  by  a 
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factor  equal  to  (1  -  C) ,  with  a  corresponding  increase  in  the  Biot  number 
by  the  factor  1/1-C.  If  the  comparison  of  the  stresses  would  take  this 
change  in  Biot  number  into  account  as  well,  the  reduction  in  thermal  stress 
level  due  to  the  spatially  varying  thermal  conductivity  is  even  greater. 

For  this  reason,  comparison  of  the  stresses  on  the  basis  of  Kq,  in  fact, 
yields  a  conservative  estimate  of  the  stress  reductions  which  can  be  ob¬ 
tained. 

As  noted  earlier,  the  spatial  variation  in  the  thermal  conductivity 
described  in  eq.  1  in  convective  cooling  leads  to  an  increase  in  the 
magnitude  cf  tensile  thermal  stresses  in  the  surface  as  compared  to  the 
stress  level  for  constant  thermal  conductivity.  Under  these  conditions, 
then,  the  spatial  variation  in  thermal  conductivity  can  lead  to  a  decrease 
in  thermal  stress  resistance.  This  problem,  however,  can  be  overcome  by 
superposing  on  the  thermal  stresses  an  internal  stress  distribution  which 
is  highly  compressive  in  the  surface  and  has  a  relatively  low  tensile 
stress  in  the  interior  of  the  cylinder.  This  type  of  internal  stress  dis¬ 
tribution,  frequently  encountered  in  glassy  materials  strengthened  by  ion- 
exchange,  should  be  particularly  effective  for  materials  which  fail  in 
tension  from  surface  flaws.  In  this  manner,  then,  by  combining  a  spatially 
varying  thermal  conductivity  (with  low  surface  values)  with  surface-com¬ 
pression  strengthening,  improvements  in  thermal  stress  resistance  for  both 
convective  heating  and  cooling  can  be  obtained. 

In  summary,  it  can  be  concluded  that  spatial  variations  in  thermal 
conductivity  can  lead  to  significant  decreases  in  the  magnitude  of  the 


tensile  thermal  stresses  generated  during  transient  convective  heating. 
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Figure  6.  Transient  maximum  tensile  thermal  stress  in  solid  circular  cylinder 
with  spatially  varying  thermal  conductivity  (('.  =  0.099,  n  -  -’O)  [or 
two  values  of  tile  specific  heat  per  unil  volume. 
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Fracture  of  brittle  materials  subjected  to  thermal  stress  can 
occur  by  unstable  or  stable  crack  propagation^  jn  terms  of  the 

strength  retained  as  a  function  of  temperature  difference  of  thermal 
shock  such  as  in  quenching  tests,  unstable  crack  propagation  is  mani¬ 
fested  by  a  precipitous  decline  in  strength  at  a  critical  temperature 
difference,  AT^.  Stable  crack  propagation,  in  contrasty  is  characterized 
by  a  monotonic  decrease  in  strength  wi*-h  increasing  temperature  dif¬ 
ference.  Generally,  if  sufficient  data  for  the  retained  strength  are 
obtained  over  the  total  range  of  temperature  difference,  for  most 
literature  data  the  stable  or  unstable  nature  of  crack  propagation 
can  be  established  unequivocally. 

Nevertheless,  the  usual  scatter  in  material  properties  relevant 
to  thermal  stress  fracture  can  introduce  a  complexity  in  the  characteri¬ 
zation  of  the  failure  mode.  In  this  respect,  due  to  variation  in  crack 
depth  and  geometry  and/or  local  variations  in  fracture  toughness, 
strength  represents  the  most  significant  variable.  For  a  given  set 
of  specimens  subjected  to  a  quenching  test,  those  with  a  higher  value 
of  tensile  strength  will  require  a  higher  temperature  difference  to 
induce  failure  than  those  with  the  lower  values  of  tensile  strength. 

As  a  direct  result,  at  values  of  AT  near  ATc  for  a  set  of  specimens 
only  a  fraction  within  the  set  may  undergo  fracture,  with  the  remainder 
surviving  intact.  For  specimens  which  undergo  relatively  large  strength 
losses  on  fracture  and  for  which  the  coefficients  of  variations  of 
the  surviving  and  failed  specimens  are  relatively  small,  the  distri¬ 
bution  of  the  strength  values  for  all  specimens  shows  a  pronounced 
bi-modal  distribution/  Clearly,  in  this  case  it  is  inappropriate  to 
judge  the  strength  loss  behavior  on  the  basis  of  the  averages  of  all 
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specimens  tested  at  a  given  value  of  AT.  Instead,  the  strength  loss 
behavior  properly  should  be  judged  on  the  averages  of  each  mode  of  the 
bi-modal  distribution. 

However,  the  situation  can  arise  that  in  unstable  fracture,  the 
distributions  of  the  two  modes  of  the  bi-modal  distribution  show  signi¬ 
ficant  overlap,  so  that  the  clear-cut  bi-modal  nature  becomes  less 
apparent.  This  is  expected  to  be  the  case  for  materials  with  moderate 
strength  which  undergo  only  a  small  relative  strength  loss.  In  this 
case,  the  determination  of  the  thermal  stress  failure  mode  becomes 
ambiguous,  and  so  could  lead  to  an  erroneous  interpretation.  The 
purpose  of  the  present  note  is  to  illustrate  these  observations  by 
data  for  the  thermal  shock  behavior  of  a  reaction  sintered  silicon 
nitride. 

Quenching  experiments  were  carried  out  with  commercial  reaction- 
sintered  silicon-nitride  with  a  density  of  2.31  g/cm^  and  an  average  size 
of  a  small  pores  of  0.05  pm,  but  with  single  pores  up  to  30  pm.  Rec¬ 
tangular  rods  6  cm  long  and  with  cross-section  of  6  mm  by  6  mm  were 

used.  The  rods  in  the  as-nitrided  form  were  quenched  from  a  range  of 

§ 

higher  temperatures  into  a  silicone  oil  bath  at  R.T. ,  followed  by  a 
strength  test  in  4-pt  bending  with  a  total  span  of  40  mm  and  20  mm  bet¬ 
ween  the  central  loading  points  at  a  cross-head  speed  of  0.05  mm/min  in 
laboratory  air  at  R.T. 

The  individual  data  points  for  all  specimens  are  shown  in  Fig.  1. 
Table  I  lists  the  average  and  standard  deviations  for  each  set  of  speci¬ 
mens  tested  at  a  given  value  of  AT.  Note  the  relatively  large  scatter 

§  Dow  Corning,  Type  705,  125  Ct. 
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in  strength  data  in  the  as-nitrided  (AT  =  0)  specimens.  Scanning  electron 
microscopy  revealed  that  the  relatively  large  scatter  most  likely  is  due 
to  large  structural  inhomogenities  in  the  near-surface  layer  of  about 
200  um  in  thickness,  shown  in  Fig.  2.  Major  decreases  in  strength  occur 
over  the  temperature  range  725  to  975°C.  In  this  range,  the  strength 
data  scatter  exceeds  that  at  R.T.  For  AT  975°C,  the  mean  values  show 
no  statistical  significant  difference  with  R.T.  data.  For  ATc  >  975°C, 
there  is  a  significant  difference  between  the  mean  values  of  these  tem¬ 
peratures  and  R.T.  at  the  99.5%  significance  level.  Note  that  for 
725°C  <  AT  <  975°C  there  is  no  visually  obvious  evidence  for  the  existence 
of  a  bi-modal  strength  distribution.  As  a  result,  the  relative  strength 
behavior  must  be  judged  by  the  averages  of  all  specimens  tested  at  given 
AT  values.  It  appears  then,  in  the  absence  of  additional  information, 
that  failure  in  these  specimens  must  have  occurred  by  crack  propagation 
in  the  stable  mode. 

The  validity  of  the  above  conclusion  can  be  examined  on  the  basis  of 

an  approximate  analysis.  As  derived  for  a  two-dimensional  model  with  N 

cracks  of  half-length  %  per  unit  area,  uniaxially  constrained  in  the 

direction  perpendicular  to  the  cracks  and  constrained  from  shrinkage 

on  cooling,  the  critical  temperature  difference  AT^  required  for  crack 

2 

instability  and  propagation  is  : 

ATc  =  (2G/a2El)5i  (1+2ttNJ2)  0) 

where  G  is  the  fracture  surface  energy  required  to  create  unit  area  of 

fracture  surface,  a  is  the  coefficient  of  thermal  expansion  and  E  is 

2 

Young's  modulus  of  the  crack-free  plate.  In  stable  fracture  the  terra 
2 

2irN£.  dominates  such  that  eq.  1  can  be  written: 

AT  -  2nN(2G/u2E)i5f3/2  (2) 

c 
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By  noting  that  strength  is  proportional  to  the  reciprocal  square 
root  of  crack  length,  from  eq.  2  the  ratio  of  strength  values  o ^ ,  and  Oj 
of  specimens  subjected  to  temperature  differences  AT^  and  AT2  can  be 
derived  to  be: 

(o2/o1)  =  (AiyAT,)173  (3) 

Note  that  eq.  3  predicts  the  relative  strength  behavior  in  stable 
crack  propagation  and  therefore  does  not  depend  on  the  specific  geometry 
of  the  test  specimen,  with  the  exception  that  the  derivation  of  eq .  1 
assumes  the  absence  of  the  effects  of  crack  interaction  and  a  specific 
dependence  of  the  effective  Young's  modulus  of  the  plate  and  crack 
length.  For  this  reason,  eq.  3  at  least  to  a  first  approximation  can 
be  applied  to  the  present  specimens. 

For  the  present  data  with  crack  propagation  occurring  over  the 

1/1 

range  of  725°C  <  AT  <  975°C,  (AT-j/AT2)  =  0.91.  The  corresponding 
observed  strength  ratio  02/0^  -  0.16,  which  is  much  smaller  than  the 
predicted  value.  This  large  discrepancy  between  predicted  and  observed 
strength  behavior  indicates  that  the  strength  loss  is  much  too  great 
over  the  temperature  difference  between  725  to  975°C  for  crack  propaga¬ 
tion  to  have  occurred  in  a  stable  manner.  For  this  reason,  in  spite  of 
the  absence  of  direct  experimental  evidence,  it  must  be  concluded  that 
at  least  the  majority  of  specimens  must  have  failed  by  unstable  crack 
propagation. 

The  correctness  of  the  latter  conclusion  was  verified  experimental¬ 
ly  by  carrying  out  a  similar  series  of  quenching  experiments  as  shown 
in  Fig.  1,  using  an  additional  set  of  specimens  of  the  same  original 
batch  which  were  sur face-ground  to  remove  the  structural  inhnmogeni t ies 
in  the  surface.  The  grinding  consisted  of  removing  a  surface  layer 


of  approximately  500  pm  with  a  diamond  wheel  with  particle  size  of  64  urn 
at  a  surface  speed  of  18  m/min  in  a  direction  parallel  to  the  length  of 
the  specimen.  The  mean  surface  roughness  of  the  ground  specimens  was 

3.5  pm. 

Figure  3  shows  the  strength  behavior  as  a  function  of  quenching 
temperature  difference  for  these  ground  specimens.  Comparisons  of  Figs.  1 
and  3  and  Table  1  show  that  the  grinding  treatment  increased  the  average 
strength  somewhat,  primarily  by  elimination  of  the  low  strength  data  of 
the  non-ground  specimens.  Since  in  unstable  fracture  the  retained 
strength  after  fracture  is  an  inverse  function  of  the  strength  before 
fracture,  removal  of  the  low  strength  specimens  prior  to  the  quench 
also  removed  the  high  strength  specimens  from  the  strength  distributions 
after  failure.  As  a  result,  the  strength  distributions  over  the  range  of 
725°C  £AT<  925°C  now  clearly  indicate  a  pronounced  bi-modal  distribution 
as  expected  from  unstable  fracture  when  only  a  fraction  of  the  specimens 
tested  at  a  given  value  of  AT  undergo  failure. 

The  data  presented  herein  suggest  that  for  materials  which  ex¬ 
hibit  significant  data  scatter,  the  thermal  shock  behavior  should  be 
analyzed  very  carefully  in  order  to  clearly  elucidate  the  mode  of  crack 
propagation.  As  indicated  by  the  present  data,  excessive  coefficients 
of  variation  within  the  range  of  AT  at  which  failure  occurs  may  be  indi¬ 
cative  of  an  unstable  failure  mode  obscured  by  data  scatter.  For  purpose 
of  engineering  design  it  is  critical  that  the  proper  failure  mode  be 
established.  The  situation  could  arise  that  a  material  is  selected  for 
its  apparent  stable  mode  of  failure,  when  in  fact  crack  propagation  will 
occur  in  the  less  desirable  unstable  mode  of  crack  propagation. 
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Strength  and  Standard  Deviation  for  As-nitrided  and  Ground  Reaction-sin tered 
on-Nitride  in  the  As-nitrided  Specimens  and  after  Gil  Quenching 


-p 

c  c 

o  o  x> 

•P  -H  <D 

O  P  TJ 

*P  (3  -H 


2  cr 

X 

X!  “ 

P  >0 

c  a; 
o  n 


P  2  O 
i»  s:  t 


p 

(D 

O'  n 

O 

C  P 

c 

■P  (T3 

<u 

JC  p 

o  a) 

a) 

c  a 

vh 

<u  s 

4-1 

3  O 

O  E-i 

Q 

CN  C-  <N 


o 

VO 

<N 

CO 

o 

CO 

r— 

r- 

CN 

VO 

in  o\ 


00 

LT) 

00 

ro 

CN 

O 

O 

00 

T-‘ 

CO 

t — 

* — 

r- 

r— 

T - 

P<  'f  O' 

O  oo  ■g' 

CN  *-  CN 


CO 

CN 

l 

+  6 

+ 

vo 

5  + 

LO 

+ 

o'  O 

00 

» 

•  , 

l n 

VO 

o  <- 

LO  ld 

T— 

in  r- 

CM 

CN 

CN 

cn 

ro  ro 

tn 

r—  CO 

ro 

T— 

i— 

in  i£>  r^- 


cn  in  cn 
r-  CN  CN 


in 

in 

uo 

lO 

in 

in 

m 

r- 

CM 

r- 

CM 

r~ 

C" 

vo 

r- 

r- 

00 

(T\ 

av 

O 

Fig.  1:  Strength  behaviur  as  a  function  of  quenching  temperature 


Fig.  2 


Fig.  3 


difference  for  the  as-nitrided  specimens  of  reaction-sintered 
silicon  nitride. 

Scanning  electron  micrographs  of  the  as-nitrided  specimens: 
structural  inhomogenit i es  in  the  near-surface  layer  which  cause 
the  relatively  large  scatter  in  strength. 

Strength  behavior  as  a  function  of  quenching  temperature 
difference  for  the  ground  specimens  of  reaction-sintered 
silicon  nitride. 


•  AVERAGE 


Quenching  Temperature  Difference  (°C ) 


Strength  behavior  as  a  function  of  quenching  tempera 
ture  difference  for  the  as-nitrided  specimens  of 
reaction-sintered  silicon  nitride. 
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the  near-surface  layer  which  cause  the  relatively  large  scatter  in  strength. 
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g.  5.  Strength  behavior  as  a  function  of  quenching 
temperature  difference  for  the  surface-ground 
specimens  of  reaction-sintered  silicon  nitride. 
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As  a  partial  solution  to  the  critical  problem  of  the  low  thermal  stress 
resistance  of  many  materials  for  high-temperature  applications,  the  results 
of  two  recent  studies1*2  showed  that  the  magnitude  of  tensile  thermal  stresses 
can  be  reduced  appreciably  by  the  introduction  of  a  spatially  varying  conducti¬ 
vity. 

An  effective  way  of  instroducing  a  spatially  varying  thermal  conductivity 
is  to  introduce  a  spatially  varying  pore  content.  Tnis  should  have  the  added 
advantage  of  lowering  Young's  modulus  as  well.  Furthermore,  if  the  pore 
phase  is  distributed  such  that  full  density  is  maintained  in  the  region  of 
maximum  tensile  stress  within  the  structure  or  component,  the  adverse  effect 
of  porosity  on  tensile  strength  is  avoided.  The  purpose  of  this  study  is 
to  establish  the  quantitative  improvement  in  thermal  stress  resistance  under 
steady-state  heat  flow  due  to  a  spatially  varying  pore  phase. 

The  geometry  selected  for  this  study  consisted  of  an  infinitely  long 
hollow  cylinder  with  inner  radius  a  and  outer  radius  b.  The  cylinder  is 
subjected  to  radially  outward  or  inward  heat  flow,  which  results  in  a  tempera¬ 
ture  difference  AT  across  the  tube  wall. 

For  radially  outward  heat  flow  the  porosity  distribution  was  arbitrarily 
taken  as: 

P  =  P  { (b-r) / (b-a)  f"  (1) 

o 

where  Pq  is  the  porosity  at  r  =  a  and  n  is  the  porosity  distribution  coefficient. 

For  radially  inward  heat  flow  the  porosity  was  assumed  to  be  distributed 
according  to: 

P  =  P  { ( r-a) /  (b-a)  }n 
o 


(2) 


The  effect  of  the  pore  phase  on  the  thermal  conductivity,  K  and  Young's 


modulus,  E  was  taken  as3: 


K,  E 


(3) 


where  the  subscript  refers  to  the  non-porous  material  and  b  is  a  numerical 
constant  taken  as  b  =  3  and  4  for  K  and  E,  resp.  Clearly,  for  a  specific 
material  the  values  of  b  must  be  established  experimentally.  Poisson's 
ratio  and  the  coefficient  of  thermal  expansion  were  assumed  to  be  unaffected 
by  the  pore  phase. 

The  thermal  stresses  were  evaluated  by  the  finite  element  method,  using 
a  computer  program  developed  by  Satyamurthy  et  al‘  .  The  validity  of  this 
numerical  program  was  established  by  verifying  the  analytical  results  of 
Hasselman  and  Youngblood1,  which  gave  excellent  agreement. 

Figures  la  and  lb,  for  radially  outward  and  inward  heat  flow,  resp., 

* 

show  the  values  of  the  non-dimensional  maximum  tensile  stress,  c  =  a(l-v)/aEoAT 
as  a  function  of  the  coefficient  n,  for  three  values  of  b/a  and  =  0.3,  which 
is  a  reasonable  value  achievable  in  practice.  To  a  first  approximation  the 
thermal  stresses  are  independent  of  n.  As  the  result  of  the  spatial  variation 
in  strength  for  value  of  n  <  1  and  a  value  of  b  =  6  (in  eq.  3)  to  express  the 
effect  of  porosity  on  relative  tensile  strength,  failure  can  occur  within  the 
cylinder  away  from  the  region  of  maximum  tensile  thermal  stress.  For  radially 
outward  heat  flow  with  b/a  =  10  and  for  rad i  illv  inward  heat  flow  with  b/a  =  1.1 
the  reduction  in  the  maximum  tensile  thermal  stress  is  about  33  and  22%,  resp. 
Since  the  failure  stress  at  the  site  of  maximum  stress  is  not  affected  by  the 
pore  phase,  these  reductions  in  thermal  stress,  result  in  corresponding  increases 
in  thermal  stress  resistance.  Similarly,  somewhat  smaller  increases  in 


thermal  stress  resistance  are  obtained  lor  the  other  values  of  b/a. 


Figures  2a  and  2b  give  the  magnitude  of  the  maximum  tensile  stress  as 
a  function  of  pore  content,  Pq  for  radially  outward  and  inward  heat  flow, 
resp. ,  corresponding  to  the  value  of  n  for  which  the  stresses  are  minimum 
as  indicated  by  Fig.  1.  These  results  indicate  that  to  a  first  approximation 
the  reduction  in  thermal  stress  is  a  linear  function  of  porosity. 

It  should  be  noted  that  the  improvement  in  thermal  stress  resistance  due 
to  the  presence  of  a  pore  phase,  is  not  expected  from  thermo-elastic  theory 
based  on  spatially  invariant  physical  properties.  In  fact,  a  uniformly 
distributed  pore  phase  is  expected  to  decrease  thermal  stress  resistance1,  > J . 
For  steady-state  heat  flow,  thermal  stress  resistance  as  measured  in  terms 
of  a  maximum  temperature  difference  which  can  be  imposed  on  a  structure  for 
a  given  coefficient  of  thermal  expansion  and  Poisson's  ratio  is  governed  by 
the  ratio  of  tensile  strength  and  Young's  modulus  of  elasticity.  For  a 
uniform  porosity,  P  =  0.3  and  b  =  4  and  6  (eq.  2)  for  Young's  modulus  and 
tensile  strength6,  resp.,  thermal  stress  resistance  would  be  reduced  by 
nearly  a  factor  of  two. 

These  results  indicate  that  for  structures  or  components  which  do 
not  permit  design  modification  or  material  substitution,  the  incidence  of 
thermal  stress  failure  may  be  reduced  by  the  incorporation  of  a  spatially 
varying  pore  phase. 
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Fig.  la.  Maximum  tensile  thermal  stress  in  hollow  cylinder  under  conditions  of 
steady-state  radially  outward  heat-flow  with  spatially  varying  pore 
content  according  to  eq.  1  with  P  (r-a)  =  0.3. 


Fig.  lb.  Maximum  tensile  thermal  stress  in  hollow  cylinder  under  conditions  of 
steady-state  radially  inward  heat  flow  with  spatially  varying  pore 


Maximum  tensile  thermal  stress  in  hollow  cylinder  under  conditions  of 
steady-state  radially  outward  heat  flow  with  spatially  varying  pore 
content  as  a  function  of  maximum  porosity,  P  (r»a). 


Maximum  tensile  thermal  stress  in  hollow  cylinder  under  conditions  of 
steady-state  radially  inward  heat  flow  with  spatially  varying  pore 
content  as  a  function  of  maximum  porosity,  P  (r=b). 
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ABSTRACT 


The  effect  of  interfacial  de-cohesion,  due  to  the  thermal  expansion 
mismatch,  on  the  thermal  diffusivity  of  a  hot-pressed  glass  matrix  with 
a  dispersed  phase  of  nickel  was  investigated  by  the  laser-flash  techni¬ 
que  over  the  temperature  range  of  25  to  600°C.  The  interfacial  gap 
formed  on  cooling  acts  as  a  barrier  to  heat  flow  and  lowers  the  thermal 
diffusivity  to  values  below  those  predicted  from  composite  theory  and 
also  creates  a  temperature  dependence  of  the  thermal  diffusivity  which 
is  strongly  positive.  Preoxidation  of  the  Ni  spheres  promotes  inter¬ 
facial  bonding  to  yield  values  of  thermal  diffusivity  higher  than  those 
for  non-oxidized  spheres  and  a  thermal  diffusivity  which  is  relatively 
temperature  independent.  The  results  of  the  present  study  also  confirm 
the  criteria  for  the  effective  thermal  diffusivity  of  composites  es¬ 
tablished  by  Lee  and  Taylor. 
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EFFECT  OF  THERMAL  EXPANSION  MISMATCH  ON  THE 
THERMAL  DIFFUSIVITY  OF  CLASS-NI  COMPOSITES 

Bob  R.  Powell,  Jr.,  G.  E.  Youngblood,  D.  P.  H.  Hasselman 
and  Larry  I).  Bentsen 

I.  INTRODUCTION 

Numerous  composite  materials  have  been  developed  to  meet  the 
special  requirements  for  many  engineering  applications.  Such  compo¬ 
sites  frequently  exhibit  favorable  properties  not  found  in  the  in¬ 
dividual  components.  The  science  of  "composite  theory"  for  the  cal¬ 
culation  of  the  properties  of  composites  from  the  properties  of  the 
components  represents  a  very  active  field. 

For  the  calculation  of  the  thermal  conductivity  of  composites 
for  a  wide  variety  of  phase  compositions  and  distributions,  many  ex¬ 
pressions  are  available  from  the  1 i tera t ure . 1 ’ ^ ’ 3 ’ 4  These  theories, 
however,  generally  assume  that  perfect  bonding  exists  between  the 
components.  For  this  reason,  the  interfacial  boundaries  offer  no 
resistance  to  the  transport  of  heat  through  the  composite.  In  practice, 
however,  many  composites  will  exhibit  less  than  perfect  bonding  between 
the  components  as  the  direct  result  of  poor  wettability  or  the  absence 
of  good  mechanical  adhesion.  This  problem  is  expected  to  be  especially 
severe  for  composites  with  a  large  mismatch  in  the  coefficients  of 
thermal  expansion  of  the  components.  Changes  in  temperature  of  such 
composites  can  lead  to  the  formation  of  high  values  of  "internal"  stress, 
which  can  lead  to  "micro-cracking"  of  the  components  or  to  interfacial 


separation,  depending  on  the  direction  of  the  mismatch. 


Such  micro- 


cracking  also  can  occur  in  polycrystalline  materials  which  exhibit  a 


high  degree  of  thermal  expansion  anisotropy  of  the  individual  grains. 

In  these  latter  materials  micro-cracking  has  been  shown  to  result  in 

5,6 

major  decreases  in  the  thermal  diffusivity,  in  accordance  with 
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analytical  results  for  the  effect  of  cracks  on  thermal  conductivity. 

Similar  effects  would  be  expected  for  composites  with  a  large  mismatch 
in  thermal  expansion.  However,  as  far  as  these  writers  are  aware  no 
such  experimental  data  exist  in  the  literature. 

The  purpose  of  the  present  study  was  to  investigate  the  effect  of 
the  mismatch  in  thermal  expansion  on  the  thermal  diffusivity  of  a  bonded 
and  non-bonded  composite  consisting  of  a  low  thermal  expansion  matrix 
containing  a  high  thermal  expansion  dispersed  phase.  On  cooling  of 
such  a  composite,  de-cohesion  between  the  phases  results  in  the  formation 
of  a  interfacial  gap,  which  in  analogy  to  micro-cracks  is  expected  to 
have  a  large  effect  on  the  transport  of  heat  through  the  composite. 

II.  EXPERIMENTAL 

A.  Materials,  Preparation  and  Characterization 

D-glass  of  the  same  composition  (16%  Na20,  14%  and  70%  Si07) 

8,9,  10 

are  used  in  previous  investigations  '  was  chosen  as  a  suitable 

matrix  phase.  The  dispersed  phase  consisted  of  near  spherical  particles 

of  Ni  with  a  diameter,  45  ^  D  <  53  pm.  Composites  contained  nickel  in 

the  "raw"  (non-oxidized)  form  or  nickel  particles  with  an  oxide  coating 

produced  by  oxidation  in  air  to  a  weight  gain  of  1  percent.  As  found  in 
10 

an  earlier  study,  the  presence  of  an  oxide  coating  promoted  the  adherence 

-  2  - 

J 


of  the  glass  and  nickel  phase.  The  composites  were  prepared  by  vacuum- 


hot-pressing  mixtures  of  the  Ni  and  glass  in  powder  form  in  graphite  dies 
at  a  pressure  of  14  MPa  at  700°C  for  10  min.  The  initial  heating  rate 
was  9  °C/min.  After  hot-pressing  the  specimens  were  permitted  to 
cool  naturally  within  the  graphite  die  within  the  hot-pressing  chamber 
maintained  at  vacuum. 

The  density  of  the  composites  was  determined  by  the  Archimedes 
principle  with  1^0  as  the  fluid  medium,  from  the  same  specimens  used 
for  the  measurement  of  the  thermal  diffusivitv.  The  experimental  and 
theoretical  densities  are  compared  in  Table  1.  The  sample  containing 
45  vol.%  Ni  exhibits  a  density  greater  than  theoretical.  This  may  be 
due  to  a  Ni  concentrat ion  somewhat  higher  than  45%  or  due  to  uncer¬ 
tainties  in  the  density  of  the  glass.  A  density  less  than  theoretical 
is  attributable  to  the  existence  of  a  pore  phase  created  during  cooling 
from  the  hot-pressing  temperature,  as  the  direct  result  of  the  expansion 
mismatch . 

Figure  1  shows  optical  micrographs  of  the  glass  with  15,  30  and  45 
percent  Ni.  Direct  Ni  partic le-to-part ic le  contact  is  evident  especially 
at  45%  Ni.  Although  most  particles  are  very  close  to  being  spherical, 
some  are  tear-shaped  or  ellipsoidal  in  shape.  As  indicated  by  the 
theoretical  results  of  Bruggeman'1  and  \icsel1,  a  deviation  from  per¬ 
fectly  spherical  shape  should  have  no  effect  on  the  thermal  conductivity 
and  thermal  diffusivitv  as  long  as  the  dispersed  particles  are  oriented 
randomly. 

Figures  2  and  3  show  scanning  electron  f raetographs  for  the  glass 
containing  raw  and  oxidized  Ni  dispersions,  respectively.  It  is  evident 
that  during  fracture  the  raw  Ni  particles  were  easily  pulled  from  the 
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glass  matrix  indicative  of  a  lack  of  adhesion.  For  the  oxidized  Ni 
the  fractographs  indicate  that  the  Ni  and  glass  show  good  adherence. 

In  these  composites  the  thermal  expansion  mismatch  appears  to  be  accom¬ 
modated  by  visco-elastic  deformation  of  the  glass,  which  resulted  in 
the  formation  of  a  highly  porous  fibrous  glass  phase  in  the  immediate 
vicinity  of  the  Ni-glass  interface.  Such  a  fibrous  glass  phase  most 
likely  occurred  at  temperatures  immediately  below  the  hot-pressing 
temperatures  at  which  such  visco-elastic  deformation  most  likely  is 
expected  to  occur. 

B.  Measurement  of  Thermal  Diffusivity 

The  thermal  diffusivity  was  measured  over  the  temperature  range  25°C 
to  600°C  by  the  laser-flash  technique.13  A  single  flash  (a,  1  msec)  from 
a  glass  Nd  laser+  irradiated  the  front  face  of  the  specimen  (1/2  inch 

diameter  by  2.0  mm  thick).  An  intrinsic  thermocouple1'*  (low  temperature) 

* 

or  an  infrared  detector  (high  temperature)  was  used  to  monitor  the  tem¬ 
perature  transient  at  the  back  face.  The  typical  steady-state  temperature 
rise  of  the  specimen  was  kept  below  3°C. 

For  low  temperature  measurements  (25°C-400°C)  the  specimens  were 
held  by  three  alignment  pins  inside  a  wire-wound  temperature  controlled 

split-tube  furnace  with  air  atmosphere.  The  specimen  front  face  was 
** 

coated  with  carbon  to  ensure  uniform  absorptivity.  A  1/4  inch  dia¬ 
meter  Ag-print  spot  was  painted  onto  the  specimen  backface  to  average 
local  heating  effects  and  to  provide  electrical  continuity  for  a  Chromel- 
Constantan  (type  F.)  intrinsic  thermocouple  spring-loaded  to  press  up 

+Korad  K-l ,  70  Joules  maximum  per  flash,  Santa  Monica,  CA. 

*Barnes  Engineering  1T-7B  (InSb)  Stamford,  CT. 

**Miracle  DCF  spray,  Miracle  Power  Products  Co.,  Cleveland,  OH. 

-  4  - 


firmly  against  the  backface  to  monitor  the  chamber  temperature  as  well 
as  the  transient  temperature  of  the  sample. 

At  higher  temperatures (200°C-600°C)  at  which  the  IR  detector  was 
used  for  monitoring  the  specimen  temperature,  the  specimens  with  each 
face  coated  with  carbon  were  held  in  a  split-tube  graphite  holder  in  a 

|  t 

carbon  resistance  furnace  with  nitrogen  atmosphere. 

The  transient  temperature  response  from  either  the  thermocouple  or 
infra-red  detector  was  fed  into  an  oscilloscope  equipped  with  signal 
storage.  The  thermal  diffusivity  (a)  of  the  composite  sample  was  cal¬ 
culated  from 


a  =  0.1388  L2/t1/2  (1) 

where  L  is  the  specimen  thickness  and  t^2  the  time  required  for  the 
temperature  of  the  back  surface  of  the  specimen  to  reach  half  its  final 
value.  The  measured  values  of  the  thermal  diffusivity  were  corrected  for 
heat  losses  in  a  manner  described  by  Heckman.15  These  corrections  amounted 
to  only  a  few  percent  at  low  temperature  and  increased  to  about  5  percent 
at  1000°C.  No  corrections  for  finite  pulse  width  were  necessary.  The 
equipment  was  calibrated  against  Armco  iron  which  resulted  in  agreement 
to  within  3%  of  standard  values16  for  the  infra-red  detector  and  within 
11%  for  the  thermo-couple.  These  discrepancies  could  be  due  to  heat  losses 
along  the  thermo-couple  and  because  above  the  Curie  point  the  thermal 
diffusivity  appears  to  be  a  function  of  thermal  history.15  Also,  intrinsic 
thermocouples  are  known  to  produce  lower  diffusivity  values  than  infra¬ 
red  detectors .  114 

++Astro  Model  1000A,  Santa  Baraba,  CA. 

^Tektronix  5111  with  plug-ins,  Beaverton,  OR. 


In  order  Co  establish  whether  the  value  of  thermal  diffusivity 


obtained  by  the  laser-flash  method  represents  the  effective  thermal 

diffusivity,  the  ratios  of  thermal  properties  and  micro-structural 

parameters  of  the  present  composites  were  compared  with  the  limits  of 

non-homogeneity  established  by  Lee  and  Taylor.17  Briefly,  Lee  and 

Taylor  showed  that  the  concept  of  a  "homogeneous"  thermal  diffusivity 

of  a  composite  is  valid  for  ratios  of  0.48-1137  for  a^/u^;  0.02-1.16 

for  (p.C.VJ/(p  C  V  );  3.8-2857  for  L/D;  0.23-0.88  for  nD/L  and  0.012- 

0.34  for  Vj,  where  a  is  the  thermal  diffusivity,  p  is  the  density,  C 

is  the  specific  heat,  V  is  the  volume  fraction,  L  is  the  sample  thick- 

3  3  1/3 

ness,  D  is  the  particle  diameter,  n  =  [V^L  /1.33n(D/2)  ]  is  the 
total  number  of  particles  in  the  sample  thickness  direction,  and  the 
subscripts  d  and  m  refer  to  the  dispersed  and  the  continuous  phase, 
repsectively .  Table  II  lists  the  thermal  properties  at  room  tempera¬ 
ture  of  the  Ni  and  glass  respectively.  Table  III  lists  the  values  of 
the  above  ratios  for  the  composites  studied. 

It  is  noted  that  all  of  the  various  ratios  fall  within  the  ex¬ 
perimental  bounds  established  by  Lee  and  Taylor  except  nD/L  for  the  45 
v/o  Ni  sample.  Thus,  with  the  possible  expception  of  interfacial  resis¬ 
tance  the  measured  values  of  thermal  diffusivity  represent  the  "homogeneous 
thermal  diffusivity  of  the  composites. 

LI.  RESULTS  AND  DISCUSSION 

The  experimental  data  for  the  thermal  diffusivity  for  D-glass  and 
D-glass  with  15,  30,  and  45  volume  percent  dispersed  nickel  spheres  are 


shown  in  Figures  4  and  5  for  the  glass  with  raw  and  oxidized  nickel 


spheres,  respectively.  The  dashed  curves  in  each  figure  show  the 
thermal  diffusivity  calculated  from  the  Bruggeman*1  variable  dispersion 
(BVD)  relation  for  the  conductivity  of  a  composite  (\^): 


X  1/3 


the  volumetric  heat  capacity: 


(2) 


p  C  =  p  C  V  +  p  ,c,v, 
cc  mrom  ddd 


and  the  relation:  a  =  k/pC 


(4) 


The  implicit  assumption  in  the  derivation  of  eq .  2  is  that  the 
components  are  perfectly  bonded  with  no  thermal  resistance  at  the 
interface.  It  may  also  be  noted  that  the  BVD  relation  gives  a  better 
prediction  of  the  "homogeneous"  thermal  diffusivity  in  the  range  of 
volume  fractions  being  considered  here  than  does  the  Rayleigh  -  Maxwell 
Equation. 1 7 

As  may  be  noted  from  Figure  4  the  thermal  diffusivity,  a,  of  the 
15  and  30  v/o  raw  Ni/glass  composites  from  room  temperature  to  n-  300°C 
matches  the  values  of  thermal  diffusivity  of  100%  D-glass  while  a 
of  the  45  v/o  Ni  composite  is  about  twice  that  of  the  D-glass.  Above 
300°C,  a  has  a  high  positive  temperature  dependence.  Values  of  a  for 
the  15  and  30  v/o  rnw-Ni  composites  increase  with  increasing  temperature 
until  they  show  excellent  agreement  with  values  calculated  by  the  BVD 
relation  at  600°C,  near  the  fabrication  temperature  of  700°C.  For  the 
45  v/o  Ni  composite,  the  thermal  diffusivity  increases  to  values  above 
that  predicted  by  the  BVD  theory  for  temperatures  above  500°C.  The 
Curie  point  transition  at  358°C  for  nickel  is  readily  observable  in  the 


7 


Ni  glass  composites,  although  the  minimum  in  a  appears  displaced  to 
slightly  lower  temperatures. 

For  the  oxidized  Ni/glass  composites  measured  a  values  are  lower 
than  those  calculated  by  the  BVD  relation  for  all  temperatures,  as  shown 
in  Figure  5.  However,  at  room  temperature  a  values  for  the  composites 
with  oxidized  nickel  spheres  are  higher  than  a  values  for  composites 
with  raw  nickel  spheres  of  corresponding  concentration.  The  positive 
temperature  dependence  of  a  for  these  composites  is  much  less  than  those 
with  non-oxidized  spheres.  Again  the  effect  of  the  Curie  transition  on 
thermal  diffusivity  is  observable  at  temperatures  slightly  below  358°C 
for  all  three  nickel  concentrations. 

At  least  qualitatively  the  above  observations  can  be  explained  on 
the  basis  of  the  mismatch  in  the  coefficients  of  thermal  expansion  and 
the  nature  of  and  differences  in  the  interfacial  bonds. 

For  the  glass  with  raw  Ni  spheres,  the  lack  of  bonding  is  clearly 
evident  from  the  SEM  fractographs  shown  in  Fig.  2.  The  lack  of  cohesion 
at  the  interface  between  the  glass  and  Ni  also  can  be  justified  from 
the  magnitude  of  the  internal  stress  which  would  have  arisen  for  per¬ 
fect  bonding.  After  Seising,18  the  magnitude  of  the  radial  stress  (o^) 
perpendicular  to  the  interface  of  a  single  spherical  particle  contained 
in  an  infinite  matrix  is: 

AFATEd 

°r  =  (1+v  )  (E./2E  )+(l-2v  ,)  (3) 

mam  a 

where  AR  is  the  mismatch  in  the  coefficients  of  linear  thermal  expansion, 
AT  is  the  temperature  difference  over  which  the  composite  is  cooled  and 
over  which  no  stress  relaxation  can  occur,  !\  is  Young's  modulus,  v  is 
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Poisson’s  ratio  and,  as  before,  subscripts  d  and  m  refer  to  the  dispersed 

phase  and  matrix  phase,  respectively.  With  Af  =  6.7x10  ^°K  \  AT  =  600°C, 

E  =  6 .0x10 5  MPa,  E  =  20.6xl05  MPa,  .  =  0.3,  v,  =  0.3  vields: 

m  d  m  d 

a  ~  345  MJ’a  (50,000  psi) 

This  magnitude  of  stress  quite  likely  exceeds  the  strength  of  the 
interfacial  bond  of  a  non-wetting  metal -glass  system  such  as  the  pre¬ 
sent,  without  any  bonding  agent.  Separation  of  the  glass-Ni  inter¬ 
face  most  likely  already  occurs  during  the  cooling  phase  immediately 
following  hot-pressing.  It  is  easily  calculated  from  the  thermal  ex¬ 
pansion  mismatch  that  for  a  50  pm  spherical  inclusion  and  a  600oC  tem¬ 
perature  range  of  cooling,  the  width  of  the  interfacial  gap  is  approximately 
0.1  pm.  Such  a  gap  is  expected  to  act  as  a  significant  barrier  to  the 
flow  of  heat.  For  this  reason,  the  measured  values  of  the  thermal 
diffusivitv  are  expected  to  fall  well  helow  those  calculated  from 
composite  theorv  even  at  the  lower  temperatures. 

Despite  the  presence  of  the  gaps,  the  nickel  metal  still  makes  a 
contribution  to  the  thermal  diffusivitv  as  indicated  by  the  unique  tem¬ 
perature  dependence  which  reflects  the  presence  of  the  Curie  point  in 
the  nickel.  Also,  the  thermal  dillusivitv  of  the  glass-raw  Ni  composites 
docs  not  fall  below  the  value  of  the  thermal  diffusivitv  ot  the  glass, 
which  would  have  been  expected  if  the  Ni-phase  made  tin  contribution  to 
the  conduction  of  heat.  The  reason  for  this  is  that  in  the  glass  in  which 
at  the  lower  temperatures  thermal  transport  occurs  primarily  hv  phonon 
transport,  the  thermal  diffusivitv  i  •  ditectlv  proportional  to  the  pro¬ 
duct  of  the  phonon  veloi  itv  and  mean  tree  path  in  the  direct  ion  of  the 
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temperature  gradient.  The  presence  of  the  gaps  increases  the  total 
distance  required  for  phonon  transport  and  accordingly  decreases  the 
thermal  diffusivitv.  Furthermore  if  the  Ni  did  not  contribute  to  the 
conduction  of  heat,  the  thermal  diffusivitv  would  be  expected  to  de¬ 
crease  with  increasing  Ni  content  (i.e.,  increasing  number  of  inter¬ 
facial  gaps).  The  contribution  of  the  Ni  to  the  conduction  of  heat 
can  occur  by  means  of  Ni  part icle-to-practice  contact  especially  ex¬ 
pected  for  the  glass  with  45  vol.%  Ni.  F.ven  isolated  particles  however 
can  contribute  to  the  thermal  transport  since  contact  with  the  glass  can 
occur  at  multiple  points  especially  for  the  non-spherieal  particles. 

With  increases  in  temperature  to  the  hot-pressing  temperature 
the  contact  between  the  glass  and  Ni  inclusions  will  improve.  At  higher 
temperatures,  therefore,  the  Ni  particles  can  make  a  greater  contribution 
to  the  thermal  transport  than  at  the  lower  temperatures.  This  results 
in  a  positive  temperature  dependence  of  the  thermal  diffusivitv,  in 
agreement  with  observations.  For  the  highest  temperatures  (600°C) 
excellent  agreement  exists  between  the  experimental  data  and  those  cal¬ 
culated  by  the  BVD  composite  theory  for  the  glass  with  15  and  30%  Ni. 

This  is  indicative  of  the  absence  of  an  interfacial  barrier  to  heat  flow. 
For  the  composite  with  45%  Ni,  the  observed  values  exceed  those  predicted 
by  the  BVD  theory,  probably  as  the  result  of  an  enhanced  thermal  dif- 
fusivity  due  to  direct  part ic le- to-part id e  contact,  not  taken  into 
account  in  the  BVD  theory.  Tt  should  also  be  noted  that  the  excellent 
agreement  between  the  observed  and  calculated  data  for  the  composites 
with  15  and  30%  Ni  confirms  the  validitv  of  the  criteria  for  the  thermal 
diffusivity  established  by  l.ee  and  Taylor,1'  when  no  significant  inter¬ 
facial  effects  are  present. 
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The  magnitude  and  temperature  dependence  of  the  thermal  diffusivity 
of  the  glass  with  the  oxidized  Ni  dispersion  is  expected  to  show  diffe¬ 
rent  behavior.  Since  the  fibrous  glass  phase  at  the  interface  provides 
physical  contact  between  the  glass  matrix  and  the  nickel,  at  least  at 
the  lower  temperatures  the  latter  can  make  a  much  larger  contribution 
to  the  thermal  diffusivity  than  the  non-oxidized  (raw)  Ni  dispersions. 

These  should  result  in  a  higher  diffusivity  in  agreement  with  observa¬ 
tions.  Because  of  the  presence  of  the  pores  at  the  interface,  the  observed 
values  are  expected  to  fall  below  those  predicted  by  composite  theory,  as 
observed. 

The  fibrous  phase  at  the  interface  will  form  only  at  those  temperatures 
at  which  visco-elastic  deformation  can  take  place.  Below  these  tempera¬ 
tures  on  further  cooling  the  thermal  expansion  mismatch  will  place  the 
fibrous  phase  in  tension  without  further  changes  in  fiber  and  pore  morpho¬ 
logy.  For  this  reason,  the  temperature  dependence  of  the  glass  with  oxidized 
Ni  particles  is  expected  to  be  less  than  for  the  glass  with  raw  Ni  dis¬ 
persions,  again  in  agreement  with  observations. 

The  results  of  the  effect  of  thermal  cycling  (i.e.  heating  followed 
by  cooling)  on  thermal  diffusivity  shown  in  Fig.  6,  provide  additional 
information  on  the  nature  of  the  inter  facial  bond  and  thermal  barrier  at 
the  glass-Ni  interface.  For  the  glass  with  oxidized  N'i  spheres  the  ex¬ 
perimental  data  on  heating  and  cooling  over  one  cycle  were  observed  to 
be  reproducible.  This  suggests  that  the  interfacial  bond  is  relatively 
permanent,  at  least  over  the  temperature  range  of  the  thermal  cycle. 

For  the  glass  with  raw  Ni  dispersions,  however,  the  experimental 
data  on  the  cooling  part  of  the  cvele  lie  well  above  the  data  obtained 
on  heating.  This  suggests  that  near  6Q0°C.  the  glass  and  nickel  form  a 


bond  of  some  sort  which  requires  some  degree  of  cooling  before  rupture 
occurs.  This  results  in  a  discontinuity  in  slope  of  the  thermal  dif- 
fusivity  as  a  function  of  temperature  at  =  400°C.  It  is  of  interest 
to  note  that  i  similar  "hysteresis”  was  observed  by  Siebeneck  et  alr’ 
for  the  effect  of  micro-cracking  on  the  thermal  diffusivity  of  iron 
titanate.  a  lack  of  reproducibility  of  the  "hysteresis"  curve  for 
the  glass-r  w  Ni  composites  in  repeated  thermal  cycling  is  further 
evidence  f  >r  the  lack  of  permanance  of  the  interfacial  bond. 

In  summary,  the  results  of  the  present  study  indicate  that  the 
nature  of  the  interfacial  bond  between  the  individual  components  of 
a  composite  can  have  a  profound  effect  on  the  thermal  transport  pro¬ 
perties  of  such  composites.  Such  effects  should  be  taken  into  account 
in  estimates  of  the  thermal  conductivity  of  composites  which  exhibit 
de-bonding  at  the  interfaces.  The  general  results  of  this  study  also 
suggest  that  controlling  the  degree  of  bonding  at  the  interface  between 
components,  may  constitute  an  effective  mechanism  for  tailoring  the 
heat  conduction  properties  of  composites. 
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Table  I. — Composition  and  Density  of  Glass-Nickel  Composites 


Sample 

p  bulk 

p  theory** 

%theoret ical 

D-glass 

2.476 

•00 

+  15  v/o 

Ni 

(raw) 

3.44 

3.44 

100 

+  30  v/o 

Ni 

(raw) 

4.  29 

4.41 

97.3 

+  45  v/o 

Ni 

(raw) 

5.45 

5.37 

101.5 

+  15  v/o 

Ni 

(ox.  ) 

3.24 

3.43 

94.5 

+  30  v/o 

Ni 

(ox.  ) 

4.24 

4.38 

96.8 

+  45  v/o 

Ni 

(ox.  ) 

5.21 

5.32 

97.9 

**  Estimate  made  using  8.8  and  8.9  gm/cc  for  raw  and  oxidized  Ni, 
respectively. 


Table  II. — Thermophysical  Properties  of  (llass  and  Nickel  at  27°C. 


Property 

Nickel* 

D-glass 

3 

Density  (Kg/m  ) 

8900 

2476 

Specific  Heat  (J/Kg-k) 

450 

+* 

840 

”6  2 

Thermal  Diffusivity  (10  m  /sec) 

22.9 

0.52 

Thermal  Conductivity  (W/m-K) 

91 

1.1 

Linear  Thermal  Expansion  (xlO^)(°K  *) 

13.4 

+* 

7.0 

Literature  values 
+Ref ,  8. 


Optical  micrographs  of  D-glass  containing,  i :  15  vol%;  b:  30 
vol.%  and  c:  45  vol.%  Ni  with  average  dianu  ;  er  of  approx.  50 


ming  electron  fractographs  of  D-glass  with 
.on-bonded  Ni  spheres  and  cavities  (bar  =  3 
'  of  interfacial  gap  between  glass  and  Ni 


Scanning  electron  micrographs  of  D-glass  with  c 
A:  bonded  Ni  spheres  and  cavities  (bar  =  30  uni 
face  (bar  =  2  pm) . 
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ABSTRACT 


By  means  of  the  inclusion  of  a  dispersed  phase  with  low  Young's 
modulus,  the  generally  low  thermal  stress  resistance  of  brittle  materials 
for  high-temperature  structures  can  be  improved  significantly.  An  analy¬ 
sis  of  this  improvement  on  the  basis  of  continuum  and  micro-mechanical 
theory  is  presented  in  this  paper. 

The  low-E  phase  is  shown  to  cause  a  significant  decrease  in  Young's 
modulus,  with  a  negligible  effect  in  the  coefficient  of  thermal  expansion. 

The  relative  change  in  thermal  conductivity  is  a  function  of  the  thermal 
conductivity  of  the  dispersed  phase. 

The  tensile  fracture  stress  is  reduced  significantly,  primarily  due  to 
the  mismatch  in  elastic  properties,  with  smaller  effects  due  to  the  mismatches 
in  the  coefficient  of  thermal  expansion  and  thermal  conductivity. 

The  relative  changes  in  Young's  modulus  and  tensile  fracture  stress 
are  such  as  to  result  in  an  increase  in  the  strain-at-f racture  and  a  simulta¬ 
neous  decrease  in  the  elastic  energy  at  fracture,  the  driving  force  for 
catastrophic  crack  propagation.  The  accompanying  increase  in  fracture 
energy  also  contributes  to  the  improvement  of  thermal  shock  resistance. 

By  changing  the  size  of  the  low-E  inclusion,  trade-offs  can  be  made  between 
strain-at-fractim  and  elastic  energy  at  fracture. 
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ABSTRACT 


The  effect  of  porosity  gradient  on  the  maximum  tensile  thermal 
stresses  during  transient  heat  flow  for  a  circular  cvliinh.r  is  calcu¬ 
lated  using  finite  element  method.  It  is  shown  that  with  appropriate 
porosity  distribution,  the  maximum  tensile  thermal  stress  in  the  cvlind 
is  substantially  reduced  without  the  loss  of  strength  in  the  critical 
region,  thus  improving  the  thermal  stress  resistance  of  the  material. 
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1.  INTRODUCTION 

Brittle  ceramic  materials  are  liighly  susceptible  to  catastrophic 

fracture  due  to  thermal  stresses  which  result  from  constrained  thermal 

expansion^.  A  detailed  analysis  of  the  magnitude  and  distribution  of 

such  stresses  is  essential  for  the  design  of  structures  or  components 

operating  in  non-ambient  thermal  environments.  The  selection  of  materials 

with  optimum  thermal  stress  resistance  can  be  based  on  thermal  stress 

2  3 

resistance  parameters  ’  .  Unfort  unate  1 v ,  even  the  best  material  selected 
on  this  basis,  may  be  unable  to  withstand  the  thermal  stresses  frequently 
encountered  in  practice.  For  such  cases,  new  materials  may  have  to  be 
developed  or  major  design  changes  may  be  called  for. 

For  a  further  solution  of  the  problem  of  thermal  stress  failure,  it 
should  be  noted  that  the  magnitude  of  thermal  stresses  in  a  given  structure 
is  a  function  of  the  imposed  heat  flux,  temperature  difference,  t't  ,  to 
which  the  structure  is  subjected  as  well  as  the  distribution  of  tempera¬ 
ture  within  the  structure.  Generally,  little  or  no  control  exists  over 
the  external  heating  conditions.  However,  by  means  of  appropriate  design 
and/or  materials  selection,  modification  in  the  temperature  distribution 
can  be  made  in  such  a  manne  -  that  the  magnitude  of  the  thermal  stress  is 
reduced.  One  method  by  which  tiiis  can  be  achieved  is  to  incorporate  within 
the  structure,  a  spatially  varying  thermal  conductivity.  The  validity 
of  this  latter  approach  was  substantiated  bv  two  recent  studies  by  Hasselraan 
and  Youngblood^*  and  Satvamurthy  et  al  .  These  studies  showed  that  tor 
both  steady-state  and  transient  convective  heat  transfer  in  a  hollow  and 
solid  cylinder,  cusp. ,  the  magnitude  of  thermal  stress  can  lie  lowered  sub¬ 
stantially  by  a  spatiallv  varying  thermal  rondur' ivitv  distributed  in  such 
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a  manner  that  at  the  hotter  parts  of  the  cylinder  the  thermal  conductivity 
is  lower  than  in  the  other  parts  of  the  cylinder. 

An  effective  method  of  incorporating  a  spatially  varying  thermal 
conductivity  in  a  material  is  by  way  of  a  spatially  varying  pore  phase. 
However,  a  spatially  varying  pore  phase  also  will  cause  spatial  varia¬ 
tions  in  Young's  modulus  and  strength  values.  If  the  pore  phase  were 
distributed  such  that  the  pore  content  is  low  (or  zero)  in  the  region  of 
maximum  tensile  thermal  stress,  to  which  brittle  materials  are  most  sus¬ 
ceptible,  full  advantage  can  be  taken  of  the  spatially  varying  thermal 
conductivity  and  decreased  Young's  modulus,  without  the  accompanying 
adverse  effects  of  porosity  on  tensile  strength.  The  validity  of  this 
concept  was  established  in  a  thermal  stress  analysis  of  a  hollow  cylinder^ 
subjected  to  steady-state  heat  flow.  This  study  showed  that  a  spatially 
varying  pore  phase  with  a  maximum  local  pore  content  of  0.3,  could  improve 
thermal  stress  resistance  by  as  much  as  30%.  This  result  is  particularly 
noteworthy,  since  a  uniformly  distributed  pore  content  of  0.3  can  decrease 
thermal  stress  resistance  by  as  much  as  a  factor  of  two  . 

The  purpose  of  the  present  study  is  to  investigate  the  improvements 
in  thermal  stress  resistance  which  can  be  obtained  by  a  spatially  varying 
pore  phase  under  conditions  of  transient  convective  heat  transfer. 

2 .  PROCEDURES 

2.1  Model  and  Heat  Transfer  Conditions 

A  solid  circular  cylinder  was  selected  as  a  geometry  appropriate  for 
the  present  study.  The  cylinder, initially  at  uniform  temperature,  is  sub¬ 
jected  to  uniform  heating  and  cooling  along  its  surface  by  a  sudden  change 


in  ambient  temperature.  Heat  transfer  is  assumed  to  be  Newtonian  with  a 
heat  transfer  coefficient  h,  described  by  the  Biot  number,  =  ah/K^, 
where  a  is  the  radius  of  the  cylinder  and  Kq  is  the  thermal  conductivity 
of  the  non-porous  material. 

2.2  Porosity  Distribution 

2.2.1  Heating 

The  distribution  of  the  pore  phase  for  convective  heating  was  se¬ 
lected  arbitrarily  as: 

P  =  PQ(r/a)n  (1) 

where  P  is  the  fractional  porosity,  Pq  is  the  maximum  fractional  pore 
content, n  is  a  constant  which  may  be  referred  to  as  the  "porosity  dis¬ 
tribution  coefficient"  and  r  is  the  radial  coordinate.  Eq.  1  indicates 
that  at  the  center  of  the  cylinder  (r=0) ,  where  the  tensile  thermal 

stresses  are  maximum,  the  pore  content  is  zero,  so  that  the  tensile 

strength  at  r=0  is  not  affected. 

2.2.2  Cooling 

Preliminary  calculations  showed  that  a  continuous  function  of  the 
general  form  of  eq.  1  would  not  give  significant  reduction  in  the  thermal 
stresses.  Instead,  the  porosity  distribution  was  chosen  as 

P  -  0  for  0.9  ^  r/a  -  1  (2) 

P  =  PQU  -  (r/0. 9a)°  !  for  0  •:  r/a  0.9  (3) 

Eqs.  2  and  3  indicate  that  the  surface,  where  the  tensile  stresses 
are  maximum,  has  a  non-porous  layer  of  thickness  0.1  r/a. 
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2.3  Effect  of  Porosity  on  Material  Properties 


The  dependence  of  the  thermal  conductivity,  K,  Young's  modulus, 

E  and  the  tensile  strength,  ,on  the  porosity  was  chosen  arbitrarily 
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to  be  expressed  by  the  well-known  exponential  equation  ’  : 


K,  E,  a,  =  K  ,  E  ,  a,  e 
f  o’  o’  fo 


-bP 


(4) 


where  the  zero  subscript  refers  to  the  value  of  the  non-porous  material 
and  b  is  a  constant  taken  as  3,  4  and  6  for  the  thermal  conductivity, 
Young's  modulus  and  tensile  strength  resp.  Poisson's  ratio  and  the  co¬ 
efficient  of  thermal  expansion  were  assumed  to  be  unaffected  by  the  pore 
phase . 


2.4  Computational  Method 

Due  to  the  spatial  variation  of  the  thermal  conductivity  and  Young's 
modulus,  the  transient  temperatures  and  thermal  stresses  required  evalua¬ 
tion  by  numerical  method.  For  this  purpose  a  computer  program  based  on 
finite  element  principles  was  utilized.  This  program  was  developed  and 
used  in  previous  studies"*’^,  and  gave  excellent  agreement  with  the  analytical 
solutions  of  Jaeger^0  for  the  transient  thermal  stresses  in  a  solid  circular 
cylinder  with  uniform  thermal  conductivity.  The  program  also  was  used  to 

it 

verify  the  analytical  results  obtained  by  Hasseltnan  and  Youngblood  for 
the  steady-state  thermal  stress  in  a  hollow  cylinder  with  spatially  varying 
thermal  conductivity. 

For  convenience,  the  numerical  results  obtained  in  the  present  study 

A 

were  reported  in  terms  of  the  non-dimensional  stress,  a  =  o(l-v)/aEo&T,  the 
non-dimensional  radial  coordinate  r/a  and  the  Biot  number,  ah/KQ,  defined 
earlier.  Since  the  calculations  of  the  transient  thermal  stresses  required 
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extensive  computer  analysis,  only  sufficient  data  points  were  generated 
to  indicate  the  positive  effect  of  the  spatially  varying  pore  phase  on 
thermal  stress  resistance  in  general,  rather  than  evaluating  thermal 
stresses  for  specific  applications. 

3.  NUMERICAL  RESULTS  AND  DISCUSSION 

3.1  Heating 

Figures  la  and  lb  show  the  spatial  variation  in  the  transient  tem¬ 
peratures  and  thermal  stresses  for  SQ  =  1,  5  and  20,  n  =  1  and  Pq  =  0.3, 
at  the  instant  of  time  the  tensile  thermal  stresses  at  the  center  of  the 
cylinder  reach  their  maximum  value.  A  value  of  Pq  =  0.3  corresponds  to 
a  value  realizable  in  practice.  Included  in  figs,  la  and  lb  are  the 
temperature  and  thermal  stress  distributions  for  the  fully  dense  material, 
Pq  =  0.  Comparisons  of  these  figures  show  that  the  spatially  varying  pore 
phase  causes  a  significant  change  in  the  temperature  distribution  and  a 
corresponding  decrease  in  the  maximum  tensile  stress  at  r  =  0,  in  agree¬ 
ment  with  the  original  hypothesis  on  which  this  study  was  based. 

Figure  2  for  Pq  =  0.3,  shows  the  values  of  the  maximum  tensile  thermal 

stress  at  r  =  0,  for  6  =  1,  5  and  20  as  a  function  of  the  porosity  dis- 

o 

tribution  function,  n.  For  nil,  the  magnitude  of  the  thermal  stress  is 
nearly  independent  of  n,  i.e.,  the  manner  in  which  the  porosity  is  dis¬ 
tributed.  A  sensitivity  analysis  showed  that  this  independence  of  stress 

on  n  is  due  to  the  counteracting  contributions  of  the  effect  of  the  spatial 

i 

variation  in  the  thermal  conductivity  and  Young's  modulus  on  the  magnitude 
of  maxinum  tensile  thermal  stress. 

Included  in  fig.  2  are  the  values  of  the  maximum  tensile  thermal 
stress  for  the  non-porous  material,  (P  =  0).  Comparison  of  these  data 
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show  that  the  spatially  varying  pore  phase  causes  a  significant  decrease 

in  the  magnitude  of  the  tensile  thermal  stress  at  r  =  0.  Since  at  this 

position  (P  =0),  the  tensile  strength  of  the  material  is  not  affected, 

the  decrease  in  tensile  thermal  stress  leads  to  a  corresponding  increase 

in  thermal  stress  resistance.  This  relative  increase  in  thermal  stress 

resistance  due  to  the  spatially  varying  pore  phase  is  shown  in  fig.  3  for 

3  =  1,  5  and  20  and  P  =  0.3.  These  increases  are  considerable.  It 

o  o 

is,  however,  important  to  note  in  figure  3  that  due  to  rapid  decrease  in 
the  tensile  strength  close  to  the  center  of  the  cylinder,  failure  will 
not  initiate  at  the  center  where  the  tensile  stress  is  maximum  but  away 
from  it  where  the  strength  to  stress  ratio  is  lower.  In  view  of  this  ob¬ 
servation  a  porosity  distribution  which  includes  a  small  core  (r/a  =  0.1) 
of  fully  dense  material  can  be  used  to  further  advantage.  This  will  not 
cause  any  significant  change  in  stress  distribution  because  of  little 
change  in  temperature  distribution  and  Young's  modulus  values.  However, 
the  strength  distribution  will  be  improved  substantially  in  the  critical  re¬ 
gion  resulting  in  a  further  gain  in  thermal  stress  resistance.  Additional 
calculations  showed  that  the  relative  increase  which  can  be  obtained  by 
a  spatially  varying  porosity  described  by  equation  1,  is  approximately 
a  linear  function  of  P  ,  as  also  found  in  a  previous  study  for  steady- 
state  heat  flow^. 

It  should  be  noted  that  the  relative  increase  in  thermal  stress  re¬ 
sistance  shown  in  figure  3,  is  reported  on  the  basis  of  constant  Biot 
number,  3Q  =  ah/KQ,  where  Kq  is  the  thermal  conductivity  of  the  non- 
porous  material.  However,  the  porosity  in  the  surface  lowers  the  thermal 
conductivity  considerably,  which  in  effect  increases  the  Biot  number. 

This  increase  in  the  Biot  number  for  spatially  varying  material  properties 
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also  Increases  the  magnitude  of  thermal  stress.  In  fact,  if  this  change 
in  Biot  number  had  been  taken  into  account  as  well,  the  relative  increase 
in  thermal  stress  resistance  obtainable  for  the  spatially  varying  pore 
phase,  would  be  even  greater.  For  this  reason,  the  data  reported  in  fig.  3 
are  based  on  the  roost  conservative  comparison. 

3.2  Cooling 

On  convective  cooling  the  tensile  thermal  stresses  are  maximum  at 
the  surface  of  the  cylinder,  (r=a).  Preliminary  calculations  indicated 
that  for  a  continuous  porosity  distribution  of  the  form: 

P  =  PQ{1  -  (r/a)n)  (5) 

no  significant  improvements  in  thermal  stress  resistance  could  be  ob¬ 
tained.  Furthermore  for  n  <  1,  the  tensile  strength  close  to  the  surface 
would  decrease  very  rapidly.  Under  these  conditions  similar  to  heating 
case,  failure  would  not  be  initiated  in  the  surface  under  the  in¬ 
fluence  of  the  maximum  tensile  stress,  but  within  the  interior  of  the 
cylinder  with  a  lower  tensile  strength  to  stress  ratio.  This  latter 
effect,  justified  selecting  a  porosity  distribution  function  of  the  type 
given  by  eqs.  2  and  3  which  avoids  the  spatially  rapid  decrease  in  tensile 
strength  in  the  immediate  surface  area. 

Figures  4a  and  4b  show  the  spatial  variation  of  transient  temperature 

and  thermal  stresses  for  n  =  1,  P_  =  3  and  for  8  =1,5  and  20.  The 

o  o 

data  for  the  case  of  constant  conductivity  (P  =  0)  overlap  with  the  ones 
for  variable  porosity  case  and  hence  have  not  been  shown  for  the  purpose 
of  clarity.  However,  fig.  5  shows  the  plots  for  thermal  stresses  relative 
to  the  maximum  stress  for  the  zero  porosity  case.  Included  in  fig.  5  is 
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also  the  strength  distribution.  It  can  be  seen  that  the  increase  in  thermal 


stress  resistance  is  less  than  that  obtainable  for  the  heating  case. 

3.3  General 

The  data  presented  in  the  present  paper  demonstrate  that  a  spatially 
varying  pore  phase  can  be  used  effectively  to  improve  the  thermal  stress 
resistance  of  brittle  ceramics.  For  the  solid  cylinder  chosen  for  the 
present  study,  such  improvements  are  particularly  significant  for  thermal 
shock  by  convective  heating.  For  other  geometries,  heat  transfer  environ¬ 
ments  and  porosity  distribution  the  improvements  in  thermal  stress  resis¬ 
tance  which  can  be  obtained  must  be  investigated  on  an  individual  basis. 

The  results  presented  in  this  paper  suggest  that  this  method  can  be  used 
to  advantage. 
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NON-DIMENSIONAL  TEMPERATURE, 


RADIAL  DISTANCE  ,  r/a 

Transient  thermal  stress  distribution  in  solid  circular  cylinde 
subjected  to  sudden  convect ive  heating  at  the  time  of  maximum 
tensile  stress  for  spatially  varying  porosity  described  by 
equation  1  with  n  =  1  and  PD  =  0.1. 
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Maximum  tensile  thermal  stress  in  solid  circular  cylinder  subjected  to  sudden  convective  heat¬ 
ing  for  spatially  varying  porosity  described  by  equation  1  with  PQ  =  0.3  for  a  range  of  values 


RADIAL  DISTANCE  ( r/a) 

Distribution  of  thermal  stress  and  tensile  strength  for  Po=0.3  and 
n=l  relative  to  the  corresponding  maximum  values  at  Po=0  in  a  solid 
circular  cylinder  subjected  to  sudden  convective  heating  for  spatial 
ly  varying  porosity  described  by  equation  1. 


NON-DIMENSIONAL  TEMPERATURE, 


RADIAL  DISTANCE, r/a 

Fig.  4a.  Transient  temperature  distribution  in  a  solid  rireular  cylinder  sub¬ 
jected  to  sudden  convective  cooling  at  tbe  time  of  maximum  tensile 
stress  for  spatially  varying  porosity  described  by  equations  2  and  3 
with  n=l  and  Po=0.3. 


NON-DIMENSIONAL  STRESS, 


RADIAL  DISTANCE,  r/o 


Fig.  4b.  Transient  thermal  stress  distribution  in  a  solid  circular  cylinder  sub¬ 
jected  to  sudden  convective  cooling  at  the  time  of  maximum  tensile 
stress  for  spatially  varying  porosity  described  by  equations  2  and  3 
with  n=l  and  Po=0.3. 
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Distribution  of  thermal  stress  and  tensile  strength  for  Po=0.3  and  n=i  relative  to  the 
corresponding  maximum  values  at  Pq“0  in  a  solid  circular  cylinder  subjected  to  sudden 
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AESTRACT 


An  analysis  is  presented  of  the  thermal  stresses  in  a  partially 
absorbing  brittle  ceramic  flat  plate  assynmetrically  heated  by  radiation 
on  the  front  surface  and  cooled  by  convection  on  the  rear  surface  with  a 
heat  transfer  coefficient,  h,  being  finite  or  zero. 

For  finite  h,  the  transient  thermal  stresses  are  a  function  of  h, 
whereas  the  steady-state  thermal  stresses  are  independent  of  h.  The 
maximum  value  of  tensile  thermal  stresses  occurs  at  an  optical  thickness 
pa  ■  2  and  equals  zero  for  pa  •=  0  or  °°.  For  an  optical  thickness  pa  <  10 
the  steady-state  thermal  stresses  exceed  the  transient  stresses,  with  the 
converse  being  true  for  pa  >  10.7.  The  maximum  tensile  thermal  stresses 
occur  in  the  front  surface  where  the  temperatures  are  highest. 

For  h  *  0,  the  tensile  thermal  stresses  increase  with  increasing 
value  of  optical  thickness.  The  role  of  the  material  properties  which 
control  the  thermal  stress  resistance  under  conditions  of  combined  radia¬ 
tion  heating  and  convection  cooling  are  discussed.  Appropriate  thermal 
stress  resistance  parameters  are  proposed. 


I.  Introduction 


Brittle  ceramics  are  excellent  candidate  materials  of  construction 
for  solar  energy  collectors  involving  concentrated  radiation  and  other 
applications  such  as  found  in  aerospace  technology,  subject  to  high- 
intensity  radiant  heat  fluxes.  Unfortunately,  in  part  due  to  high  brit¬ 
tleness  and  an  unfavorable  combination  of  other  pertinent  properties, 
ceramic  materials  are  susceptible  to  catastrophic  failure  under  the  in¬ 
fluence  of  the  high  thermal  stress  which  inevitably  accompany  transient 
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or  steady  state  heat  conduction.  *  For  the  purpose  of  selection  of  the 
optimum  material  of  construction  and  reliable  engineering  design,  it  is 
imperative  that  the  variables  which  control  the  magnitude  of  thermal 
stresses  arising  from  radiation  heating  are  well  understood. 

An  earlier  analysis  of  thermal  stresses  by  radiation  heating  con- 

3 

sidered  materials  opaque  to  incident  black-body  radiation.  This  analysis 

was  extended  to  include  materials  opaque  above  a  given  wavelength  of  the 

black-body  radiation  and  completely  transparent  below  this  value  of  wave- 
4  5 

length.  A  more  recent  analysis  concentrated  on  partially  absorbing 
ceramics  symmetrically  heated  by  radiation  which  is  absorbed  as  it  is 
being  transmitted  through  the  material.  This  analysis  showed  that  for 
ceramic  plates  with  small  optical  thickness,  the  thermal  stresses  are 
proportional  to  the  first  or  third  power  of  the  absorption  coefficient, 
for  convective  cooling  corresponding  to  heat  transfer  coefficients  h  =  0 
and  h  -  »,  respectively. 

The  purpose  of  the  present  paper  is  to  analyze  the  role  of  the  perti¬ 
nent  material  properties  which  control  the  magnitude  of  thermal  stresses 
in  a  partially  absorbing  ceramic  plate  subjected  to  thermal  radiation  on 
one  side  and  cooled  by  convection  with  values  of  the  heat  transfer  coef¬ 
ficient  h  ranging  from  zero  to  infinity  at  the  opposite  side. 


II.  Analysis 


A.  Geometry  and  boundary  conditions 

The  geometry  selected  for  the  analysis  consisted  of  an  unconstrained 
infinite  flat  plate  of  thickness  2a,  located  parallel  to  the  y,z-plane 
and  -a  <  x  <  a.  The  plate,  initially  at  uniform  temperature  T  ,  is  sub¬ 
jected  to  a  constant  normally  Incident  uniform  thermal  radiation  atx  *  -a 
beginning  at  time  t  ■  0.  At  x  -  a,  the  plate  is  cooled  by  Newtonian  con¬ 
vection  with  a  heat  transfer  coefficient  h.  The  optical  properties  of 
the  plate  were  assumed  to  be  independent  of  wavelength.  Also,  the  reflec¬ 
tivity  of  the  plate  was  assumed  to  be  sufficiently  low  so  that  the  effect 
of  multiple  reflections  within  the  plate  can  be  neglected.  The  other 
relevant  material  properties  which  control  the  magnitude  of  thermal 
stresses,  such  as  the  coefficient  of  thermal  expansion,  Young's  modulus, 
Poisson's  ratio,  the  thermal  conductivity  and  thermal  diffusivity  were 
assumed  to  be  independent  of  temperature. 
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Following  the  approach  of  the  previous  analyses  ’  it  will  be 

assumed  that  over  the  time  period  required  for  the  stresses  to  reach  their 

maximum  value,  the  temperature  of  the  plate  remains  sufficiently  low  so  that 

re-emission  of  the  absorbed  energy  is  negligible  in  comparison  to  the 

incident  radiation.  The  validity  of  this  assumption  was  demonstrated 

3  5 

in  the  earlier  analyses  ’  by  means  of  a  numerical  example,  as  will  also 
be  done  in  the  present  study. 

B.  Derivation  of  transient  temperatures 

For  normally  incident  radiation  of  Intensity  qQ  on  the  surface  of 

the  plate,  the  flux  transmitted  through  the  surface  equals  eq  ,  where  e 

o 

is  the  emlssivity  or  e  ■  1  -  r,  where  r  is  the  reflectivity. 
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The  intensity  q  of  the  radiation  within  the  plate  at  a  distance  x  Is: 


q  =  eqo  e 


■u(a  +  x) 


(1) 


where  y  is  the  absorption  coefficient. 

The  rate  of  heat  absorption  (g ' *  * )  at  x  is: 


tii  -y (a  +  x) 

g  =  yeqQe 


(2) 


For  the  appropriate  boundary  conditions  the  solutions  for  the  temperature 


(T)  are  obtained  by  solving  the  heat  conduction  equation  : 

32T/3  x2  +  g"'(x)/k  =  (1/k)  3T/3x 


(3) 


where  k  is  the  thermal  conductivity  and  k  is  the  thermal  diffusivity. 

Two  cooling  conditions  of  the  plate  at  x  -  a  were  considered,  namely 
a  heat  transfer  coefficient  0  <  h  £  °°  and  h  =  0.  In  this  respect,  it 
should  be  noted  that  a  finite  heat  transfer  coefficient  h  -*•  0,  does  not 
correspond  to  the  same  condition  obtained  for  h  =  0. 

The  solutions  for  the  transient  temperatures  are  as  follows: 


a.  Finite  heat  transfer  coefficient,  h. 

For  this  case  the  initial  and  boundary  conditions  are: 

T(x,0)  =  T  ;  3T/3x(-a , t)  =  0 
o 

-k  3T/3x(a, t)  =  h{T(a,t)  -  T  \ 

o 

Solution  of  Eq.  (3)  for  these  boundary  conditions  yields: 

00 

T(x,t)  ■  T  +  B  (t)  cosA  (x  +  a) 

o  c _ i  n  n 


n  =  0 


-<X2t. 


where 


B  (t)  =  G  [l  -  e  "  ]M2 
n  nt  J  n 


„  _  Kq0^ 

^n  N  k 
n 


y  +  e  sin(2A  a)  -  ucos(2X  a)] 

_ n _ n _ n 

2  .  .2 
y  +  X 

n 


(4a) 

(4b) 

(3a) 

(3b) 

(5c) 
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N  =  t— ~[sin(2A  a)  cos(2A  a)  +  2A  a]  ( 

n  2A  n  n  n 

n 

and  X  are  the  roots  of  the  transcendental  equation: 
n 

A  tan(2A  a)  =  h/k  ( 

n  n 

A  value  of  heat  transfer  coefficient  h  =  <=°  represents  the  limiting 
case  for  the  finite  value  of  heat  transfer  coefficient. 

For  h  *  00 ,  the  initial  and  boundary  conditions  are: 


T(x,0)  =  T(a,t)  «  Tq;  9T/9x(-a,t)  =  0 


Solution  of  Eq.  (3)  for  these  boundary  conditions  yields: 

oo  .  /nit  i  2 

Z16q  epa  u+  h— 1  “2Pa  -*X  t  , 

— — — —  — - -  1  -  e  )  cos{ A  (x  +  a)}  (7a) 

n  =  odd  n  *  k  u  +(nirMa) 


where 


A  *  nir/4a  with  n  =  1,  3,  5,  7  .  .  . 
n 


b.  Heat  transfer  coefficient  h  *  0. 


The  initial  and  boundary  conditions  are: 


T(x,0)  -  T 


3T/9x(-a, t)  -  3T/3x(a,t)  =  0 


Solution  of  Eq.  (3)  for  the  boundary  conditions  of  Eq.  (8)  yields: 


-xA2t 


C  e  n  cos{ A  (x  +  a) } 
n  n 


n  *  0 


q  e  r  -p(x  +  a)  ,  ,  .2  0  n 

o  e  (x  +  a)  ,,  _-2yaN  .  ,  .  , 

_  _ - - - - -  (1  -  e  )  +  (x  +  a) 

.  « 


tcq  e  - 

+  ST  (1  -  e  3)t 


where: 


V  ri  -  (-l)ne~2tJa  (-l)n 


n  ak 


~1  -  (-D 
2  . 

L  w  + 


(i .  e-2“a>  + 


n  +  0  (10a) 


(10b) 


and 


C 

o 


V 

2ak 


'(1  - 

+  2a2 

W  i 

■ 

A  =  mr/2a  with  n  =  0,  1,  2,  3  .  . 
n 


(10c) 


C.  Derivation  of  the  thermal  stresses 


From  the  known  temperature  distributions,  the  thermal  stresses  can 
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a  a 


be  obtained  from:' 


a  =  [  -T  +  ~  f  Tdx  +  — _  f 

y,z  1  -  v  |  2a  J  ^3  J 

-a  -a 


Txdx 


(ID 


where  a  is  the  coefficient  of  thermal  expansion,  E  is  Young's  modulus  and 
v  is  Poisson's  ratio. 

For  the  two  cooling  conditions  the  stresses  are  as  follows: 


a.  Finite  heat  transfer  coefficient  h. 

For  finite  values  of  the  heat  transfer  coefficient,  from  Eqs.  (5) 
and  (11),  the  thermal  stresses  can  be  derived  to  be: 

oo 

o  =  -  —  V  B  (t)  cos{A  (x  +  a)  } 

y,z  1  -  v  L_,  n  n 

n  =  0 


+ 


ctE 

2a(l  -  v) 


sin(2A  a) 
n 


3xaE 


2a(l  -  v) 


z 

n  =  0 


B  (t) 
n 


sin(2A  a) 
_ n 

aA 


■t-t-  cos  (2A  a) 
2,2  n 


a  A 


2  2 
a  A 

n 


(12) 


with  B^(t)  and  A^  defined  by  Eq.  (5). 

For  the  special  case  of  h  =  ®  from  Eqs.  (7)  and  (11),  the  thermal 


stresses  can  be  calculated  to  be: 


z 


aE 


1  -  v 


I 


odd 


16q  epa  .  ,  .  -2pa 

o  V  +  (nu/4a)e 

2  2  2  2 
n  it  k  p  +  (nn/4a) 


-kA2t 

e  n  I  cos{ A  (x  +  a) } 
I  n 


aE 


8q  £p 


1  -  v 


V"1  H°  U  +  (mr/4a) 

/  2  2  2 

n  ii  k  p  +  (nr/4a) 

n  =  odd  *• 


-2pa 


,2.  n  -  1 

-K)  t  — - — 

r  I1  - e  "j*-1* 


3xaE 


8q  ep 
o 

(1  -  v)  *  n27r2k 

n  =  odd 


E 


p  +  (nTi/4a)e 


■2pa 


p2  +  (nn/4a)2 
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The  stresses  attain  their  maximum  value  at  t  =  For  the  maximum 
tensile  stresses  which  occur  at  x  =  -a,  Eq.  (13)  can  be  simplified  for  the 
optically  thick  (pa  -*•  °°)  and  the  optically  thin  (pa<<l)  plates  as  follows: 


o  (-a,“)  =  0 

y,z 

(pa  ® ) 

(14a) 

0.2719aEq  epa2 

o  (— a,00)  =  - r - 2 - 

y,z  (1  -  v)k 

(pa<<l) 

(14b) 

Setting  o  in  Eq.  (14)  equal  to  the  tensile  strength, 

y  * z 

St>  yields 

the  maximum  possible  heat  flux,  q  ,  to  which  the  plate  can  be  subjected 

max  J 

without  risk  of  failure  under  steady-state  conditions  (t  -+  °°) : 


’max 


3.677S  (1  -  v)k 


max 


(pa  -*•  «) 


(pa<<l) 


aEepa 

b.  Heat  transfer  coefficient  h  =  0. 

For  h  ■  0,  the  thermal  stresses  from  Eqs.  (9)  and  (11)  can  be  derived 


(15a) 

(15b) 


to  be: 
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=  V  C  e"V  coS(A  (x  +  a) }  - 

,Z  1  -  vl  L  n  n 


V  [  e~u(~*  +  a) 


n  =  1 


(1  -  e"2v,a)  +  (x  +  a) 
4a 


aE  [  V  -  -2ua  ,  Va  .  -2tja.  %ea 

+  r tt  <e  - 1> +  —  (1  - e  >  -  — 

L  2V  ak 


2(1  -  v)a' 


-<x2t 

P  2C  e  n  q  £  r  1 

-  2_,  — ^~2 —  +  ~ir  e  pa(,ja  +  i)  +  (va  - 1) 

n  =  odd  Xn  P  *  L  J 


3  3  1 

q  ea  „  2q  ea 

+  ir-a-e  >--5ir-|  (1< 

Eq.  (16)  can  be  simplified  for  an  optically  thick  plate  (ya  ■+•  °°)  and 
an  optically  thin  plate  (ya<<l).  The  maximum  tensile  stresses,  which  occur 
at  t  =  are  of  primary  interest. 

For  the  optically  thick  plate  the  tensile  thermal  stresses  have  their 
maximum  value  at  x  =  0  and  are  given  by: 


a  (0,°°)  =  aEq  ea/12(l  -  v)k 
y,z  o 


(ya  -*■  °°)  (17a) 


For  the  optically  thin  plate  the  tensile  thermal  stresses  reach  their 


maximum  value  at  x  =  -a  and  are  given  by: 


o  (-a,°°)  =  aEq  cya  /(I  -  v)k 
y,z  o 


(ya<<l) 


(17b) 


In  analogy  to  Eqs.  (15a)  and  (13b),  Eqs.  (17a)  and  (17b),  setting 
the  maximum  tensile  stress,  o  equal  to  the  tensile  strength  yields 

y.z 

an  expression  for  the  maximum  permissible  radiation  heat  flux  to  which 
the  plate  can  be  subjected  without  risk  of  failure. 


12S  (1  -  v)k 
t 


St  (1  -  v)k 


(ya  -*•  »)  (18a) 


(y a<<l) 


(18b) 


aEcya 
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III.  Numerical  results  and  discussion 


The  temperatures,  thermal  stresses  and  their  time  dependence  for  con¬ 
venience  are  reported  in  terms  of  the  non-dimensional  temperature: 

T*  =  (T  -  T  )k/eq  a  (19) 

o  o 

the  non-dimensional  thermal  stress: 

a*  =  a(l  -  v)k/aEeq  a  (20) 

o 

and  the  non-dimensional  time: 

t*  =  Kt/a2  (21) 

For  simplicity,  the  quantities  T*,  a*,  and  t*  hereafter  will  be  referred 
to  as  temperature,  stress,  and  time,  respectively.  Because  brittle  ceramics 
generally  fail  in  tension,  the  maximum  tensile  stresses  are  of  primary 
interest. 

The  validity  of  the  assumption  that  the  stresses  reach  a  value  near 
or  equal  to  their  maximum  value  before  the  plate  becomes  hot  enough  that 
the  re-emitted  radiant  energy  must  be  taken  into  account,  must  be  examined 
first.  Most  conveniently,  this  can  be  done  for  the  extremes  in  the 
values  of  the  heat  transfer  coefficient,  i.e.,  h  =  0  and  h  =  °°. 

For  h  =  0,  Figs,  la  and  lb  show  the  maximum  temperatures,  which 

occur  at  x  =  -a,  and  maximum  tensile  thermal  stresses,  respectively,  as 

a  function  of  time,  for  several  values  of  the  optical  thickness  ua.  The 
corresponding  values  for  h  =  00  are  shown  in  Figs.  2a  and  2b.  Table  I 
lists  the  calculated  surface  temperatures  for  a  value  of  incident  heat 
flux,  9max>  for  a  material  with  properties  approximately  corresponding 
to  those  of  aluminum  oxide.  As  indicated  by  the  data  in  Table  I,  the 
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re-emitted  heat  flux  is  negligible  compared  to  the  incident  heat  flux 
at  a  time  when  the  stresses  are  within  5  percent  of  their  maximum  value. 

Examination  of  Figs,  la  and  lb  show  that  for  h  =  0,  after  an  initial 
transient,  during  which  the  stress  reaches  a  maximum,  the  temperature 
increases  linearly  with  time,  because  of  the  absence  of  cooling.  In  con¬ 
trast,  for  h  =  °°  as  shown  in  Figs.  2a  and  2b,  both  the  temperatures  and 
stresses  after  an  initial  transient  reach  steady-state  values  as  t* 

A  number  of  interesting  features  of  the  time  dependence  and  the 
magnitude  of  the  tensile  stresses  may  be  noted  from  Fig.  2b.  First,  the 
stresses  increase  with  time  to  pass  through  a  maximum  followed  by  a 
decrease  and  then  to  rise  again  to  a  steady-state  value  as  t  -+  00 .  This 
behavior  is  the  direct  result  of  a  corresponding  time  dependence  of  the 
position  within  the  plate  at  which  the  tensile  stresses  are  a  maximum. 

Figure  2b  also  shows  that  the  steady-state  stresses  (t  -*  <*>)  increase 
with  increasing  value  of  optical  thickness  to  reach  a  maximum  at  pa  =  2, 
followed  by  a  decrease  as  pa  -+  °°.  In  fact,  it  can  be  shown  that  for 
pa  =  °°,  the  steady-state  thermal  stresses  are  identically  equal  to  zero. 
Furthermore,  it  may  be  noted  that  the  steady-state  stresses  exceed  the 
transient  stresses  or  vice-versa,  depending  on  the  value  of  pa. 

Figure  3  shows  the  time  dependence  of  the  stresses  for  pa  =  3  and  a 
range  of  values  of  the  heat  transfer  coefficient.  It  can  be  noted  that 
after  an  initial  transient,  the  maximum  tensile  stress  is  independent  of 
the  value  of  h.  As  may  be  determined  from  the  temperature  solutions,  the 
final  steady-state  temperatures  at  which  the  h-independent  stress  value 
occurs  is  inversely  proportional  to  the  heat  transfer  coefficient.  With 
the  exception  of  h  =  0,  these  observations  and  conclusions  can  be  shown 
to  be  valid  for  any  value  of  pa. 
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Figure  4  summarizes  the  values  of  maximum  tensile  thermal  stress  as 
a  function  of  optical  thickness.  For  finite  values  of  h  or  h  =  »,  the 
steady-state  stresses  (t  ■+  <°)  exceed  the  transient  stresses  for  0  <  pa  <  '10.7. 
For  10.7  <  pa  <  °°,  the  maximum  transient  stresses  exceed  the  steady-state 
stresses.  As  discussed  earlier,  for  h  =  <=°  (or  finite)  the  stresses  pass 
through  a  maximum  at  pa  -  2  and  vanish  for  pa  =  0,°°.  The  data  in  Fig.  4 

qualitatively  are  quite  similar  to  the  corresponding  data  obtained  in  the 

earlier  study^  for  a  symmetrically  heated  and  cooled  plate.  The  magni¬ 
tude  of  maximum  tensile  stress  found  in  the  present  study  for  the  assym- 
metrically  heated  and  cooled  flat  plate  is  less  than  that  for  the  sym¬ 
metrically  heated  and  cooled  plate.  At  least  one  reason  for  this  difference 
is  that  in  the  assymmetrically  heated  and  cooled  plate  the  differential 
thermal  strains  within  the  plate  in  part  can  be  accommodated  by  bending. 

This  is  not  the  case  for  symmetric  heating  and  cooling. 

Figure  5  shows  the  spatial  distribution  of  the  stresses  at  t  -*•  °°  for 

h  =  °°  (or  finite)  for  a  number  of  values  of  pa.  It  is  of  interest  to  note 

that  the  maximum  tensile  stresses  are  formed  in  the  front  surface  of  the 
plate  which  is  at  the  highest  temperature.  These  results  were  unexpected 
and  initially  thought  to  be  incorrect.  However,  a  separate  study  devoted 

g 

to  this  effect  to  be  published  elsewhere  showed  that  these  results 
indeed  were  correct  and  could  be  attributed  to  the  nature  of  the  concavity 
of  the  temperature  distribution.  In  principle,  the  present  problem  involves 
an  assymmetrically  non-uniform  internal  heat  generation,  which  results  in 
a  temperature  distribution  which  is  concave  downward.  This  temperature 
distribution  contrasts  with  the  usual  concave  upward  temperature  distri¬ 
butions  found  in  many  cases  of  thermal  stresses  with  the  heat  exchanged 
only  at  the  outside  boundary  of  the  structure. 
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A  similar  effect  is  noted  for  the  stress  distributions  for  h  =  0 
shown  in  Fig.  6,  which  indicated  that  at  least  for  the  small  values  of 
ya,  the  maximum  tensile  stress  also  is  found  at  the  position  of  maximum 
temperature  at  x  =  -a,  i.e.,  the  front  surface  of  the  plate. 

The  material  properties  which  control  the  thermal  stress  resistance 
for  the  assymmetrically  heated  and  cooled  flat  plate  can  be  determined 
most  conveniently  from  the  equations  for  q^  for  h  =  0,“>  and  the  two 
limiting  cases  of  the  optically  thin  and  thick  plates.  From  these  equa¬ 
tions  the  following  thermal  stress  resistance  parameters  can  be  obtained: 

f S  (1  -  v)k/ctEye}  (h  =  0,°°;  pa<<l)  (22) 

and 

{ S ^  ( 1  -  v)k/aEe}  (h  =  0;  ya  -*•  °°)  (23) 

These  parameters  indicate  that  high  thermal  stress  resistance  for 
heating  and  cooling  requires  high  values  of  tensile  strength  and  thermal 
conductivity  in  combination  with  low  values  of  the  coefficient  of  thermal 
expansion.  Young's  modulus,  absorption  coefficient  and  emissivity  (high 
values  of  reflectivity).  It  should  be  noted  that  for  h  =  0  and  ya<<l, 
thermal  stress  resistance  is  inversely  proportional  to  the  first  power 
of  the  absorption  coefficient.  This  contrasts  with  the  findings  for  the 
symmetrically  heated  plate  with  h  =  0  and  ua<<l,  which  showed  that  the 
thermal  stress  resistance  was  an  inverse  function  of  the  cube  of  the 
absorption  coefficient.  This  suggests  that  in  general  the  relative  role 
of  the  material  properties  which  effect  the  thermal  stress  resistance  of 
brittle  ceramics  is  a  function  of  the  specific  heating  and  cooling  con¬ 
ditions  and  will  require  analysis  for  each  individual  case. 
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It  should  be  noted  that  the  material  requirements  for  high  thermal 


stress  resistance  may  be  incompatible  with  the  material  requirements  for 
other  performance  criteria.  For  high  efficiency  of  flat  plate  solar 
collectors,  it  will  be  required  to  have  high  values  of  the  absorption 
coefficient  and  emissivity  to  keep  the  amount  of  transmitted  and  reflected 
radiation  to  a  minimum.  Clearly,  for  the  selection  of  materials  for 
solar  energy  collectors  trade-offs  between  maximum  thermal  stress  resis¬ 
tance  and  collector  efficiency  will  have  to  be  made.  Similar  trade-offs 
of  perhaps  another  nature  may  have  to  be  made  for  aerospace  structures 
of  components  subjected  to  intense  radiation  from  solar  or  other  sources. 

Focussing  attention  on  the  absorption  coefficient,  Fig.  4  suggests 
that  for  a  non-zero  heat  transfer  coefficient  and  steady-state  operation, 
intermediate  values  of  the  optical  thickness,  especially  near  ya  =  2, 
should  be  avoided.  The  incidence  of  thermal  stress  fracture  can  be 
minimized  by  either  setting  ya  -*•  0  or  ua  -*■  °°.  Zero  stresses  are  obtained 
for  the  latter  case  because  at  steady-state  (t  ■+  00 )  the  temperature  dis¬ 
tribution  within  the  plate  becomes  linear  and  therefore  stress-free. 

High  values  of  the  transient  stresses  can  be  minimized  by  gradually 
increasing  the  heat  flux  towards  the  final  steady-state  value  rather 
than  subjecting  the  plate  to  the  full  value  of  heat  flux  at  t  =  0. 

The  observation  that  the  maximum  tensile  stresses  occur  at  the  hot¬ 
test  part  of  the  plate  (i.e.,  the  front  surface)  could  be  critical  for 
environmental  conditions  which  promote  fatigue  by  stress-corrosion  or 
other  mechanisms.  Stress-corrosion  generally  is  a  thermally  activated 
process  and  is  expected  to  occur  at  a  maximum  rate  for  the  temperature 
and  stress  distributions  prevalent  for  the  heating  and  cooling  conditions 
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of  the  present  study.  The  magnitude  of  temperature  at  which  the  plate 
operates,  however,  can  be  minimized  by  having  a  high  value  of  the  heat 
transfer  coefficient.  As  noted  earlier,  the  steady-state  stress  value 
is  independent  of  the  heat  transfer  coefficient,  but  the  steady-state 
temperature  level  is  inversely  proportional  to  the  heat  transfer  coef¬ 
ficient.  It  appears  then  that,  in  general,  the  incidence  of  thermal 
stress  fracture  in  brittle  ceramic  flat  plates  subjected  to  radiation 
heating  can  be  kept  to  a  minimum  by  a  judicious  selection  of  materials, 
in  combination  with  the  optimum  choice  of  optical  thickness  and  cooling 
conditions . 
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ABSTRACT 


The  degree  of  micro-cracking  in  BeO-SiC  composites  under  the  influence 
of  the  internal  stresses  which  arise  due  to  mismatch  in  the  coefficients  of 
thermal  expansion  was  studied  by  measuring  thermal  diffusivity  of  the  com¬ 
posite  as  a  function  of  SiC  content.  The  experimental  results  on  diffusivity 
measurement  indicated  that  the  micro-cracking  was  extensive  at  approx.  30  and 
80%  SiC  content  and  minimum  at  ~  50%  SiC  content.  This  was  in  agreement 
with  similar  observation  of  the  effect  of  micro-cracking  on  strength,  and 
fracture  energy  of  other  investigators. 

A  stress  analysis  indicated  that  these  observations  cannot  be  attributed 
to  the  difference  in  magnitude  of  internal  stresses,  but  instead  are  thought 
to  be  due  to  statistical  effects  controlled  by  difference  in  stress-state 
in  accordance  with  suggestions  by  Evans. 
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1.  Introduction 

The  development  of  composites  has  led  to  materials  with  advantageous 
properties  not  found  in  single-phase  materials.  Mismatches  in  the  coefficient 
of  thermal  expansion  of  the  individual  component  within  the  composite  can  lead 
to  the  development  of  high  levels  of  internal  stresses  on  cooling  the  compo¬ 
site  from  its  manufacturing  temperature1-7.  In  composites  in  which  one  of 
the  components  is  a  polymer  or  a  metal,  relaxation  of  such  internal  stresses 
can  occur  by  viscous  or  plastic  flow  even  at  relatively  low  temperatures.  In 
brittle  composites,  internal  stress  relaxation  can  occur  by  diffusional  creep, 
but  only  at  high  temperatures  near  the  manufacturing  temperature.  Below  these 
temperatures  no  stress  relaxation  can  occur  so  that  in  brittle  composites,  for 
a  given  degree  of  mismatch  of  the  coefficients  of  thermal  expansion,  the  magni¬ 
tude  of  the  internal  stresses  will  be  particularly  high.  Such  stresses,  if 
sufficiently  large,  can  lead  to  extensive  "micro-cracking. "  Such  micro-cracking 
can  occur  at  the  boundaries  between  the  components  or  within  the  component  with 
the  higher  coefficient  of  thermal  expansion  which  generally  is  subjected  to  the 
high  value  of  tensile  stress.  It  should  be  noted  that  extensive  micro-cracking 
can  also  occur  in  single  phase  non-cubic  polycrystalline  materials  which  ex- 

8,9 

hibit  a  large  degree  of  anisotropy  in  thermal  expansion  behavior 

Micro-cracking  in  brittle  composites  or  polycrystall ine  materials  is  not 
necessarily  disadvantageous.  In  fact,  it  can  lead  to  major  improvements  in 

many  engineering  properties.  As  observed  experimentally  micro-cracking  can 

10 

lead  to  siginficant  increase  in  fracture  toughness  ,  presumably  due  to  the 

formation  of  a  process-zone  ahead  of  a  propagating  crack.  Micro-cracking  can 

10,11,12 

also  lead  to  major  decrease  in  elastic  behavior  and  the  thermal  con¬ 

ductivity  and  diffusivity1 3’lu.  Due  to  a  large  increase  in  the  strain-at- 


fracture  in  combination  with  a  stable  mode  of  crack  propagation,  micro-cracked 
solids  appear  to  be  excellent  candidate  materials  for  applications  at  high 
temperature  involving  severe  thermal  shock10’10.  In  particular,  micro-cracked 
materials  probably  are  the  only  class  of  engineering  materials  appropriate  for 
conditions  which  require  high  thermal  shock  resistance  in  combination  with 
good  thermal  insulating  ability10.  Finally,  micro-cracking  can  enhance  the 
machinabil ity  of  brittle  materials  substantially  to  an  extent  that  they  can 
be  shaped  by  conventional  metal  working  techniques17. 

As  indicated  by  a  large  number  of  literature  studies1-9,  methods  for  the 
calculation  of  the  magnitude  and  distribution  of  internal  stresses  which  arise 
from  non-uniformities  of  the  coefficients  of  thermal  expansion  in  composites 
or  poly-crystalline  materials  appear  to  be  well  developed.  In  general, 
these  calculations  show  that  the  magnitude  of  these  internal  stresses  is  a 
function  of  the  mismatch  or  degree  of  anisotropy  of  the  thermal  expansion  co¬ 
efficients,  the  range  of  temperature  over  which  the  composite  or  polycrystal¬ 
line  material  is  cooled  without  any  stress  relaxation,  as  well  as  the  phase 
distribution  of  the  individual  components  or  the  relative  orientation  of  the 
individual  grains  within  the  poly-crystal.  Of  significance  to  the  present 
study  is  that  these  studies  show  that  the  magnitude  of  the  internal  stresses 
for  a  given  composite  or  non-cubic  polycrystalline  material  is  not  a  function 
of  the  size  of  the  components  or  the  grains.  For  this  reason,  if  micro¬ 
cracking  were  controlled  solely  by  the  magnitude  of  internal  stress,  the 
degree  of  micro-cracking  in  terms  of  the  number  of  micro-cracks  formed  should 
be  independent  of  the  scale  of  the  microstructure. 

This  latter  conclusion,  however,  is  in  contradiction  with  experimental 
observations.  As  observed  by  Davidge  and  Green1R,  for  a  continuous  matrix 
dispersed  phase  system,  a  minimum  dispersion  size  was  necessary  for  micro- 
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cracking  to  occur.  Also,  in  polycrystall ine  materials  with  high  thermal 
expansion  anisotropy,  a  critical  grain  size  is  required  before  micro-crack¬ 
ing  can  take  place.  It  appears  then  that  the  presence  of  internal  stresses 
is  a  necessary  but  not  sufficient  condition  for  micro-cracking  and  that  other 
factors  need  to  be  considered. 

In  providing  an  explanation  for  their  observations,  Davidge  and  Green18 
noted  that  the  elastic  energy  available  for  crack  propagation  was  proportional 
to  the  third  power  of  the  inclusion  size.  However,  the  total  energy  required 
to  propagate  a  radial  crack  is  proportional  to  the  crack  size.  For  this 
reason,  below  a  given  inclusion  size  insufficient  elastic  energy  was  avail¬ 
able  for  the  formation  of  a  micro-crack  of  sufficient  size  to  give  sufficient 
stress  relief.  A  similar  explanation  was  offered  by  Kuszyk  and  Bradt10  to 
explain  the  grain  size  effect  on  micro-cracking  in  polycrystalline  single  phase 
materials.  It  should  be  noted,  however,  that  this  explanation  considers  only 
the  energy  conditions  after  microcrack-formation,  but  does  not  address  the 
condition  of  the  onset  of  formation  of  micro-cracks. 

Lange19  based  his  analysis  of  micro-crack  formation  on  an  energy  balance 
approach  and  found  that  below  a  given  value  of  the  product  of  the  square  of 
the  stress  and  the  particle  radius,  micro-crack  formation  cannot  occur  of 
the  size  of  the  pre-existing  flaw  from  which  the  micro-crack  originated. 

More  recently,  Evans"1  pointed  out  that  micro-crack  formation  most  likely 
occurs  from  precursor  flaws  in  the  form  of  pores  at  grain-boundaries  or 
triple  points.  The  size  of  such  pore  frequently  is  directly  related  to  the 
grain  size.  For  this  reason,  the  larger  the  grain  for  a  given  level  of  in¬ 
ternal  stress,  the  larger  will  be  the  stress  intensity  factor  (Kj)  for  the 
precursor  flaw.  Micro-cracking  will  occur  whenever  Kj  3  K j c ,  the  critical 
stress  intensity  factor.  This  latter  mechanism  provides  a  ready  explanation 
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for  the  observed  component  or  grain  size  effect  on  the  incidence  of  micro- 
cracking  in  brittle  materials.  It  is  expected  that  micro-cracking  originating 
from  grain-boundary  pores  would  be  governed  by  statistical  variables  such  as 
the  size  distribution  and  orientation  of  the  pores  as  well  as  the  magnitude 
and  distribution  of  the  internal  stresses.  Experimental  data  in  support  of 
this  hypothesis  are  presented  in  the  present  paper. 

2.  Experimental 

2.1  Materials  and  Experimental  Method 

Nearly  fully-dense  samples  of  composites  of  BeO  and  SiC  made  by  hot- 

★ 

pressing  were  obtained  from  a  commercial  source.  Specimens  of  fully  dense 
SiC  made  by  conventional  sintering  were  obtained  as  well.  Figure  1  shows 
scanning  electron  micrographs  of  composites  containing  30,  SO  and  80  wt.% 

SiC.  The  grain  size  of  BeO  and  SiC  was  below  5  pm.  Due  to  the  relatively 
large  mismatch  in  their  coefficients  of  thermal  expansion,  BeO-SiC  composites 
exhibit  extensive  micro-cracking.  As  shown  in  Fig.  1,  the  isolated  inter- 
grannular  cracking  is  frequently  observed  to  initiate  from  a  triple  point  and 
propagate  along  one  or  two  adjacent  grain  faces  in  agreement  with  Evans’1. 

As  shown  analytically70,  micro-cracking  has  a  pronounced  effect  on  ther¬ 
mal  conductivity.  For  this  reason,  the  degree  of  micro-cracking  was  monitored 
by  measurements  of  the  thermal  diffusivity  by  the  laser-flash  method21,  using 
equipment  described  elsewhere13.  Measurements  were  made  over  the  range  of 
temperature  from  about  300  to  1400°C.  One  particular  advantage  of  using  the 
measurement  of  thermal  diffusivity  for  monitoring  the  extent  of  micro-cracking 
is  that  the  results  are  not  affected  by  additional  micro-crack  formation 
due  to  mechanical  loading  as  is  expected  to  be  the  case  during  measurements 
of  strength  and  fracture  toughness. 

*Ceradyne,  Inc.,  Santa  Ara,  CA,  USA;  Code  Ceralloy  2700  Series. 
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2.2  Experimental  Results 

Figure  2  shows  the  experimental  data  for  the  thermal  diffusivity  as  a 
function  of  temperature  for  BeO,  SiC  and  two  of  the  BeO-SiC  composites. 

The  data  shown  in  Fig.  2  were  obtained  by  heating  the  specimens  slowly  from 
room  temperature  to  1400°C.  On  cooling  a  hysteresis  effect  was  noted.  A 
detailed  nalysis  of  this  latter  observation  is  deferred  to  a  later  publica¬ 
tion. 

Figure  3  shows  the  thermal  diffusivity  as  a  function  of  SiC  content 
for  a  number  of  temperatures  as  obtained  by  smooth  curves  drawn  through  the 
experimental  data. 

Figure  4  shows  the  value  of  the  therma1  conductivity  (K)  at  400°C  cal¬ 
culated  from  the  corresponding  curve  for  the  thermal  diffusivity  in  Fig.  3 
using  the  relation: 

K  =  t.  c  (1) 

/ 

where  a  is  the  thermal  diffusivity,  p  is  the  density  and  c  is  the  specific 
heat.  Latter  two  values  as  a  function  of  SiC  content  were  assumed  to  vary 
linearly  between  the  end  values  for  SiC  and  BeO. 

Included  in  Fig.  4  is  the  value  for  the  thermal  conductivity  as  a 
function  of  SiC  content  for  a  non-microcracked  composite  calculated  from  the 
Rayleigh-Maxwel  1  relation'’'’’71 

1+2V(1 -K  /K  )/(2K  /K  +1 ) 

Kc  =  Km  1  -  V(l"-K/Kp?7TKm/Kp+?y  (2) 

where  K  is  the  thermal  conductivity,  the  subscripts  c,  m  and  p  refer  to  the 
composite,  matrix  and  dispersed  phase,  resp.  For  the  numerical  calculation  of 
conductivity  by  eq.  2  the  following  data  were  used:  SiC,  p  =  3.21  gm/cc, 
c  =  0.255  cal  gm"^nK  ^  and  BeO,  p  =  3.008  gm/cc,  c  =  0.408  cal  gm  ^ °K  ^ . 
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3.  Discussion  and  Conclusions 

The  discussion  will  primarily  concentrate  on  the  mi crostructural  and 
mechanical  variables  which  affect  the  formation  of  micro-cracks.  A  compara¬ 
tive  study  of  the  effect  of  micro-cracking  on  heat  transport  properties  for 
a  number  of  different  micro-cracked  composites  is  deferred  to  a  subsequent 
publ ication. 

The  thermal  conductivity  for  the  BeO-SiC  composites  calculated  from 
the  Rayleigh-Maxwell  theory  varies  monotonically  between  the  values  of  the 
thermal  conductivity  of  the  BeO  and  SiC.  In  fact  regardless  of  the  phase- 
distribution  of  each  component,  the  composite  thermal  conductivity  should 
fall  between  the  corresponding  values  of  the  two  components  For  this 
reason  the  value  for  K  for  the  SiC  represents  the  lower  bound  on  the  thermal 
conductivity  of  the  composite  regardless  of  volume  fraction.  However,  the 
variation  of  the  thermal  conductivity  as  a  function  of  silicon  carbide  content 
inferred  from  the  experimental  data  for  the  thermal  diffusivity  lies  below 
the  value  for  the  SiC.  Since  these  composites  are  fully-dense  (or  nearly  so), 
this  effect  cannot  be  attributed  to  the  presence  of  a  pore  phase.  Instead, 
the  decrease  of  the  thermal  conductivity  must  tie  attributed  to  the  presence 
of  the  micro-cracks. 

It  is  of  particular  interest  to  note  that  the  decrease  in  thermal  conduc¬ 
tivity  is  most  pronounced  at  approximately  30  and  BO  SiC.  In  contrast, 
at  approximately  50'  SiC  the  decrease  in  thermal  conductivity  is  rather  small, 
if  existent  at  all.  This  observation  suggests  that  micro-cracking  is  most 
pronounced  at  30  and  803  content  of  SiC  and  appears  to  be  absent  at  approxi¬ 
mately  50%  SiC.  Support  for  this  latter  conclusion  is  provided  by  experimental 
data?s  made  available  to  the  present  authors  prior  to  future  publication  in  detail 
at  a  later  date,  for  the  strength  and  fracture  energy  for  these  composites. 
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At  approximately  20  and  75%  SiC,  these  composites  have  higher  values  for 
tensile  strength  (measured  in  bending)  and  fracture  energy  as  compared  to  the 
values  for  single  phase  SiC  and  BeO.  At  40%  SiC,  however,  the  values  of 
strength  are  comparable.  In  other  words,  at  these  compositions  which  exhi¬ 
bit  a  decrease  in  thermal  conductivity,  a  corresponding  increase  in  strength 
and  fracture  energy  is  found.  The  minima  in  the  strength  and  fracture  energy 
plots,  however,  occur  at  a  slightly  lower  SiC  content  than  the  corresponding 
SiC  content  for  the  maxima  in  the  diffusivity  plot.  This  is  believed  to  be  due  to 
the  fact  that  during  the  strength  and  fracture  energy  measurements  additional 
micro-cracks  could  have  been  formed  during  the  loading  period.  These  opposite 
effects  are  entirely  consistent  with  the  presence  of  micro-cracks  which  are 
known  to  enhance  fracture  toughness10  and  to  decrease  the  thermal  conductivity*- 1 . 
Further  supporting  evidence  for  the  low  degree  of  mi cro-cracking  can  be  obtained 
from  the  experimental  data  for  the  coefficient  of  thermal  expansion  of  BeO-SiC 
composites  measured  by  Rossi76.  At  50  vol.%  SiC  the  coefficient  of  thermal  ex¬ 
pansion  lies  relatively  much  closer  to  the  coefficient  of  thermal  expansion  cal- 

★ 

culated  from  composite  theory  than  at  25  and  75'  SiC. 

As  shown  analytically,  the  effect  of  randomly  oriented  circular  cracks 
of  equal  size  on  thermal  conductivity  is  given  by16: 

K  =  K  [1  +  8  Nb3/9]_1  (3) 

where  N  is  the  number  of  cracks  per  unit  volume,  and  b  is  the  crack  radius. 
Unfortunately,  the  density,  size  and  geometry  of  the  micro-cracks  cannot  be 

*It  should  be  noted  that  Rossi  draws  a  smooth  curve  through  his  data.  Close 
examination  of  the  experimental  data,  however,  suggests  that  a  better  curve 
would  have  the  same  general  form  as  the  curves  in  Figure  3. 
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determined  with  any  accuracy  by  metal lographi c  methods.  For  this  reason, 
the  decrease  in  thermal  conductivity  cannot  be  estimated  from  prior  in- 
formation.  Instead,  from  eq.  3,  the  value  of  Nb  can  be  assessed  from  the 
experimental  data  for  thermal  diffusivity.  These  results  are  shown  in 
Fig.  5  for  the  case  of  uniformly  distributed  micro-cracks  within  the  com¬ 
posite  as  a  whole.  Since  the  micro-cracks  are  more  likely  to  be  found  within 
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the  BeO  phase,  Nb  was  also  evaluated  by  means  of  eq.  2  in  which  the  thermal 
conductivity  of  the  BeO  phase  was  assumed  to  be  given  by  eq.  3.  As  indicated 
in  Fig.  5,  the  value  of  Nb  for  micro-cracks  in  the  BeO  phase  can  go  to  high 
values,  as  expected,  since  at  the  higher  fractions  or  SiC  the  BeO  will  exhibit 
the  highest  degree  of  micro-cracking.  A  value  of  Nb  ~  0.5  for  uniformly  dis¬ 
tributed  micro-cracking  at  -  80  SiC  content  suggests  that  for  a  micro-crack  rad 
equal  to  the  grain  size  5um)  approximately  one  half  of  the  BeO  grains 
contain  a  micro-crack.  It  should  be  noted  that  an  estimate  of  the  thermal 
conductivity  of  a  micro-cracked  composite  obtained  by  substitution  of  eq.  3 
in  eq.  2  for  the  conductivity  of  the  micro-cracked  phase  assumes  that  the 
effect  of  the  micro-cracking  and  the  presence  of  the  dispersions  on  the 
thermal  conductivity,  in  fact,  are  uncoupled.  This  latter  assumption,  however 
is  not  expected  to  be  valid  since  the  micro-cracking  in  the  BeO  phase  is  ex¬ 
pected  to  occur  in  the  immediate  vicinity  of  the  SiC  matrix.  An  analysis  of 
the  coupling  effect  of  dispersions  and  micro-cracking  on  thermal  conductivity 
is  recommended  as  a  fruitful  area  for  further  theoretical  and  experimental 
work . 

The  question  as  to  why  at  the  intermediate  volume  fractions  of  the  SiC 
phase  micro-cracking  appears  to  be  largely  absent,  still  remains  to  be 
answered.  To  find  a  solution  to  this  problem,  the  internal  stresses  at  the 
BeO-SiC  interface  were  calculated  for  a  composite  sphere  model  in  which  either 


the  BeO  or  the  Si C  represents  the  outer  layer.  With  SiC  on  the  inside 
represents  the  SiC-BeO  composite  with  low  volume  fractions  of  SiC.  With 
the  BeO  on  the  inside  of  the  composite,  represents  the  microstructure  at 
high  volume  fractions  SiC  at  which  the  BeO  takes  the  form  of  the  dispersed 
phase.  At  infinite  dilute  volume  fractions  the  answers  obtained  agreed  with 
those  obtained  from  solutions  presented  by  Seising1.  Figure  6  shows  the 
tensile  radial  and  tangential  stresses  as  a  function  of  volume  fraction  SiC 
for  the  BeO  or  SiC  as  the  outside  layer  of  the  composite  sphere.  It  may 
be  noted  that  the  stresses  increase  linearly  with  increasing  volume  fraction 
SiC.  Clearly,  no  reason  exists  for  the  absence  of  micro-cracking  at  the 
intermediate  volume  fractions  of  SiC. 

A  spherical  composite  model  considered  above,  however,  assumes  that 
either  the  SiC  or  BeO  phase  is  the  dispersed  phase  over  the  total  volume 
fraction  of  SiC.  Near  50"  SiC,  however,  BeO  and  SiC  are  both  expected  to 
be  continuous.  For  this  reason,  a  composite  cylinder  model  for  the  assess¬ 
ment  of  the  internal  stresses  appears  more  appropriate. 

For  such  a  composite  cylinder  model,  solutions  for  the  internal  stresses 
due  to  the  mismatch  of  coefficients  of  thermal  expansion  were  derived  by 
Gatewood?/.  For  the  SiC-BeO  system  the  magnitude  of  these  tensile  internal 
stresses  (in  the  BeO  phase)  are  shown  in  Fig.  7.  Again,  as  for  the  composite 
sphere  model,  these  stresses  increase  with  increasing  volume  fraction  SiC. 

At  first  sight,  from  the  point  of  view  of  the  magnitude  of  stress,  no  reason 
appears  to  exist  as  to  why  at  approx.  50  vol.'  SiC  the  micro-cracking  seems 
suppressed.  However,  it  is  well  known  that  fracture  of  brittle  materials 
is  governed  by  statistical  variables  affected  by  flaw  size  distribution, 
and  orientation  as  well  as  by  the  stress  state.  In  respect  to  this  effect  of 
stress  state  on  failure  probability,  it  should  be  noted  that  in  the  composite 
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sphere  model  the  tensile  stress  in  the  BeO  phase  at  the  BeO-SiC  interface 
is  uniform  bi-axial  in  nature.  Therefore  at  a  radially  oriented  pore  at  a 
triple  point  at  the  BeO-SiC  interface,  micro-crack  formation  can  occur  re¬ 
gardless  of  the  pore's  other  orientation.  For  the  composite  cylinder,  however, 
which  more  closely  represents  the  microstructure  of  the  BeO-SiC  composite 
near  50  %  SiC,  the  maximum  tensile  stress  in  the  BeO  appears  in  the  axial 
direction.  This  axial  stress  is  of  a  uniaxial  nature.  For  this  reason,  at 
a  flat  pore  oriented  radially  at  the  BeO-SiC  interface,  micro-crack  formation 
can  occur  only  if  the  plane  of  the  pore  is  oriented  perpendicul arly  to  the 
axial  direction.  For  this  reason,  the  probability  of  micro-crack  formation 
in  the  composite  cylinder  model  is  less  than  the  corresponding  probability 
in  the  composite  sphere  model.  Since  the  composite  cylinder  model  approached 
the  phase  distribution  of  the  present  composites  near  50.  SiC,  the  observations 
reported  herein  are  consistent  with  expected  behavior  governed  by  statistical 
variables . 

In  general,  then,  the  present  results  and  discussion  suggest  that  theories 
for  the  formation  of  micro-cracks  based  on  the  magnitude  of  stress,  elastic 
energy  principles  or  stress  intensity  approaches  also  should  include  the 
statistical  aspects  of  brittle  fracture  such  as  the  size  distribution  of  the 
precursor  flaws,  their  orientation  as  well  as  the  stress-state  of  the  residual 
stress  fields. 
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Fig.  1.  Scanning  electron  fractographs  of  Beo-SiC 
composites  showing  the  triple  points  and 
micro-cracks  initiating  from  triple  point. 
Composition:  (a)  30*  SiC  and  (b)  50  SiC  by 
weight. 


Scanning  electron  fractograph  of  BeO-SiC 
composite  (containing  80%  SiC  by  weight)  show¬ 
ing  the  triple  points  and  micro-cracks  ini¬ 
tiating  from  triple  point. 
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The  values  of  thermal  conductivity  of  BeO-SiC  composites  as  a  function  of  SiC  content  cal 
culated  from  the  diffusivity  values  in  Figure  3  at  400°C.  Included  in  the  figure  is  also 
the  calculated  values  for  thermal  conductivity  based  on  Rayleigh  Maxwell  theory. 
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ABSTRACT 


The  transient  thermal  stresses  in  a  long  cylinder  with  square 
cross-section  resulting  from  convective  heating  or  cooling  by  a  sudden 
change  in  ambient  temperature  were  calculated  by  the  finite  element 
method.  On  heating,  the  maximum  tensile  stresses  occur  in  the  center 
of  the  rod,  whereas  the  maximum  principle  stresses  occur  along  the 
diagonal  of  the  cross-section.  On  cooling,  the  maximum  tensile 
stresses  occur  at  the  center  of  the  faces.  The  results  obtained  are 
discussed  in  terms  of  experimental  measurement  of  the  thermal  stress 
resistance  of  brittle  ceramics.  Possible  failure  modes  during  heating 
and  cooling  are  discussed.  Comparison  with  literature  data  shows  that 
for  a  given  value  of  the  Biot  number  the  magnitude  of  the  maximum  ten¬ 
sile  stresses  on  cooling,  in  a  rod  with  square  cross-section  exceeds 
those  in  a  rod  with  circular  cross-section. 


I.  INTRODUCTION 

Brittle  ceramic  materials  are  highly  susceptible  to  catastrophic  failure 
under  non-isothermal  conditions  involving  high  levels  of  thermal  stress.^- 
For  purposes  of  reliable  engineering  design,  characterization  of  the  thermal 
stress  resistance  of  ceramics  for  high-temperature  structural  application  is 
essential . 

A  favorite  method  for  measuring  the  thermal  stress  resistance  of  ceramics 

consists  of  quenching  appropriate  specimens  into  a  quenching  medium  which  is 

at  a  temperature  different  from  the  initial  temperature  of  the  specimens. 

2  3  4  5 

Such  quenching  media  may  consist  of  water,  ’  liquid  metals,  oils,  fluidized 

beds*3  or  eutectic  mixtures  of  salts.'7  The  initial  temperature  of  the  specimens 

may  be  above  or  below  the  temperature  of  the  quenching  medium.  The  specimen 

4  7  2 

geometry  can  consist  of  spheres,  ’  solid  cylinders  with  circular  or  square 
3 

cross-section. 

g 

For  interpretation  of  the  results  for  the  spherical  specimen  or  the  solid 
9 

circular  cylinders,  analytical  solutions  for  the  transient  temperatures  and 
thermal  stresses  are  available  in  the  literature.  Unfortunately  for  cylindri¬ 
cal  specimens  with  square  cross-section  no  such  solutions  are  available.  For 
this  reason,  experimental  quenching  data  for  such  specimen  can  be  used  to 
establish  the  relative  thermal  stress  resistance  of  different  materials. 

But  a  quantitative  estimate  of  the  magnitude  of  thermal  stress  at  which  a  square 
specimen  fail  cannot  be  made.  A  brief  analytical  study  carried  bv  one 
of  the  present  authors  some  years  ago  revealed  that  the  lack  of  availability 
of  analytical  solutions  for  the  thermal  stresses  in  square  cylinders,  under 
convective  heat  transfer,  appears  to  be  attributable  to  the  lack  of  rapid  con¬ 
vergence  of  the  product  of  two  infinite  series  involved  in  the  general  solu- 
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tion.  It  may  be  noted  that  under  conditions  of  a  constant  heat  flux,  no 

,  ,  ,  ,10 

such  problems  are  encountered. 

An  alternative  approach  to  evaluation  of  thermal  stresses  in  square  cy¬ 
linders  is  to  solve  them  by  numerical  methods,  such  as  the  finite  element 
method,  which  formed  the  basis  of  the  present  study. 

II.  PROCEDURE 

The  cylinder  with  square  cross-section  located  at  -a  <  x  <  a  and 
-a  <  y  <  a  was  assumed  to  be  infinite  in  extent  in  the  Z-direction  in  order 
to  avoid  end-effects.  The  cylinder  initially  at  thermal  equilibrium  is 
subjected  to  thermal  stresses  by  heating  or  cooling  by  an  instantaneous 
change  in  ambient  temperature.  Heat  transfer  was  assumed  to  occur  by  Newtonian 
convection,  i.e.,  a  rate  of  heat  transfer  proportional  to  the  instantaneous 
temperature  difference  between  the  surface  and  the  ambient  environment.  Heat 
transfer  was  assumed  to  be  uniform  along  the  length  in  order  to  avoid  an 
axial  component  of  the  heat  flow.  The  severity  of  the  heat  flux  was  ex¬ 
pressed  in  terms  of  the  non-dimensional  Biot  number,  B  =  ah/k  where  a  is  the 
half-thickness  of  the  cylinder,  h  is  the  heat  transfer  coefficient  and  k  is 
the  thermal  conductivity  of  the  specimen. 

The  transient  temperatures  and  thermal  stresses  were  evaluated  with 
the  finite  element  computer  program  developed  for  a  number  of  previous  studies. 
The  reliability  of  this  program  was  checked  out  by  verifying  analytical  re¬ 
sults  reported  in  the  literature  for  the  thermal  stresses  in  a  solid  circular 
9 

cylinder  subjected  to  convective  heat  transfer  and  for  the  thermal  stresses 

in  a  hollow  cylinder  with  spatially  varying  thermal  conductivity  under  con- 

19 

ditions  of  steady-state  heat  transfer.  Excellent  agreement  between  the 
analytical  and  numericai  results  was  obtained. 


In  view  of  the  symmetry  of  the  temperature  and  thermal  stress  distri¬ 
bution  around  the  x-  and  y-  axes  of  the  square  cross-section  of  the  cylinder, 
the  finite  element  computations  were  simplified  by  considering  only  one  of 
the  four  quadrants  of  the  cross-section  located  at  0  <  x  <  a  and  0  <  y  <  a. 
This  quadrant  was  subdivided  into  64  elements  of  equal  size  with  the  element 
at  the  comer  (a, a)  further  subdivided  into  4x4  equal  size  element.  For 
values  of  the  Biot  number  as  high  as  20,  this  finite  element  mesh  gave  identi¬ 
cal  results  with  a  mesh  for  which  the  corner  element  was  divided  into  6x6 
sub-elements.  The  calculations  were  not  extended  to  values  of  the  Biot  number 
8  >  20,  as  the  high  temperature  gradients  for  these  high  values  of  8,  required 
an  even  finer  mesh  size,  requiring  excessive  computer  time.  All  pertinent 
material  properties  were  assumed  to  be  independent  of  temperature. 

For  convenience  all  the  results  for  the  transient  temperature,  T  are 

*  * 

reported  in  terms  of  the  non-dimensional  temperature,  T  defined  by  T  = 

(T  -  T)/T  -  Tj.)  for  convective  cooling  and  T  =  (T  -  T  ^ )  /  ( T  ^  -  T^)  for 

convective  heating,  where  T^  is  the  initial  temperature  of  the  rod,  and  T^  is 

the  temperature  of  the  ambient  to  which  the*  rod  is  subjected. 

The  transient  thermal  stresses,  are  reported  in  terms  of  the  non- 
* 

dimensional  stress,  o  =  u(l- .>)  h  KAT,  where  v  is  Poisson  s  ratio,  a  is  the 

coefficient  of  thermal  expansion,  E  is  Young's  modulus  of  elasticity  and 

AT  =  [ T ^  -  T.|.  The  time  dependence  of  the  temperatures  and  stresses  is 

*  2 

expressed  in  terms  of  the  non-dimensional  t ime  l  =  ut/a“  where  a  is  the 
thermal  diffusivity  and  t  is  the  time.  For  simplicity  the  term  "non-dimen¬ 
sional"  will  be  deleted  from  the  further  discussion. 

III.  RESULTS  AND  DISCUSSION 

Typical  temperature  and  stress  distributions  will  he  examined  first. 
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Because  of  the  symmetry  of  the  square  cross-section  about  both  the  x-  and 
y-  axis  as  well  as  the  diagonal,  it  is  sufficient  to  present  the  temperature 
and  stresses  within  an  octant  only. 

Figure  la  shows  the  distribution  of  the  temperature  and  the  transverse 
stress,  o  at  the  time  it  reaches  its  maximum  value  for  convective  heating 

y 

with  6=5.  The  transverse  stress  at  (o,  o)  reaches  its  maximum  value  at 
the  same  instant  as  the  axial  stress,  (o,  o) .  The  maximum  value  of  the 
compressive  transverse  stress  (o  )  is  obtained  at  the  center  of  the  face 
(a,  o) .  As  expected,  the  edge  of  the  rod  has  attained  the  highest  value  of 
temperature  at  the  time  of  maximum  transverse  stress. 

Figure  lb  shows  the  axial  and  principal  stresses  at  the  time  of  their 
respective  maximum  values  on  heating  for  Biot  number,  S  =  5.  The  maximum 
value  of  the  axial  tensile  stress  occurs  at  the  center  of  the  rod.  Apart 
from  a  slight  stress  gradient  near  the  edge,  the  axial  stresses  are  fairly 
uniform  along  the  faces  of  the  rod.  The  maximum  values  of  the  principal 
stress  occur  along  the  diagnoal  of  the  cross-section,  in  general  agreement 
with  the  results  of  Cheung  et.  al.  ^  for  heating  of  a  rectangular  prism  by 
a  constant  heat  flux. 

Figure  lc  shows  the  temperatures  and  axial  stresses  at  the  instant 
of  maximum  axial  stress,  (a,  a)  for  convective  cooling  with  P  =  5.  The 
maximum  value  of  tensile  axial  stress  occurs  at  the  edge  of  the  rod.  The 
maximum  compressive  stress  occurs  at  the  center  of  the  rod.  Again,  with  the 
exception  of  the  immediate  edge  region,  the  axial  stresses  are  nearly  uniform 
along  the  rod  face. 

Figures  2a-d  give  the  time  dependence  of  the  transverse  (a,  o) ,  the 
axial  stresses  o z  (o,  o),  n  (a,  a)  and  the  principal  stress  for  a  range  of 
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values  of  Blot  number,  (3.  Comparison  of  these  stress  values  show  that 


on  cooling,  the  transverse  stress  a (a,  o)  reaches  the  highest  value.  It 
should  also  be  noted,  however,  that  the  value  of  the  axial  stress,  a (a,  a) 
exceeds  the  value  of  transverse  stress  o  (a,  o)  at  short  times. 

y 

Figure  3a  gives  the  values  of  the  maximum  tensile  axial  (°z)>  trans¬ 
verse  (Oy) ,  and  principal  stress  (c^)  f°r  convective  heating  and  cooling 
as  a  function  of  the  Biot  number,  8.  In  general,  these  results  indicate 
that  the  tensile  stresses  developed  during  cooling  exceed  the  corresponding 
values  obtained  during  heating. 

In  Fig.  3b,  the  curves  marked  (C)  indicate  the  temperatures  at  various 
positions  in  the  rod  at  the  time  of  maximum  tensile  axial  o^(a,  o)  obtained 
during  convective  cooling  as  a  function  of  Biot  number.  Similarly,  the 
curves  marked  (H)  indicate  the  temperature  at  various  positions  at  the  time 
of  maximum  axial  stress,  a (o,  o)  at  the  center  of  the  rod.  These  results 
indicate,  as  expceted  that  the  temperature  differences  within  the  rod  at 
the  time  of  maximum  stress  increase  with  increasing  Biot  number. 

For  purposes  of  comparing  the  results  of  various  thermal  shock  tests, 
it  is  of  interest  to  compare  the  present  results  for  the  magnitude  of  the 


maximum  tensile  stresses  developed  during  cooling  in  a  cylinder  with  a 


square  cross-section  with  those  in  a  cylinder  with  circular  cross-section. 
These  results  are  shown  in  Fig.  4.  For  the  limiting  cases  of  8  =  o,  00 
these  stresses  should  be  identical.  At  intermediate  values  the  magnitude 


of  the  thermal  stresses  in  the  square  rod  exceeds  the  values  of  stress  in 
the  circular  rod.  Particularly  for  values  of  8  frequently  encountered 
in  laboratory  thermal  shock  test,  (0.5  <  8  5.0),  the  relative  difference 
in  stress  can  amount  to  as  much  as  50  percent. 
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From  the  point  of  the  interpretation  of  thermal  shock  tests,  it  is 
of  interest  to  speculate  on  the  mode  of  failure  of  the  rods  with  square 
cross-section.  Thermal  shock  by  cooling,  as  found  in  a  quench  test,  will 
be  considered  first.  If  it  is  assumed  that  failure  occurs  at  a  predeter¬ 
mined  value  of  tensile  stress,  failure  would  be  expected  to  occur  at  the 
center  of  the  face  under  the  influence  of  the  transverse  stress.  Because 
of  the  axial  stress,  the  stress  state  at  the  center  of  the  face  is  nearly 
uniform  biaxial.  This  makes  the  center  of  the  face  an  even  more  likely 
site  of  failure  because  the  biaxial  state  of  stress  increases  the  probability 
of  failure  for  those  brittle  materials  for  which  fracture  is  controlled  by 
statistical  considerations. 

In  assessing  the  site  of  failure,  however,  it  should  be  recognized 
that  during  the  initial  stages  of  the  thermal  shock  (i.e.  at  times  less  than 
the  time  of  maximum  stress),  the  axial  stress  at  the  specimen  edge  exceeds 
the  transverse  stress  at  the  face  center,  as  pointed  out  earlier.  For 
this  reason,  in  case  of  a  thermal  shock  overload  in  which  the  maximum 
stresses  exceed  the  fracture  stress  significantly,  failure  more  likely 
will  be  initiated  at  the  specimen  edge  than  at  the  face  center. 

During  preparation  of  thermal  shock  specimens  with  square  cross-section, 
such  defects  as  nicks,  scratches  and  other  flaws  more  likely  are  introduced 
in  the  specimen  edges  than  in  the  specimen  face.  This  provides  a  further 
reason,  why  thermal  stress  failure  in  square  rods  will  be  initiated  at 
the  specimen  edges  rather  than  at  the  center  of  the  faces.  This  effect  will 
be  enhanced  further  in  stress-corrosive  environments  which  will  promote  more 
sub-critical  crack  growth  in  flaws  at  the  edges  because  of  the  more  rapid 


rise  in  the  magnitude  of  thermal  stress  than  in  the  case  of  the  face  centers. 


For  these  combined  reasons,  it  is  speculated  here  that  for  most  brittle 
engineering  materials  tested  in  practice,  the  specimen  edge  may  well  be 
the  preferred  site  of  failure. 

On  heating,  failure  could  occur  under  the  influence  of  the  tensile 

stress  in  the  specimen  center.  Alternatively,  failure  can  be  initiated 

by  the  principal  stresses  along  the  diagonals  of  the  cross-section.  The 
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observations  of  Chueng  et.  al.  suggest  that  the  latter  may  be  the  pre¬ 
ferred  mode  of  failure  for  brittle  materials. 

A  final  comment  should  be  made  about  the  finding  that  for  a  given 
value  of  the  Biot  number,  on  cooling  the  maximum  value  of  tensile  stress  in 
the  square  rod  exceeds  that  for  a  circular  rod.  This  implies  that  for  a 
given  quenching  medium,  square  rods  will  require  a  lower  critical  temperature 
difference  (ATc)  for  the  initiation  of  thermal  stress  failure.  Firstly, 
this  suggests  that  for  a  given  quenching  environment,  the  use  of  square 
rod  specimens  will  permit  comparative  testing  of  a  wider  range  of  materials. 
Secondly,  the  use  of  square  rod  specimens  with  their  lower  values  of  AT^ 
will  suppress  changes  in  the  heat  transfer  coefficient  due  to  nucleate  boil¬ 
ing  or  vapor  formation  in  quenching  media  with  high  vapor  pressures.  This 
will  facilitate  the  Interpretation  of  data  for  the  comparative  thermal 
shock  resistance  of  materials  with  large  differences  in  ATc> 
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Maximum  tensile  axial  (oz) ,  transverse  (oy)  and  principal  stress  (op),  as  a  function  of  the 
Biot  number  in  a  long  cylinder  with  square  cross-section  subjected  to  sudden  convective  heat 
ing  (H)  or  cooling  (C) . 


Figure  3b.  Temperatures  at  various  positions  in  a  long  cylinder  with  square  cross-section  at  the  time 
of  maximum  tensile  axial,  az  (o,o)  and  transverse  oy  (o,o)  stresses  for  heating  (H) ,  and  at 
the  time  of  maximum  axial  tensile  stress  oy  (a,o)  for  cooling  (C)  as  a  function  of  Biot  num 
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ABSTRACT 


High  densities  of  micro-cracks  can  have  a  profound  effect  on  the  conduction 
of  heat  through  materials.  An  analysis  of  this  effect  on  thermal  conductivity 
is  presented.  The  variables  which  affect  the  formation  of  micro-cracks  in 
brittle  materials  due  to  thermal  expansion  mismatches  are  discussed.  It  is 
shown  that  micro-cracking  results  in  a  "thermo-mechanicr lly"  coupled  thermal 
conductivity  which  is  a  function  of  elastic  properties,  the  coefficient  of 
thermal  expansion,  fracture  toughness,  microstructural  variables  such  as 
grain  size,  as  well  as  the  closure  and  healing  of  cracks  at  higher  temperatures. 

These  effects  are  demonstrated  by  experimental  data  for  the  temperature 
dependence  of  the  thermal  diffusivity  of  polycrystalline  iron  titanate,  magne¬ 
sium  dititanate,  composites  of  alumina  with  silicon  carbide,  and  a  glass  con¬ 
taining  a  dispersed  phase  of  nickel. 


INTRODUCTION 


Materials  with  high  densities  of  micro-cracks  are  of  current  or  potential 
interest  for  engineering  applications  involving  high  temperature.  Such  cracks 
can  arise  from  creep-deformation,  fatigue  or  mechanical  failure  due  to  thermal 
shock.  High  densities  of  micro-cracks  can  result  from  the  thermal  expansion 
anisotropy  of  the  individual  grains  of  polycrystalline  aggregates.  Mismatches 
in  the  coefficient  of  thermal  expansion  of  the  component  materials  in  a  composite 
also  represents  an  important  source  of  micro-crack  formation. 

High  densities  of  cracks  are  expected  to  have  a  large  effect  on  material 
properties.  A  quantitative  understanding  of  such  effects  is  essential  for 
purposes  of  material  selection,  prediction  of  performance  and  failure  analysis. 
Previous  experimental  and  analytical  studies  have  shown  that  extensive  micro¬ 
cracking  can  promote  significant  increases  in  thermal  shock  resistance [ 1 . 2  1 
with  corresponding  large  decreases  in  the  moduli  of  elasticity  and  Poisson's 
ratio[3,4]  and  large  increases  in  Fracture  toughness [ 5 ) .  Of  significance  also 
is  that  micro-cracking  in  brittle  materials  changes  the  usual  highly  catastrophic 
(unstable)  failure  mode  to  a  stable,  far  more  controllable,  mode  of  failure[2]. 
This  latter  effect  is  very  important  for  materials  of  construction  of  many  high- 
temperature  structures  or  components,  which  frequently  are  expected  to  provide 
satisfactory  performance  even  while  in  a  fractured  condition.  For  these  various 
reasons,  micro-cracking  constitutes  a  very  effective  mechanism  for  obtaining 
significant  improvements  in  the  engineering  properties  of  brittle  materials. 


The  thermal  conductivity  of  materials  also  is  expected  to  be  affected 
strongly  by  the  presence  of  micro-cracks  especially  at  relatively  low  tem¬ 
peratures  at  which  heat  transfer  by  radiation  across  the  cracks  is  minimal. 

The  purpose  of  this  paper  is  to  review  previous  analytical  results  for 
the  effect  of  micro-cracking  on  thermal  conductivity.  The  thermo-mechanical 
properties  which  affect  micro-cracking  to  result  in  a  "thermo-mechanically" 
coupled  thermal  conductivity  will  be  discussed.  Experimental  data  which  illus¬ 
trate  this  effect  will  be  presented. 


ANALYSIS 

A.  Effect  of  Cracks  on  Thermal  Conductivity 

Solutions  for  the  effect  of  cracks  on  thermal  conductivity  were  presented 
by  Willis[6]  for  crack  oriented  perpendicularly  to  the  direction  of  heat  flow. 
The  present  writer[7]  analyzed  the  effect  on  thermal  conductivity  of  cracks 
oriented  randomly,  perpendicularly  and  parallel  to  the  direction  of  heat  flow. 
These  latter  analyses  were  based  on  the  theoretical  work  of  Fricke[8]  and  Powers 
[9]  on  the  thermal  conductivity  of  composites  consisting  of  a  continuous  matrix 
with  dilute  dispersions  in  the  shape  of  oblate  ellipsoids  of  revolution  with 
major  axis  b  and  minor  axis  a.  By  setting  the  thermal  conductivity  and  minor 
axis  equal  to  zero,  solutions  are  obtained  for  the  effect  of  infinitely  thin 
pores  in  the  form  of  circular  (penny-shaped)  cracks  on  thermal  conductivity, 
as  follows: 


a.  Randomly  oriented  cracks 

For  the  solution  of  Fricke[8]  and  Powers[9]  the  thermal  conductivity  (k) 
of  a  matrix  with  randomly  oriented  circular  cracks  of  zero  thickness  and  zero 
thermal  conductivity  can  be  written: 

K  =  K  (l+2bV/3ira)-1  (1) 

o 

where  K0  is  the  thermal  conductivity  of  the  matrix  without  cracks  and  V  is  the 
volume  fraction  occupied  by  the  cracks.  However,  for  infinitely  thin  cracks, 
a  =  0  and  V  =  0,  which  on  substitution,  renders  eq.  1  indeterminate.  Fortu¬ 
nately,  this  indeterminacy  can  be  removed  by  noting  that  for  an  oblate  ellipsoid 
the  volume  equals  4irab3/3.  For  N  cracks  per  unit  volume,  the  total  volume  V 
occupied  by  the  cracks  equals: 

V  =  47rNab2/3  (2) 

which  upon  substitution  into  eq.  1  yields: 

K  =  K  (l+8Nb3/9)_1 

o 


b.  Cracks  oriented  perpendicular  to  the  direction  of  heat  flow. 


Following  the  same  approach  as  for  randomly  oriented  cracks,  the  solutions 
of  Fricke  and  Powers  with  the  aid  of  eq.  2  for  the  effect  of  cracks  oriented 
perpendicular  to  the  heat  flow  on  thermal  conductivity,  becomes: 

K,  =  K  (l+8NbV3)  (4) 

i  0 


c.  Cracks  parallel  to  heat  flow. 

For  this  case  the  solutions  of  Fricke  and  Powers  indicate  that  the  cracks 
do  not  inhibit  heat  flow  such  that 


K|,  -Ko  (5) 

Comparison  of  eqs.  3,  4  and  5  indicates  that  cracks  oriented  perpendicularly 
to  the  direction  of  heat  flow  are  most  effective-  in  reducing  the  thermal  con¬ 
ductivity,  followed  by  randomly  oriented  cracks  in  turn  followed  by  cracks  parallel 
to  the  heat  flow,  which  have  no  effect  at  all.  It  is  also  of  interest  to  note 
that  the  dependence  of  cracks  on  thermal  conductivity  expressed  by  the  quantity 
Nb3  was  also  found  for  the  effect  of  cracks  on  elastic  behavior.  Furthermore, 
for  a  given  value  of  Nb3,  the  relative  effect  of  cracks  on  elastic  moduli  is 
about  twice  the  relative  effect  on  thermal  conductivity. 

Figure  1  shows  the  relative  thermal  conductivity  of  a  matrix  with  cracks 
of  various  orientations  as  a  function  of  Nb3. 


B.  Material  Variables  Which  Control  Micro-Crack  Formation 

To  provide  additional  background  for  the  experimental  data  to  be  presented, 
the  material  properties  and  other  variables  which  control  the  formation  and 
nature  of  micro-cracks  should  be  examined.  The  formation  of  micro-cracks  in 
brittle  materials  is  due  primarily  to  mismatches  in  the  coefficients  of  thermal 
expansion.  For  single-phase  polycrystalline  materials  such  a  mismatch  can  arise 
from  the  thermal  expansion  anisotropy  of  the  individual  grains.  For  composites, 
differences  in  the  coefficients  of  thermal  expansion  of  the  components  give 
rise  to  a  thermal  expansion  mismatch. 

Such  thermal  expansion  mismatches  result  in  the  generation  of  "internal" 
stresses,  whenever  the  temperature  is  changed  from  the  value  at  which  the 
differences  in  thermal  strains  are  zero.  This  latter  temperature  usually 
corresponds  to  the  temperature  at  which  the  polycrystalline  material  or  com¬ 
posite  was  manufactured. 

The  material  properties  which  affect  the  magnitude  of  the  "internal  stress" 
can  be  determined  from  the  solution  of  tiie  radial  stress  component  (or)  at  the 
interface  between  a  single  spherical  isotropic  inclusion  contained  within  an 
infinite  matrix,  derived  by  Selsing[10]: 

AaATF.d 

°r  =  0+^“T(EH/2F.  J+(L-2v'd')  (6) 

m  cl  m 


where  Aa  is  the  difference  in  the  coefficients  of  thermal  expansion,  E  is 
Young's  modulus,  v  is  Poisson's  ratio  and  the  subscripts  d  and  m  refer  to  the 
inclusion  and  matrix,  resp. 

Eq.  6  indicates  that  in  addition  to  the  thermal  expansion  mismatch,  the 
magnitude  of  internal  stress  is  a  function  of  the  elastic  properties  as  well. 

More  complex  solutions  for  the  internal  stresses  in  composite  circular  cy¬ 
linders  [11]  and  composite  spheres [12]  show  that  the  magnitude  of  the  internal 
stresses  is  a  function  of  the  volume  fractions  of  each  component  as  well. 

Very  important  to  note  is  that  all  literature  solutions  indicate  that  the 
magnitude  of  internal  stress  is  independent  of  the  dimensions  of  the  grains 
in  polycrystalline  materials  or  inclusions  in  composites. 

If  brittle  materials  failed  at  some  pre-determined  stress  level,  micro¬ 
cracks  could  form  whenever  the  magnitude  of  internal  stress  exceeds  the  value 
of  strength.  The  latter  value  usually  ( orresponds  to  the  tensile  strength, 
since  brittle  materials  are  most  prone  to  failure  in  tension.  Generally,  how¬ 
ever,  fracture  in  brittle  materials  does  not  occur  at  a  predetermined  stress 
level,  but  depends  on  the  size  and  geometry  of  defects  or  flaws  at  which  failure 
is  initiated.  The  same  should  hold  for  failure  resulting  in  micro-cracks. 

The  variables  and  their  relationships  which  control  micro-crack  formation 
should  be  examined.  An  isotropic  inclusion  will  be  considered,  which  is  placed 
in  a  state  of  uniform  triaxial  tension  due  to  a  mismatch  in  thermal  expansion 
with  a  surrounding  rigid  matrix.  The  inclusion  contains  a  precursor  flaw  in 
the  form  of  a  small  circular  (penny-shaped)  crack.  By  a  slight  modification  of 
the  analytical  theory  for  crack  stability  under  conditions  of  thermal  stress[2] 
the  critical  temperature  difference  ( ATC )  required  for  instability  of  the  pre¬ 
cursor  flaw  can  be  shown  to  be: 

AT  =  {nGd-.’v  )2/2E.(Aa)2(l-v  )r.}1/2  (7) 

c  mn  m 

where  G  is  the  surface  fracture  energy  required  to  create  unit  area  of  new 
crack  surface  and  l  is  the  radius  of  the  precursor  flaw. 

Eq.  7  indicates  that  if  for  a  precursor  flaw  of  given  size  AT  does  not 
exceed  ATC,  micro-cracks  will  not  form.  Conversely,  for  a  given  value  of  AT, 
the  size  of  the  precursor  flaw  must  exceed  some  minimum  size  before  micro¬ 
cracking  will  occur. 

The  defects  thought  to  be  primarily  responsible  for  micro-crack  formation 
in  brittle  materials  are  residual  pores  1  >cated  at  grain  boundaries  or  at  triple- 
points[13].  Such  pores  generally  have  a  flat  shape  and  can  act  as  effective 
cracks.  The  size  of  such  pores  general lv  is  proportional  to  the  grain  size. 

For  this  reason  it  is  expected  that,  although  the  magnitude  of  internal  stress 
is  independent  of  grain  size,  micro-crack  formation  is  expected  to  be  related  to 
the  grain  size  in  view  of  the  inverse  square  root  dependence  of  AT^  on  pre¬ 
cursor  crack  size.  For  this  reason,  the  effect  of  micro-cracking  on  thermal 
conductivity  is  expected  to  be  grain-size  dependent. 

Once  the  condition  for  micro-crack  formation  as  expressed  by  eq.  7  is 
satisfied,  strain-relief  and  eventual  crack  arrest  will  fol low.  On  reheating 
the  material  (i.e.,  reducing  AT)  two  effects  will  occur.  First,  depending  on 


crack  geometry  crack  closure  can  take  place.  Such  closure  will  provide  direct 
contact  across  the  crack  surfaces,  which  should  result  in  a  positive  dependence 
of  thermal  conductivity  on  temperature.  Secondly,  at  sufficiently  high  tem¬ 
peratures,  "crack-healing"  can  occur  by  diffusional  or  other  transport  processes. 
This  latter  effect  will  reduce  the  crack  size  and  lead  to  the  partial  or  com¬ 
plete  recovery  of  the  thermal  conductivity  to  its  original  value  before  micro- 
crack  formation.  For  this  reason,  the  thermal  conductivity  of  micro-cracked 
materials  is  expected  to  be  a  function  of  thermal  history. 

A  number  of  the  above  effects  are  demonstrated  by  experimental  data 
described  below. 


EXPERIMENTAL 


A.  Materials 

Experimental  results  for  the  effect  of  micro-cracking  on  thermal  transport 
behavior  are  reported  for  four  different  materials,  namely  polycrystalline  iron 
titanate[14] ,  polycrystalline  magnesium  dititanate[ 15] ,  composites  of  beryllium 
oxide  and  silicon  carbide[12]  and  a  glass  containing  a  dispersed  phase  of 
spherical  nickel  inclusions [16 ] .  In  the  iron  titanate  and  magnesium  titanate 
micro-cracking  results  from  the  thermal  expansion  anisotropy  of  the  individual 
grains.  In  the  BeO-SiC  composites,  micro-cracking  is  due  to  the  mismatch  between 
the  coefficients  of  thermal  expansion  of  the  BeO  and  SiC,  with  crack  formation 
occurring  in  the  component  with  the  higher  coefficient  of  thermal  expansion, 
i.e.  the  BeO.  In  the  glass-nickel  composites  micro-cracking  takes  the  form  of 
de-cohesion  at  the  glass-Ni  interface,  due  to  poor  mechanical  bonding.  All 
materials,  which  were  nearly  fully  dense,  were  made  by  hot-pressing.  For  the 
details  of  their  synthesis  and  micro-structures,  the  reader  is  referred  to 
the  original  papers  (12,14,15,16). 

B.  Experimental  Method. 

The  effect  of  micro-cracks  on  thermal  transport  was  determined  experi¬ 
mentally  by  measurement  of  the  thermal  difl’usivity  by  the  laser-flash  technique 
[17].  Details  of  the  equipment  are  described  in  the  original  studies  [12,14,15, 
16].  The  specimen  temperature  was  monitored  by  either  an  IR-detector  or  by  a 
thermocouple . 

Micro-cracking  has  little  or  no  effect  on  density.  Furthermore  the  specific 
heat  is  uneffected  by  the  presence  of  micro-cracks.  For  this  reason,  the  rela¬ 
tive  effect  of  micro-cracking  on  thermal  diffusivity  at  a  given  temperature  is 
expected  to  be  nearly  identical  to  the  relative  effect  of  micro-cracking  on 
thermal  conductivity. 

C.  Experimental  Results  and  Discussion 
1.  Iron  titanate 

Figure  2  shows  the  thermal  diffusivity  of  iron  titanate  as  a  function  of 
temperature  for  a  complete  heating  and  cooling  cycle  for  three  thermal  anneals. 


These  annealing  treatments  promoted  grain  growth  to  result  in  a  bi-modal  grain 
size  distribution,  the  fraction  of  large  grains  increasing  with  increasing  time 
of  anneal. 

The  experimental  data  for  the  thermal  diffusivity  indicate  a  number  of 
effects.  Firstly,  the  values  of  the  thermal  diffusivity  at  a  given  temperature 
decrease  with  increasing  time  of  anneal,  i.e.,  increasing  fraction  of  large  grains. 
Secondly,  the  negative  temperature  dependence  of  the  thermal  diffusivity  on  the 
heating  part  of  the  cycle  decreases  with  increasing  fraction  of  large  grains. 
Thirdly,  the  two  materials  subjected  to  grain-growth  show  a  pronounced  hysteresis 
on  heating  and  cooling. 

These  observations  are  in  accordance  with  the  expected  effect  of  the  micro- 
structural  variables  on  micro-cracking  and  thermal  transport  behavior,  discussed 
earlier.  The  fine-grained  (non-annealed)  material  shows  the  typical  temperature 
dependence  for  a  dielectric  material  in  which  heat  conduction  occurs  primarily 
by  phonon  conduction.  The  reproducibility  of  the  experimental  data  on  heating 
and  cooling  also  is  typical  for  dielectric  materials  which  do  not  undergo  a 
microstructural  change  on  thermal  cycling.  In  contrast,  the  thermal  diffusivity 
of  the  sample  which  underwent  grain-growth  is  less  than  the  value  for  the  fine¬ 
grained  material  by  as  much  as  a  factor  of  three.  This  effect  can  be  attributed 
to  the  micro-cracking  which  occurred  as  the  result  of  the  grain-growth.  Clearly, 
in  the  fine-grained  material  the  size  of  the  precursor  flaws  was  too  small  for 
micro-cracks  to  form.  In  the  coarse-grained  materials  the  number  and  size  of 
the  micro-cracks  was  difficult  to  ascertain  quantitatively.  However,  if  it  is 
assumed  that  the  micro-cracks  were  of  uniform  size  and  randomly  oriented,  the 
quantity  of  Nb3  in  eq.  3  for  the  coarse-grained  material  is  of  the  order  of  about 
two. 

The  more  positive  temperature  dependence  of  the  coarse-grained  material, 
compared  to  the  fine-grained  material,  is  the  direct  result  of  crack  closure  with 
increasing  temperature.  Possibly  the  effect  may  also  be  attributed  in  part  to 
in  increased  contribution  of  radiation  across  the  crack  to  the  heat  conduction. 

Of  particular  interest  is  the  counter-clockwise  hysteresis  observed  on  cool¬ 
ing  after  initial  heating.  This  behavior  is  indicative  of  permanent  crack  closure 
due  to  recombination  of  atomic  bonds  across  the  crack  surface  or  due  to  diffu- 
sional  processes.  After  such  crack  healing,  the  material  must  be  cooled  over 
some  range  of  temperature  before  the  micro-crack  can  form  again,  automatically 
leading  to  a  hysteresis  effect.  At  least  qualitatively,  all  the  above  pheno¬ 
mena  can  be  explained  in  terms  of  the  variables  which  control  micro-cracking 
and  in  terras  of  the  effect  of  temperature  on  the  size  of  the  micro-cracks  and 
thermal  diffusivity. 

2.  Magnesium  dititanate 

The  data  for  the  temperature  dependence  of  the  thermal  diffusivity  for  mag¬ 
nesium  dititanate  with  two  different  grain  sizes  are  shown  in  Fig.  3.  For  the 
finer-grained  material,  the  hysteresis  on  heating  and  cooling  is  clock-wise 
in  contrast  to  the  counter-clockwise  behavior  observed  for  the  iron-titanate. 

The  coarse-grained  magnesium  dititanate  exhibits  a  combination  of  clockwise 
as  well  as  counter-clockwise  behavior.  Comparing  these  results  with  those  in 
Fig.  2,  it  can  be  concluded  that  the  mechanism  of  micro-crack  formation  in 
magnesium  dititanate  differs  considerably  from  that  in  iron  titanate. 


It  is  suggested  here  that  micro-crack  formation  in  the  magnesium  diti- 
tanate  involves  substantial  subcritical  crack  growth  in  contrast  to  iron 
titanate  which  involves  catastrophic  crack  propagation.  Crack  extension  by 
sub-critical  crack-growth  requires  time.  For  this  reason,  it  is  conceivable 
that  during  the  time  period  of  cooling  the  specimen  from  the  hot-pressing 
temperature  the  total  extent  of  crack  growth  was  relatively  significant. 

For  this  reason,  on  reheating  the  thermal  diffusivity  will  exhibit  the  typical 
temperature  dependence  of  a  dielectric  material.  On  reaching  sufficiently 
high  temperatures,  however,  additional  crack  growth  can  occur  leading  to  a  de¬ 
crease  in  the  thermal  diffusivity  and  a  clock-wise  hysteresis.  On  the  other 
hand,  substantial  micro-crack  formation  by  catastrophic  crack  propagation  must 
have  occurred  during  the  initial  cooling  from  the  hot-pressing  temperature,  as 
indicated  by  the  lower  value  for  the  thermal  diffusivity  at  the  beginning  of  the 
thermal  cycle  in  comparison  to  the  fine-grained  material.  Such  cracks  are  ex¬ 
pected  to  close  or  heal,  leading  to  the  counter-clockwise  hysteresis,  at  least 
at  the  higher  temperature  range  of  the  thermal  cycle.  The  decrease  in  thermal 
diffusivity  to  a  value  less  than  the  original  value  observed  during  heating  can 
be  attributed  to  additional  micro-crack  formation  to  yield  a  crack  density  in 
excess  of  the  original  value.  For  this  reason,  the  clockwise  and  counter¬ 
clockwise  hysteresis  appears  to  be  the  result  of  the  competitive  effects  bet¬ 
ween  crack  closure  and  healing  and  additional  crack  formation.  This  explanation, 
if  correct,  suggests  that  in  assessing  the  effect  of  micro-cracks  on  thermal 
conductivity  information  on  the  kinetics  of  subcritical  c rack  growth  may  be 
required  as  well. 

3.  Beryllium  oxide-silicon  carbide 

The  thermal  diffusivity  of  these  composites  with  a  range  of  values  for 
the  volume  fraction  silicon  carbide  was  measured  over  the  temperature  range 
from  approximately  300°C  to  1400°C,  as  shown  in  Fig.  4.  Figure  5  shows  the 
calculated  thermal  conductivity  at  400°C  as  i  function  of  silicon  carbide  content. 
Included  in  the  figure  is  the  thermal  conductivity  predicted  by  the  Rayleigh- 
Maxwell  theory  for  the  crack-free  composites.  The  experimental  data  show  a 
large  decrease  near  20  and  80  percent  silicon  carbide,  whereas  near  40%  SiC, 
the  value  of  the  thermal  diffusivity  is  near  the  crack-free  material.  For 
an  interpretation  of  these  results  it  should  be  noted  that  detailed  calculations 
have  shown  that  the  magnitude  of  the  internal  stress  in  the  BeO  phase  increases 
linearly  with  increasing  volume  fraction  SiC. [12]  For  this  reason,  the  dependence 
of  the  thermal  diffusivity  on  SiC  content  cannot  be  explained  by  corresponding 
variations  in  the  magnitude  of  the  internal  stress.  A  more  likely  explanation 
is  based  on  the  changes  in  stress  distribution  with  increasing  silicon  carbide 
content.  At  the  low  and  high  concentraitons  of  SiC,  beryllium  oxide  is  the 
continuous  and  dispersed  phase  respectively.  However,  near  50%  SiC  both 
phases  are  expected  to  be  continuous.  For  one  of  the  phases  to  be  continuous 
the  internal  stress  distribution  will  approximate  that  in  a  composite  sphere, 
with  a  multi-axial  stress  state.  For  both  phases  to  be  continuous  near  50%  SiC. 
the  stresses  are  more  closely  described  by  a  composite  cylinder.  For  such  a 
model  the  stresses  in  the  BeO  are  uniaxial.  For  a  precursor  flaw  of  given  size, 
micro-crack  formation  is  more  likely  to  occur  under  a  multi-axial  stress  state 
for  which  the  orientation  of  the  flaw  is  less  important,  than  for  a  uniaxial 
stress  state  under  which  failure  will  occur  only  for  a  flaw  oriented  perpen¬ 
dicular  (or  nearly  so)  to  the  stress.  If  this  explanation  for  the  data  in 


Fig.  4  is  indeed  proven  to  be  correct,  it  indicates  that  estimates  of  micro¬ 
crack  formation  and  its  effect  on  heat  transport  behavior,  in  addition  to 
the  variables  discussed  above,  should  take  the  statistical  nature  of  brittle 
fracture  into  account  as  well.  Also,  this  explanation  illustrates  that  the 
existence  of  an  internal  stress  of  given  magnitude  is  a  necessary  but  not 
sufficient  condition  of  micro-crack  formation. 


4.  Glass-Ni  Composites 

These  composites  were  made  by  vacuum  hot-pressing  at  700 °C  glass  powders 
and  Ni  spheres  with  a  diameter  of  approximately  50  pm.  The  coefficient  of 
thermal  expansion  of  the  Ni  (a  -  13.4  x  10-6K-*)  exceeded  the  value  of  the 
glass  (~7xl0~6K.-* )  by  nearly  a  factor  of  two.  As  a  result,  on  cooling  these 
composites  the  relative  shrinkage  of  the  Ni  is  larger  than  for  the  glass. 

Nickel  and  glass  generally  show  very  poor  mutual  adhesion.  As  a  result,  on 
cooling  a  gap  will  form  at  the  interface  between  the  glass  and  Ni.  In  effect, 
at  the  lower  temperatures,  the  nickel  will  occupy  a  cavity  in  the  glass  larger 
than  the  size  of  the  nickel  particles.  In  analogy  to  the  effect  of  cracks  on 
thermal  conductivity,  such  gaps  should  have  an  effect  on  the  transport  of 
heat  through  the  composite. 

Figure  6  shows  the  data  for  the  thermal  diffusivity  as  a  function  of 
temperature  for  the  glass  with  0,  15,  30  and  45  vol.%  Ni.  Included  in  the 
figure  is  the  thermal  diffusivity  calculated  by  means  of  the  Bruggeman  variable 
dispersion  relation,  which  assumes  no  resistance  to  heat  flow  at  the  inter- 
face[18].  The  large  negative  deviation  of  the  experimental  data  from  the  calcu¬ 
lated  is  indicative  of  the  effect  of  the  resistance  of  the  interfacial  barrier 
to  heat  flow.  The  positive  temperature  dependence  of  the  thermal  diffusivity  is 
expected  from  the  closure  of  the  gap  with  increasing  temperature. 

In  spite  of  the  existence  of  the  inter  facial  gap,  the  data  indicate  that 
the  Ni  still  contributes  to  the  conduction  of  heat  as  indicated  by  the  presence 
of  the  Curie  point  and  also  by  the  fact  that  the  thermal  diffusivity  increases 
with  increasing  Ni  content.  This  latter  observation  would  not  be  expected 
if  the  Ni  did  not  aid  conduct  heat.  Heat  conduction  through  the  nickel  presumably 
occurs  as  the  result  of  multiple  contact  between  the  Ni  particles  and  the  wall 
of  the  cavity.  As  a  final  observation  it  should  be  noted  that  for  the  glass 
with  45%  Ni,  at  the  highest  temperature,  the  measured  value  for  the  thermal 
diffusivity  exceeds  the  calculated  value.  This  effect  is  thought  to  be  the 
result  of  Ni  particle-to-particle  contact,  not  accounted  for  in  the  Bruggeman 
theory. 


DISCUSSION 

The  experimental  data  presented  above  qualitatively  illustrate  many  of  the 
expected  effects  of  microcracking  on  thermal  diffusivity.  Because  of  the  quan¬ 
titative  uncertainties  in  calculations  of  microcrack  formation,  it  appears  unlikely 
that  quantitatively  reliable  estimates  of  thermal  diffusivity  of  micro-cracked 
materials  can  be  made  at  this  time.  From  a  practical  point  of  view,  however,  the 
present  observations  suggest  that  micro-( ra< king  represents  an  effective  mechanism 
by  which  the  thermal  diffusivity  of  solid  materials  can  be  tailored  to  specific 
values . 

1 

1 


The  present  data  also  indicate  that  although  the  basic  understanding  of 
heat  transport  phenomena  appears  to  be  well  developed,  many  uncertainties  still 
exist  at  the  microstruc. tural  level  which  offer  excellent  opportunities  for  further 
theoretical  analysis. 
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Fig.  2.  Temperature  dependence  of  the  thermal  diffusivity 
of  Iron  tit. mate  for  three  thermal  anneals  to  pro 
mote  grain-growth.  a:  As-prepared;  h:  Annealed 
16  hrs;  c:  Annealed  32  hrs.  (After  Tiebenerk,  et 
a  1  [16]). 
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Temperature  dependence  of  thermal  diffusivity 
of  magnesium  dititanate  for  two  different  mean 
grain  sizes,  a:  1.0  m<b  and  b:  70  ym  (After 
Siebeneck  et  al  (15]. 
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Thermal  diffusivity  of  SiC-BeO  composites 
as  a  function  of  Si C  content  for  a  range 
of  temperatures  (After  Singh,  et  al  [12], 


Thermal  conductivity  of  SiC-BeO  composites 
as  a  function  of  SiC  content  at  400°C  (After 
Singh,  et  al  [12]. 
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Fig.  6.  Comparison  between  theory  and  experimental  data 
for  the  thermal  diffusivity  of  composites  of  a 
glass  with  various  volume  fractions  of  Ni  spheri 
cal  inclusions.  (After  Powell,  et  al.,  [16]. 
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ABSTRACT 

An  analysis  is  presented  of  the  thermal  stresses  resulting 
from  "thermal  trapping"  in  semi-absorbing  brittle  ceramic  materi¬ 
als  in  the  form  of  a  thin  flat  plate  subjected  to  intense  thermal 
radiation.  Solutions  for  the  thermal  stresses  are  presented  for 
symmetric  and  assymmetric  radiation  heating  and  convective  cooling 
with  limiting  values  of  the  heat  transfer  coefficient,  h  =  0  and 
“*•  For  h  =  °°,  the  stresses  are  identically  equal  to  zero  for 
values  of  the  optical  thickness  ya  of  0  and  °°,  and  reach  their 
maximum  value  at  ya  =  1.3  and  2.0  for  symmetric  and  assymmetric 
heating  respectively.  For  h  *  0,  the  magnitude  of  thermal  stress 
increases  with  increasing  optical  thickness. 

Expressions  are  derived  for  the  maximum  radiation  heat  flux 
to  which  the  plate  can  be  subjected  without  failure,  in  terms  of 
the  pertinent  material  properties.  These  properties  are  combined 
in  "figures-of-merit"  for  the  selection  of  material  with  the  opti¬ 
mum  thermal  stress  resistance,  appropriate  for  given  heating  and 
cooling  conditions.  These  figures-of-merit,  which  may  contain  as 
many  as  seven  material  properties,  indicate  that  high  thermal 
stress  resistance  requires  low  values  of  the  coefficient  of  ther¬ 
mal  expansion.  Young's  modulus,  Poisson's  ratio,  emissivity,  and 
the  absorption  coefficient  in  combination  with  high  values  of 
tensile  strength  and  thermal  conductivity. 
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I.  INTRODUCTION 


Non-uniform  thermal  expansions  in  materials  or  structures 
subjected  to  non-linear  temperature  distributions  result  in  ther¬ 
mal  stresses.  The  magnitude  of  these  thermal  stresses  can  be 
sufficiently  high  to  result  in  catastrophic  failure.  This  prob¬ 
lem  is  particularly  severe  for  brittle  materials  which  do  not 
permit  thermal  stress  relief  by  plastic  flow.  Thermal  stress 
analysis  is  (or  should  be)  an  essential  element  in  the  design  of 
engineering  structures  or  components  operating  under  non-uniform 
temperature  conditions. 

In  practice,  the  problem  of  thermal  stress  failure  can  be 
solved  by  a  structural  re-design.  An  alternative  solution  is  to 
select  a  material  with  the  optimum  combination  of  properties 
which  keep  the  magnitude  of  the  thermal  stresses  well  below  the 
failure  stress.  The  magnitude  of  thermal  stress  is  not  only  a 
function  of  the  coefficient  of  thermal  expansion,  but  also  depends 
on  the  elastic  properties,  the  thermal  conductivity  and  diffusiv- 
ity  and  many  other  properties  depending  on  the  nature  of  heat 
transfer,  performance  criteria  and  mechanisms  of  failure  [1]. 

The  choice  of  the  optimum  material  to  avoid  thermal  stress  failure 
for  a  given  thermal  environment  can  be  based  on  "thermal  stress 
resistance  parameters"  or  "figures-of -merit."  The  two  best  known 
thermal  stress  resistance  parameters  are  [1,2]: 


R  =  St(l  -  v)/aE 

and 

R'  =  St(l  -  v)k/aE, 

in  which  S  is  the  tensile  strength,  v  is  Poisson's  ratio,  k  is 
the  thermal  conductivity,  a  is  the  coefficient  of  thermal  expan¬ 
sion  and  E  is  Young's  modulus.  Materials  with  higher  values  of 
the  parameters  R  and  R'  will  have  higher  thermal  stress  resis¬ 
tance  than  those  with  lower  values  of  these  parameters.  A  review 
of  some  thirty  other  thermal  stress  resistance  parameters  was 
presented  some  time  ago.  The  tensile  strength,  St,  is  present  in 
R  and  R'  because  of  the  high  ratio  of  compressive  to  tensile 
strength  of  brittle  materials,  which  makes  them  most  susceptible 
to  failure  in  tension. 

Radiation  represents  an  important  mechanism  of  heat  transfer, 
especially  at  high  temperatures.  The  incidence  of  thermal  stress 
failure  under  conditions  of  radiative  heat  transfer  is  receiving 
increased  attention  lately  as  the  direct  result  of  the  importance 
of  solar  radiation  as  an  energy  source.  An  analysis  of  the  ther¬ 
mal  stress  failure  of  brittle  materials  opaque  to  incident  black- 


body  radiation  was  presented  some  time  ago  [3],  This  analysis 
was  extended  to  include  materials  opaque  above  a  given  wavelength 
of  incident  black-body  radiation  and  completely  transparent  below 
this  wavelength  [4]. 

The  absorption  properties  of  many  materials  are  such  that 
incident  radiation  transmitted  through  the  surface  is  absorbed 
throughout  the  thickness.  This  results  in  the  so-called  "thermal 
trapping"  effect,  in  which  the  temperature  in  the  interior  of  a 
material  can  greatly  exceed  the  temperature  of  the  ambient  envi¬ 
ronment.  An  analysis  of  the  thermal  stresses  which  result  from 
this  type  of  heat  transfer  was  performed  recently  [5]  for  a  flat 
plate  symmetrically  heated  by  normally  incident  radiation  and 
cooled  by  convection  with  heat  transfer  coefficients  h  =  0  and  °°. 

The  purpose  of  this  paper  is  to  report  these  latter  results, 
together  with  results  for  the  thermal  stresses  in  a  plate  assym- 
metrically  heated  by  radiation  and  cooled  by  convection  again  for 
h  values  of  0  and  °°.  In  this  paper,  special  attention  is  paid  to 
the  interaction  between  the  coefficient  of  thermal  expansion  and 
the  other  material  properties  which  control  the  magnitude  of  ther¬ 
mal  stress.  Only  the  principal  analytical  results  will  be  high¬ 
lighted.  For  the  mathematical  details  the  reader  is  referred  to 
the  original  report  [5,6], 


II.  ANALYSIS 

The  analysis  of  thermal  stresses  first  requires  obtaining 
the  solutions  for  the  temperature  distributions.  From  these, 
expressions  are  derived  for  the  thermal  stresses.  By  equating 
the  maximum  tensile  thermal  stress  to  the  tensile  strength,  an 
expression  can  be  obtained  for  the  maximum  incident  radiation 
heat  flux  to  which  the  material  can  be  subjected  without  causing 
failure. 

For  the  analysis,  the  following  simplifying  assumptions  were 
made.  First,  reflectivity  of  the  material  of  the  plate  is  suffi¬ 
ciently  low  that  the  effect  of  multiple  internal  reflections  can 
be  neglected.  Secondly,  the  plate  will  be  assumed  to  be  at  low 
enough  temperature  that  the  maximum  value  of  thermal  stress  is 
reached  before  the  plate  becomes  sufficiently  hot  that  re- 
emission  of  the  absorbed  radiation  must  be  taken  into  account. 
Indeed,  the  validity  of  this  latter  assumption  is  easily  demon¬ 
strated  with  numerical  examples  [5].  Finally,  it  was  also 
assumed  that  the  emissivity  and  absorptivity  of  the  plate  are 
independent  of  wavelength.  If  required,  the  spectral  dependence 
of  these  properties  can  easily  be  taken  into  account  in  the 
analysis. 
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For  the  various  types  of  radiative  heating  and  convective 
cooling  conditions  and  the  rates  of  internal  heat  generation,  the 
general  approach  used  for  the  solutions  and  the  expressions  for 
the  maximum  value  of  tensile  thermal  stresses  are  as  follows. 

A.  Symmetric  radiative  heating  and  convective  cooling 

For  normally  incident  symmetric  thermal  radiation,  the  rate 
of  Internal  heat  generation  (g'1')  within  the  plate  can  be  derived 
to  be: 

g,,,«2lieqoe  Ua  cosh(px)  (1) 

where  u  is  the  absorption  coefficient,  e  is  the  emissivity 
(e  ■  1  -  r,  where  r  is  the  reflectivity),  q  is  the  intensity  of 
the  incident  heat  flux  at  each  side  of  the  plate,  a  is  the  half¬ 
thickness  of  the  plate  and  x  is  the  through- the- thickness  coordi¬ 
nate  with  x  ■  0  at  the  center  of  the  plate. 

Derivation  of  the  transient  temperature  (T)  requires  solution 
of  the  differential  equation  [7] 

32T/3 x2  +  g'"(x)/k  -  (1/6) 3T/ 3t  (2) 

where  k  is  the  thermal  conductivity,  t  is  the  time,  and  6  is  the 
thermal  diffusivity. 

The  thermal  stresses  in  the  plane  of  the  plate  are  obtained 
by  substitution  of  the  solution  for  the  temperature  in  [8]: 

a  a 

V* '  '-T  +  /“*  +  Txdx|  <3) 

-a  -a 

For  the  two  cooling  conditions,  the  initial  and  boundary  condi¬ 
tions  used  for  the  solution  of  Eq.  (2)  and  the  maximum  tensile 
stresses  and  maximum  permissible  heat  flux  are  as  follows: 

1.  Heat  transfer  coefficient  h  =  0. 

Initial  and  boundary  conditions: 

T(x,0)  =  Tq;  3T/3x(0,t)  -  3T/3x(a,t)  =  0  (4) 

For  these  conditions  the  maximum  tensile  thermal  stresses 
which  occur  at  x  ■  0  and  t  ■  00  can  be  derived  to  be: 
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— ua 

aEe<J«e  ?  a  0 
Jy»z  ‘  {u  +  <3  ' 


(5) 


In  general,  by  equating  the  maximum  tensile  therm*!  stress 
to  the  tensile  strength,  Eq.  (5)  can  be  rearranged  into  an  expres¬ 
sion  for  the  maximum  permissible  heat  flux,  q  ,  such  that  thermal 
stress  failure  is  avoided.  For  the  limiting  cases  of  an  optically 
thick  (pa  «°)  and  optically  thin  (pa  <<  1)  plate,  simple  analyti¬ 
cal  expressions  for  can  be  obtained: 


q  =  6S  (1  -  v)k/aEea 
’max  t 


3  4 


q  *  180S  (1  -  v)k/7oEep  a 
’max  t 


(pa  “) 
(pa  <<  1) 


(6a) 

(6b) 


2.  Heat  transfer  coefficient  h  =  “. 

For  this  case  the  Initial  and  boundary  conditions  are: 


T(x,Q)  *  Tq;  3T/3x  (0,t)  =  0 


T(a,t)  -  T(-a,t)  *  T 


(7a) 

(7b) 


The  tensile  thermal  stresses  are  a  maximum  at  t  =  00  and 
x  =»  -a, a  and  can  be  derived  to  be: 


o  (a,“)  = 

y.z 


4aE(pa)£qQe 
(1  -  v)ka‘ 


-pa  » 

■x —  cosh  (pa)  V*  {(p2  +  X2)(aX2)}  1  (8) 
*  /  -  n  n 


n=0 


where  X  =  (n  +  1/2) w/a 
n 


In  analogy  to  Eq.  (6),  the  maximum  permissible  radiation 
heat  flux  qmax  for  the  optically  thick  and  optically  thin  plate 
can  be  derived  to  be: 


max 


nax 


S  (1  -  v)krr 


(pa  -*•  »>) 
(pa  <<  1) 


64ciEpea 

B.  Assymmetric  radiation  heating  and  convective  cooling 


(9a) 

(9b) 


The  thermal  stresses  will  be  given  specifically  for  a  plate 
subjected  to  incident  radiation  on  one  face  (x  *  -a)  and  convec¬ 
tive  cooling  on  the  opposite  face  (x  *  a). 
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For  this  case  the  rate  of  internal  heat  generation  is: 


g'"  =  q^pe 


-p(x  +  a) 


1.  Heat  transfer  coefficient  h  =  0. 

The  initial  and  boundary  conditions  are 


(10) 


T<x,0)  =  Tq;  3T/3x(-a,t)  =  3T/3x(a,t)  =  0  (11) 


The  solutions  for  the  stresses  indicate  that  the  tensile 
thermal  stresses  reach  their  maximum  value  at  t  ■  »,  but  that  the 
position  of  the  maximum  tensile  stress  can  occur  anywhere  within 
the  plate,  depending  on  the  value  of  optical  thickness.  The  com¬ 
plete  solution  for  the  thermal  stresses  at  t  =  ®  as  a  function  of 
the  coordinate  x  is: 


aEq  e 

, _ 9 - 

y,z  (1  -  v)k 

aEq 


e  2,Ja)  +  (x  +  a) 


(1 


fe~“(x  +  •>  U  ±  «>2., 

[  M  4a  ^ 

fV_  f.-.2"  -  1  +  *(1  -  .  1 

'  L  2a2.  3  J 


3aEq  ex 
o 


2(1  -  v)a  k  (p 
3 


^j[e_2ua(pa  +  1)  +  (ya  -  1) ] 


+  f  (i  -  e-2va)  -  is: 


(12) 


For  the  optically  very  thick  plate  the  maximum  tensile 
stresses  occur  at  x  =  0  and  can  be  written: 


aEq  ea 

o  (0,®)  *  ° — r -r 

y,z  12(1  -  v)k 


(pa  -*■  ®) 


(13a) 


For  the  optically  very  thin  plate  the  maximum  tensile  ther¬ 
mal  stress  occurs  at  x  *  -a,  and  is  given  by: 

aEq  epa^ 

■  o^jr  <“*  ~  *>  «3b> 

The  expressions  for  the  maximum  permissible  heat  flux  become 


q  *  12S  (1  -  v)k/aEea  (pa  ->  ®)  (14a) 

IudX  c 

=  St  (1  -  v)k/aEepa2  (pa  <<  1)  (14b) 
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2.  Heat  transfer  coefficient  h  =  °». 


The  initial  and  boundary  conditions  are: 


T(x,0)  =  T(a, t)  =  Tq;  3T/3x(-a,t)  =  0 

The  tensile  thermal  stresses  have  their  maximum  value  at 
x  ”  -a  and  t  =  °°,  and  are  given  by 


y.z 


16Vpa  p  +  (n,Ma)e~2tJa 
1  -  v  /  .  2  2,  2  .  ,  , ,  v  2 


aE 


CQ 

I 


n=  odd 


n  j  k  p  +  (mr/4a) 


aE  V*'  8Vy  p  +  (mr/4a)e  j 

1  -  V  L  _2„2,.  .  2  A  /a„\2  X„ 

00 

3aE  V-1 

a(l  -  v) 


n  -  1 


1) 


n=odd 
8q  ep 


n  tt  k  p  +  (nir/4a) 
P  +  (mr/4a)e  r 


n 

n  -1 


n=odd 


n^w^k  p^  +  (mr/4a)^ 


t(~1)  2  -  h 

n  X 

n 


(15) 


where  X^  =  mr/4a  with  n  =  1,  3,  5,  7  . 


For  the  optically  thick  and  optically  thin  plates  Eq.  (15) 
can  be  simplified  to: 


a  (-a,”)  =  0 

y.z 


(pa  -*■  °°) 


o  (-a,00) 
y.z 


0 . 2719aEqoepa 


(1  ~  v)k 

which  yield  for  the  maximum  incident  heat  flux: 


(pa  <<  1) 


max 


‘max 


3.677St(l  -  v)k 
-  - 

aEepa 


(pa  -*•  «) 
(pa  <<  1) 


(16a) 

(16b) 

(17a) 

(17b) 


III.  NUMERICAL  RESULTS  AND  DISCUSSION 

For  the  symmetrically  heated  plate,  Fig.  1  shows  the  maxi¬ 
mum  tensile  thermal  stress  as  a  function  of  the  optical  thickness, 
pa.  For  h  ■  ”,  the  stresses  are  identically  equal  to  zero  for 
pa  ■  0  and  <*>,  with  a  maximum  at  pa  =  1.34.  For  pa  =  0,  no  ther¬ 
mal  stress  will  arise,  because  no  heat  is  absorbed.  For  pa  *  ®, 
the  radiation  is  absorbed  in  the  immediate  surface  and  removed  by 
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the  convective  medium.  Again,  no  heating  of  the  plate  occurs.  At 
pa  =  1.34,  the  thermal  trapping  effect  is  a  maximum  with  the  cor¬ 
responding  highest  values  of  temperature  increase  and  thermal 
stresses  within  the  plate.  These  results  suggest  that  in  general 
for  a  highly  efficient  solar  collector  (pa  -*■  ») ,  the  heat  should 
be  removed  from  the  plate  at  the  same  side  as  the  incident  radia¬ 
tion.  For  h  =  0,  the  stresses  are  zero  for  pa  =  0  and  increase 
monotonically  with  increasing  value  of  pa. 

Figure  2  shows  the  maximum  values  of  tensile  stress  as  a 
function  of  the  optical  thickness  of  the  plate  heated  on  one  side 
by  the  incident  radiation  and  cooled  on  the  other  side  by  convec¬ 
tion.  Qualitatively,  Fig.  2  looks  similar  to  Fig.  1.  The  magni¬ 
tude  of  the  stresses  for  the  assymmetrically  heated  plate,  how¬ 
ever,  is  much  less  than  that  for  the  symmetrically  heated  plate. 
The  principal  reason  for  this  is  that  in  the  assymmetrically 
heated  flat  plate,  the  non-uniform  thermal  expansions  in  part  can 
be  accommodated  by  bending. 

For  h  =  00  and  pa  =  0,  no  stresses  occur  because  no  heat  is 
absorbed.  For  pa  =  °°  and  h  =  °°,  however,  although  the  heat  is 
transmitted  through  the  plate,  no  stresses  occur  because  at  t  =  “ 
the  temperature  distribution  is  linear,  which  results  in  zero 
thermal  stress.  For  h  =  0,  the  stress  is  zero  for  pa  =  0  and 
increases  monotonically  with  increasing  pa.  Included  in  Fig.  2 
are  the  values  of  the  maximum  transient  stresses  for  h  =  °°.  For 
these,  the  reader  is  referred  to  the  original  study  [6]  for 
further  details. 

Attention  can  now  be  focused  on  the  role  of  the  coefficient 
of  thermal  expansion  and  its  interaction  with  the  other  material 
properties  establishing  the  magnitude  of  the  thermal  stress  resis¬ 
tance  of  semi-absorbing  ceramics.  It  may  be  noted  that  the  vari¬ 
ous  expressions  for  qmax  are  a  function  of  a  numerical  constant, 
the  thickness  of  the  plate,  and  a  number  of  material  properties. 
For  a  given  heating  and  cooling  condition  and  plate  thickness, 
thermal  stress  resistance  is  governed  only  by  the  material  proper¬ 
ties.  For  the  purpose  of  selection  of  the  material  with  optimum 
thermal  stress  resistance,  so-called  "thermal-stress-resistance" 
parameters  or  figures-of-merit  can  be  separated  from  the  expres¬ 
sions  for  q  as  follows: 
hnax 

S  (1  -  v)k  S  (1  -  v)k  S  (1  -  v)k 

-E— s - ;  — - ~ — ;  -E— = — —  (18) 

aEep  „  3  aEc 

«Ecp 

Materials  with  higher  values  of  these  figures-of-merit  have  better 
thermal  stress  resistance  than  materials  with  lower  values.  For 
this  reason,  high  thermal  stress  resistance  is  associated  with 
high  values  of  tensile  strength  and  thermal  condvic t ivi ty ,  coupled 


Fig.  1  Magnitude  of  maximum  tensile  thermal  stresses  as  a  func¬ 
tion  of  optical  thickness  in  flat  plate  symmetrically  heated  by 
thermal  radiation  and  cooled  by  convection. 


Fig.  2  Magnitude  of  maximum  tensile  thermal  stresses  as  a  func¬ 
tion  of  optical  thickness  in  flat  plate  assymmetrically  heated  by 
thermal  radiation  on  one  side  and  cooled  by  convection  on  the 
opposite  side. 


with  low  values  of  the  coefficient  of  thermal  expansion.  Young's 
modulus,  emissivity  and  absorption  coefficient.  For  optically 
thin  materials  (pa  <<  1),  thermal  stress  resistance  is  inversely 
proportional  to  p  or  p3  depending  on  the  heating  and  cooling  con¬ 
dition,  but  is  independent  of  p  for  optically  thick  materials  and 
h  *  0.  Of  course,  since  the  objective  of  solar  collectors  is  to 
absorb  as  much  radiation  as  possible,  the  values  of  e  and  p  should 
be  non-zero.  For  this  reason,  trade-offs  in  design  and  material 
properties  may  be  required. 

In  terms  of  the  objective  of  the  present  conference,  it  is 
of  interest  to  note  that  in  establishing  the  values  for  the 
figures-of-merit  in  Eq.  (18),  the  coefficient  of  thermal  expan¬ 
sion  interacts  with  as  many  as  six  other  material  properties, 
including  optical,  thermal  as  well  as  mechanical  properties.  As 
far  as  the  present  writers  are  aware,  no  f igure-of-merit  for 
materials  in  such  engineering  fields  as  aerospace,  thermo¬ 
electricity,  etc.  contains  as  many  dissimilar  material  properties 
such  as  those  given  by  Eq.  (18).  Clearly,  the  coefficient  of 
thermal  expansion  will  play  a  significant  role  in  determining  the 
feasibility  of  high-intensity  solar  energy  collectors. 
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